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LXXXTIT. Crystal Dislocations.—Elementary Concepts and Definitions. 


By F. C. Franx, 
H. H. Wills Physical Laboratory, University of Bristol*, 


[Received April 23, 1951.} 


SUMMARY. 


The elementary concepts of crystal dislocation theory—in particular 
the Burgers vector, perfect and imperfect dislocations, sessile dislocations, 
extended dislocations, twin and twinning dislocations—are defined in 
a manner independent of particular crystal lattices or models or special 
orientations of the dislocation: and it is shown*how their essential 
properties arise, in so far as they are independent of the particular crystal 
in which they occur. 


f 


§1. INTRODUCTION. 


Durine the past three or four years, those accustomed to deal with 
the theory of crystal dislocations have to some extent developed a private 
language and a number of special concepts, each simple enough in itself 
but together forming a substantial self-contained body of theory. The 
definitions of these concepts are to be found scattered in various papers 
and have been propagated in conferences and summer schools. A formal 
collection seems desirable. The concepts and definitions have largely 
originated in discussions at Bristol, in which, among others, 
F. R. N. Nabarro, G. Wyllie and the writer have been involved. We 
attribute the basis of our approach to a paper by J. M. Burgers (1939). 
We have also made use of important contributions by W. Shockley. 


§2. THe IpHAL CRYSTAL. 


The ideal crystal is built on a mathematical lattice, such that ions of 
the pth kind are to be found at the points 


rp—la+mb+ne+ry q; 


where 1, m, n, are all positive or negative integers, including zero ; p takes 
the values 1, 2, . . . up to uw which is the number of kinds of ion in the 
ideal chemical formula of the substance; and q takes the values 
1,2... up to x, which is the number of ions of the pth kind per unit 
cell. The unit cell is the parallelipiped whose three contiguous edges 


* Communicated by the Author. 
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are the lattice vectors a, b,c. Any three independent linear combinations 
of a, b and c make alternative lattice vectors. When a, b and c are so 
chosen that the v, are as small as possible, they are said to be primitive 
lattice vectors, and the lattice (la-+-mb-+ne) i is called the Bravais lattice 
of the crystal. The choice is not unique, and all vectors which are 
members of a set of three satisfying this condition will be called primitive 
vectors of the lattice. It is these which determine the perfect dislocations 
(to be defined below) of which the lattice is capable. Classified according 
to symmetry, there are fourteen different Bravais lattices. For example, 
if a, b and c are equal in length and make equal angles with each other, 
the lattice is in general rhombohedral : if the angles are 60°, 90° or 1093° the 
lattice has cubic symmetry, these three being the face-centred cubic 
(f.c.c.), simple cubic (s.c.) and the body-centred cubic (b.c.c.) lattices. 
Sodium chloride and aluminium have the same (f.c.c.) Bravais lattice 
and are equivalent in respect of their perfect dislocations. 

The further classifications by symmetry which depend on the vectors 
Yy,q, and which give 32 crystal classes and 230 space-groups, do not 
appear to be useful in classifying the further possibilities for dislocations : 
the latter are associated more with the twinning possibilities of the lattices, 
which have never been brought into a formal framework in an entirely 
successful way. 

The conventional unit cell of a crystal is one displaying the symmetry 
properties, which is not necessarily the primitive cell. Thus, for the 
three members of the cubic system mentioned above, one uses cubic 
cells containing respectively 4, 1 and 2 lattice points of the Bravais 
lattice. 


§ 3. ActuAL CRYSTALS: GooD AND Bap REGIONS. 


The ions in actual crystals have thermal vibrations, and quantum- 
mechanical uncertainty of position which may be described as 
“zero-point vibrations”. Both of these may be eliminated by taking 
a time average. There may also be vacant lattice points, interstitial 
ions, and wrong ions substituted for the right ones at lattice points : 
all these atomic defects may be notionally rectified, unless present to 
an excessive extent (e.g., when there are so many interstitial ions that 
it becomes impossible to recognize which is the principal lattice and which 
the interstitial one). Finally, the crystal may be continuously deformed, 
e.g. elastically, and yet retain an unmistakable local correspondence 
to the ideal lattice, provided the strains are not too large. Following 
an admirably simple use of words suggested by Slater at the request of 
Shockley (private communication), we call all regions in which the 
strains, and density of atomic defects, are small enough to allow this 
unmistakable correspondence, regions of good ‘crystal : otherwise, regions 
of bad crystal. We shall sharpen the definition of good and bad crystal 
regions presently, 
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§4. ASSOCIATED Partus. 


Consider a good region of real crystal, with its vibrations and atomic 
defects notionally eliminated, and let a tetrahedron of four neighbouring 
ions ABCD in this region be associated, one for one, with four 
corresponding ions A’B’C’D’ of an ideal “ reference lattice’, forming 
an almost congruent tetrahedron. A fifth neighbouring ion, E, in the 
real crystal can then be associated unambiguously with one, E’, in the 
reference lattice, the angles AEB and A’E’B’, BEC and B’E'C’, etc., 
taken in corresponding pairs, being almost equal or, at the worst, more 
nearly equal (say, in the sense of least squares) than if any alternative 
ion of the reference lattice is taken for the associated ion E’. Continuing 
this process, using angles made with BCD and E to fix a sixth ion F, 
and so on, we may extend a path from ion to ion through any region 
of good crystal, with an unambiguously associated path in the reference 
lattice. Moreover, a path may be displaced to a neighbouring path 


Fig. 1. 


‘Moebius Crystals ” 


between the same end points, with an unambiguous corresponding 
displacement of the associated path. The whole region within which 
such displacements can be made will now be called “ good ”. 


§5. THe BurGEerRs CIRCUIT. 


A path which closes on itself in the real crystal will be called a Burgers 
circuit. A Burgers circuit which lies in a simply connected region of 
good crystal may be formed by expansion from the initial tetrahedron 
hence its associated path is also a closed circuit. The associated path 
of a Burgers circuit which links a topological torus (or infinite rod) of 
bad crystal does not necessarily close on itself. In this case the continuity 
of good crystal round the Burgers circuit ensures only that the 
associated path closes within a lattice vector. 


3K2 
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§6. ExcLuDED CASES. 


We shall extend the definition of bad crystal so as to avoid certain 
topological oddities, exemplified by fig. (1), which may be termed 
“Moebius crystals’ by analogy with “Moebius strips”. Such 
geometrical situations may possibly occur in pathological crystals like 
chrysotile. If unlimited strains are impossible, they necessarily have 
relatively large cores of emptiness, or of bad crystal in the sense already 
defined, and are of limited extension in at least one direction. A region 
extending from the core to a boundary is declared “ bad” so that the 
core cannot be encircled by a Burgers circuit. We do this wherever, 
otherwise, we should fail to have a single-valued correspondence between 
directions in the real crystal and the reference lattice. 3 


§7. CrystaL DISLOCATIONS. 


If the closed Burgers circuit in the real crystal has an associated path 
which fails to close on itself, the circuit is said to encircle one or more 
lines of crystal dislocation (or, briefly, dislocations). We shall call the 
vector required to complete the associated circuit in the reference lattice 
the resultant Burgers vector of the dislocations enclosed. Burgers 
(1939) gave the name “ dislocation strength’ to the modulus of this 
vector. 

When the Burgers circuit passes continuously through good crystal 
we shall say we have perfect dislocations. (Other cases will be considered 
later). Then the Burgers vector (regarded as a free vector) is independent 
of the starting point on the circuit, and invariant against displacements 
of the circuit through good crystal regions. It is necessarily a lattice 
vector of the Bravais lattice. 

If two Burgers circuits are separated only by good crystal, they may 
be united into one circuit, whose Burgers vector is the vector sum of the 
Burgers vectors of its component circuits. Conversely, one circuit may 
be dissected into two, the sum of whose Burgers vectors is that of the 
one, by making a path through good crystal connecting two points on 
its circumference. When the circuit cannot be further dissected in this 
way, it may be shrunk on to its central core of bad crystal, which contains 
a single dislocation. 

Since the Burgers circuit can be continuously displaced along the 
dislocation line, without change in its Burgers vector, a dislocation 
cannot terminate within the crystal. 


§8. Distocation Nopss. 


A dislocation line may branch. The Burgers circuit enclosing the one 
line may be displaced so that it encloses the two (or more) dislocation 
lines into which it divides : hence the sum of their Burgers vectors must 
equal that of the first line. A more symmetrical statement of the 
situation is that if the Burgers vectors of all dislocation lines meeting 
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at a node are defined by right-handed circuits, when looking outwards 
from the node, the sum of these Burgers vectors is zero. This corresponds 
to a vectorial version of Kirchhoff’s law. 


§9. Tae Kinematic Properties or PERFECT DISLOCATIONS. 


Dislocations can move in many types of crystal. Indeed, the principal 
role of the dislocation concept in the theory of plastic deformation is to 
describe briefly the main essentials of the collective atomic motions 
involved in crystal slip. The possible motions of dislocations include 
those which do and those which do not involve the production or 
disappearance of lattice vacancies or interstitial atoms. Consider 
a Burgers circuit described anywhere in the crystal and let a dislocation 
line move from outside the circuit to inside it, passing between two 
adjacent ions P and Q on the circuit. The closure residue of the 
associated path is thereby changed by the Burgers vector, b, of the 
dislocation. Hence the ion P is moved relative to Q by a distance, and 
in a direction, which locally corresponds to b in the reference lattice. 
If the real crystal is, locally, nearly perfect apart from the dislocation 
considered, we may make the orientation and spacing of the reference 
lattice correspond to the local region and say that the displacement 
of P relative to Q is b, to an approximation of the same order as is usually 
involved in small-strain elastic theory. In order that our statements 
may remain valid in severely distorted crystals, however, let us denote 
by b’ the local vector in the real crystal which corresponds to b which 
is an invariant vector associated with the reference lattice. We may call 
b’ the local Burgers vector when we need a name which distinguishes 
it from b. If the line element do of the dislocation moves so as to sweep 
a planar element of area, the unit vector normal to which is n, and there 
is an elastic stress represented by the stress tensor P so that the traction 
on the element of area is nP, the work done by the elastic field in this 
motion of the dislocation, per unit area swept, is (nP).b’. By virtue of 
the symmetry of the stress tensor this may also be written (b’P).n. 
Hither of these expressions therefore gives the magnitude of the force, 
due to an elastic field, per unit length of dislocation line, tending to move 
the dislocation in the specified plane. This corresponds simply to the 
existence of a force (b’P) do on each element of the dislocation line. 
Fuller derivations of these formulae have been given by Peach and 
Koehler (1950) and Nabarro (1951) (though n has a less general significance 
in Nabarro’s derivation). 

There is a change of volume (b’.n) per unit area when the dislocation 
moves in the plane specified by n. When this is positive, it must occur 
by production of lattice vacancies, or transfer of ions from interstitial 
sites to lattice sites: when it is negative it must occur by removal of 
lattice vacancies, or transfer of ions to interstitial sites. 

Such motions are accompanied by a change in free energy if the local 
concentration of atomic defects differs from its equilibrium value. In 
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the case, probably appropriate to metals in a suitable temperature range, 
that the only diffusing atomic defects we need consider are vacancies, 
of one kind of ion only, each occupying a lattice volume v’, and their 
concentration exceeds the equilibrium value by a factor «, this change in 
free energy per unit area swept in the plane whose normal is n is equal 
to (b’.n/v’)kT Ina. This may be expressed by saying that there is an 
osmotic force (b’ x do/v')kT In « on the dislocation. 

The change of volume is zero if, and only if (b’. nm) is zero: 7. €. if the 
dislocation line moves in a surface containing itself and its Burgers 
vector. This is called the glide surface of the dislocation. The only 
motions (if any) of a dislocation which can occur rapidly under small 
stresses are glide motions. 

In the special case that the dislocation line is parallel to its* Burgers 
vector, all motions are glide motions. A dislocation in this orientation 
is called a screw dislocation. 


Fig. 2. 


The glide surface of a dislocation. 


Disregarding local changes of b’ we may say that a dislocation is 
free to glide on the cylinder generated by lines parallel to its Burgers 
vector, passing through all points along an instantaneous position of the 
dislocation line, and also, whenever it becomes coincident with one of 
these generators, on any surface formed by moving the generator parallel 
to itself : so that its projection, along the direction of the Burgers vector, 
on to a fixed plane, is invariant in area and in shape, but may acquire any 
line extensions enclosing no area (fig. 2), 

Motion of a dislocation (other than a screw dislocation) normal to 
its glide surface, by removal or generation of vacancies (alternatively, 


by generation or removal of interstitialcies) will be called positive or 
negative climb. 
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A node of three or more dislocation lines cannot move by glide alone 
_ unless their glide surfaces have at least a line in common. In general, 
they have only a common point. For example, at a triple node the 
Burgers vectors necessarily have a common plane: unless all three 
dislocation lines also lie in this plane the node cannot move without 
climb. 

The motion of a dislocation node, changing the lengths of the dislocation 
lines which meet at it, can effect the progressive union of two or more 
dislocations into one, or the dissociation of one into a larger number. 


§10. EpGE anp Screw Distocations. 


The screw dislocation is one parallel to its Burgers vector, and is so 
called because lattice planes normal to this vector are transformed 
by the presence of the dislocation into one or more helicoidal surfaces. 
A dislocation lying normal to its Burgers vector is called an edge 
dislocation, since it may be thought of as lying at the edge of an extra 
plane inserted into the lattice. It provides a second case whose elastic 
field is relatively easy to calculate (Burgers 1939), and the elastic field 
of any straight dislocation may be found by superposing the fields of 
screw and edge dislocation, in proportion to the cosine and sine of the 
angle between the Burgers vector and the dislocation line. 


§11. TRANSLATION TWINS. 


It occurs commonly that two like or mirror-image crystals differing 
in orientation by a particular rotation or reflection will fit on to each other 
on a particular plane, crystallographically the same for each, with an 
energy very little greater than the continuing single crystal. This is 
ordinary twinning*. When two crystal lattices which are parallel, and 
differ only by a translation, will fit together on a particular plane with 
low energy, we shall call the situation “ translation-twinning”’. It is 
not recognized in classical crystallography because it cannot be detected 
with a goniometer ; nor is it visible by X-ray diffraction unless it repeats 
at short intervals in the specimen. The best-known examples of 
translation twinning are “stacking faults” in the sphere-packing 
lattices, i.e. breaches of the stacking rules which lead to face-centred 
cubic or hexagonal close-packing. But other examples should exist. 
For example, in the crystal of a long-chain paraffin it is probable that a 
displacement of one-half crystal parallel to the long axis of the molecules 
will lead to a number of positions at which the energy shows subsidiary 
minima, not much higher than the absolute minimum characteristic of 


the perfect crystal. 

ot a ae Le 
* The converse of this statement is not necessarily true. It probably is 

true for growth-twins and solid recrystallization twins : but the composition 

plane of some deformation twins.may not be a surface of particularly low 


energy. 
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We may distinguish two types of translation twin : type I, or intrinsic, 
in which the two different lattice patterns are followed right up to the 
twinning surface (irrespective of whether this occurs on a plane of atoms 
or no); and type II., or extrinsic, in which a layer is sandwiched in 
between, not belonging to the continuing patterns of either. 


§ 12. ImpERFECT DISLOCATIONS. 


Unlike an ordinary twin-boundary, a translation-twin-boundary 
need not go right through the crystal. When it does not, its edge in the 
interior of the crystal is a dislocation. Assuming that ; there is no 
appreciable singularity of elastic properties at the translation-twin surface, 
it has the same stress-field as an ordinary dislocation. It will be called 
an imperfect dislocation (alternatively, a partial dislocation) in contrast 
with perfect dislocations which are surrounded entirely by good crystal. 


Fig. 3. 


Burgers circuits for imperfect dislocations. 


The Burgers vector of the imperfect dislocation bounding a translation 
twinning surface of type I may be defined in just the same way as that 
of a perfect dislocation, except that the Burgers circuit must start and 
finish on the twinning surface instead of at a completely arbitrary point. 

The same prescription must be followed when the translation-twinning 
surface is of type I., with the additional rule that the interpolated layer 
must be crossed normally, taking an equal distance normal to the 
corresponding plane when repeating the circuit in the reference lattice. 

If the lattice allows translation-twin surfaces of more than one kind ; 
e.g. two kinds, A and B, we may have an imperfect dislocation at the 
line where an A surface joins a B surface. In this case we make one 
Burgers circuit round the edge of the A surface, and a second one enclosing 
both this and the AB junction (fig. 3) : the difference of the residues of the 
associated paths is the Burgers vector of the AB junction. 

The Burgers vector of an imperfect dislocation is not a lattice vector. 
The Burgers vectors of the various different dislocations which may 
bound a given translation-twin surface (or provide the junction between 
two) all differ by lattice vectors. - 


§ 13. Kinematics oF ImpERrEcT DisLocaTIons, 


Compared with perfect dislocations, imperfect dislocations are subject 
to an additional kinematic restraint. As the imperfect dislocation moves 
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it extends or contracts the bounding surface between crystal regions 
in translation-twin relationship. This boundary has low energy only when 
it is parallel to a particular crystallographic plane—the conformity plane. 
In any other orientation the boundary will have high energy, like an 
ordinary intercrystalline boundary. Estimates for the metals indicate 
that the stress required to make the dislocation move into a non- 
conformity plane will always exceed the practical ultimate breakage 
stress of the work-hardened metal. This motion might therefore occur 
on the microscopic scale in regions where stress amplifications due to 
other dislocations raise the stress to values near to the theoretical yield 
stress of the perfect crystal. For ordinary stresses, however, we may 
consider such motions forbidden. Motion which extends or contracts 
the translation twin surface in the conformity plane is resisted or assisted 
by a force per unit length equal to the surface energy of the translation- 
twin surface. The estimates which can be made indicate that this is a 
force of the same order of magnitude as that exerted by the highest 
ordinarily applicable stress, higher in some metals and lower in others. 

Briefly, under ordinary stresses an imperfect dislocation is confined 
to its conformity plane. If its Burgers vector does not lie in that plane, 
it cannot glide: if its Burgers vector does lie in that plane, it’ cannot 
climb. Frank (1949) proposed that imperfect dislocations for which 
glide is forbidden should be called “ sessile ”’. 


§14. ExtenpED DISLOCATIONS. 


Two imperfect dislocations forming the right and left boundaries 
of a strip of translation-twin surface may be the opposites of each other. 
In this case (except in a highly inhomogeneous stress field, or in a 
concentration of vacancies or interstitial atoms far from its equilibrium 
value) free energy will be released when they move together, annihilating 

the translation-twin surface and the dislocation pair. If they are not 
~ the opposites of each other they cannot annihilate, and if the scalar 
product of their Burgers vectors is positive (expressing the condition 
approximately) they repel each other with a force inversely proportional 
to their separation. This force is balanced against a constant force arising 
from the surface energy of the translation-twin surface. Hence they 
assume an equlibrium separation, forming an extended dislocation, as 
described for the face-centred cubic or close-packed hexagonal lattices 
by Heidenreich and Shockley (1947). The combination, being surrounded 
by good crystal, is necessarily a perfect dislocation. Conversely, a 
perfect dislocation, when brought into the conformity plane of a possible 
translation-twin, should spontaneously dissociate into a pair of imperfect 
dislocations with a ribbon of translation-twin surface between them, so 
lowering its energy. 

If the Burgers vectors of both imperfect dislocations lie in the conformity 
plane, the combination is free to glide in that plane ; but in no other, 
even when it or one of its components becomes a screw dislocation 
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And it cannot climb. Both statements apply for ordinary stresses, 
not exceeding the practical ultimate strength. 

We shall use the term extended dislocation, in general, for any stable 
set of imperfect dislocations, connected together by translation-twin 
surfaces. An extended dislocation is free to glide only if these surfaces, 
and the Burgers vectors of its components, are coplanar. 


§ 15. DisLocaTIONS IN INTERCRYSTALLINE BOUNDARIES. 


Intercrystalline boundaries in general constitute continuous “ bad ” 
regions separating the ‘‘ good” regions on either side. It is then 
impossible to make a Burgers circuit, runing continuously in good crystal, 
and encircling a portion of the boundary. When this is the case, we may 
prescribe that the boundary region should be crossed to the nearest ion 
in the good region on the other side, while in the associated path we make 
an equivalent displacement, and a further displacement to the nearest 
equivalent ion. This introduces an uncertainty of the order of magnitude 
of one lattice spacing into the resultant Burgers vector, at each crossing. 
These are generally negligible errors. However, to continue the path in 
the second good region we must re-establish the cardinal directions by 
choosing a further three ions, making a tetrahedron with the first, to 
associate with a tetrahedron of ions in the reference lattice. The number 
of ways in which this can be done depends on lattice symmetry—for 
example, in the cubic system it may be done in twenty-four different ways. 
For each of these choices we shall find a different value for the resultant 
Burgers vector of dislocations enclosed in the circuit. We can show 
(Frank 1950) that if r is the vector connecting the two points at which 
the circuit crosses the boundary, and the one crystal is brought parallel 
to the other by a rotation « about an axis | (a unit vector), this resultant 
Burgers vector is rXx1.2 sin }x. We have the same multiplicity as before 
in the choice of 1 and «, and can adopt the convention that we choose 
the 1 which makes « least. 

When « is small (say, less than 10° to 20°, according to the orientation 
of the boundary) the structure of the boundary which has lowest energy 
will probably be one consisting of discrete dislocation lines, a distance 
of the order of magnitude |b|/« apart, separated by “ good’’ regions 
(cf. Burgers (1939, 1940), van der Merwe (1950), Shockley and Read 
(1950) ). A portion of the boundary can then be encircled by a Burgers 
circuit made entirely in good crystal, and its dislocation content can be 
unambiguously defined. 

When « is large (it may take values up to 635° for a cubic crystal), the 
analysis of a grain boundary into line dislocations is of formal value only. 


§16. Twin anv Twinntne DisLocatrons. 
When the real crystal approximates to a perfect twin, the perfect 
twin may be taken as reference lattice ; the whole twinned real crystal, 
including the twin-boundary, being considered good. Ordinary perfect 
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and imperfect dislocations may be present, but there are additional 
possibilities. For, since the twin boundary may be on any one of an 
infinite succession of parallel equivalent lattice planes called composition 
planes, the reference lattice is multivalued. The possible Burgers vectors 
include all the vectors joining a point of the lattice to the points of its 
twin lattice interpenetrating it. Such dislocations correspond to steps 
in the twinning boundary, and are therefore confined to the boundary. 
They may be called twin dislocations. 

If the twinning law is one which can be represented by a shearing 
deformation, with all displacements parallel to the composition plane, 
the Burgers vector of the dislocation may be parallel to the composition 
plane (being then equal to the height of the step in the boundary, 
multiplied by the shear in twinning). If this is so the dislocation may 
glide under mechanical stress. This motion is the growth of a mechanical 
twin, and the dislocation involved will be called a twinning dislocation. 
Fig. (4) illustrates one (cf. Frank and van der Merwe, 1949). 


Fig. 4. 


Burgers circuit and its associated path for a twinning dislocation. 
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LXXXIV. The Axis of Distortion of a Twisted Elastic Prism 


By D. G. ASHWELL, 
Engineering Laboratory, Cambridge*. 


[Received April 11, 1951.] 


SUMMARY. 


Duncan, Ellis and Scruton (1933) have shown that if an elastic cantilever 
is twisted, the displacement of its tip relative to its support consists of a 
rotation about the flexural centre of its cross-section. In this paper is 
considered the distortion of a twisted prism without reference to its support, 
and it is found that the axis of distortion of such a prism is in general quite 
different from the axis of displacement found by Duncan, Ellis and Scruton— 
an effect that has been previously noticed by Cullimore (1949). General 
expressions are derived for the determination of the axis of distortion, and 
its effect on practical cantilevers subjected to finite twisting is estimated. 
Experiments designed to test the analysis are described, and their results 
are found to confirm it satisfactorily. 


§ 1. INTRODUCTION. 


THE stresses and strains given by the St. Venant theory for the torsion 
of a prismatic elastic cantilever, held at one end and twisted at the other, 
are independent of the axis about which the prism is assumed to twist. 
Duncan, Ellis and Scruton (1933), discussing this fact, remark that it 
would appear as though the position of the axis of twist were actually 
indeterminate—* a conclusion quite at variance with physical intuition ”’ 
—and go on to show that there can, in fact, be obtained from the St. 
Venant theory a position of the axis of twist for any specified relationship 
between the system of coordinate axes chosen and the root cross-section of 
the cantilever. This is because the root cross-section, like all other cross- 
sections, suffers warping out of its plane, and in order to specify the 
relationship between the cantilever and the frame of reference defined by 
the coordinate axes, three points in this cross-section must be regarded as 
fixed. The choice of these points determines the direction of the plane at 
the root of the cantilever relative to which its displacement is measured, 
and this determines the position of the axis of twist relative to the 
cantilever. In practical cases the direction of this plane will depend on 
the local stress distribution at the root—that is, on the system of local 
stresses which must be added to the St. Venant stresses in order to produce 
continuity between the cantilever root and its support. 

Thus the axis of twist is essentially a means of specifying the displacement 
of the cantilever as opposed to its distortion ; for while the displacement 
may be expected to depend on the frame of reference chosen, the distortion 


* Communicated by Professor J. F. Baker.’ 
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—that is, the displacements of various parts of the cantilever relative to 
each other—will be quite independent of the external frame of reference— 
local distortion at the root being, for the moment, ignored. Now physical 
intuition, while accepting the validity of the axis of twist as an axis of 
displacement, suggests that there must also be an axis of twist defining the 
distortion of the prism, and it is with the determination of this axis that the 
following four sections are concerned. The effect has been previously 
noticed by Cullimore (1949), who has given a tentative solution based on 
energy considerations. (See § 4, below.) 

The distinction between the axes of twist and distortion is made clear 
in §6, where is considered the effect of each of them on practical 
cantilevers. 


\ 


Figs: 


to) 


§ 2. DETERMINATION OF THE DISTORTION AXIS. 


First will be considered an elastic prism twisted uniformly in such a way 
that its distortion is the same everywhere along its length, and that there is 
no constraint against warping of its cross-sections anywhere. In general, 
longitudinal filaments that were straight before the prism was strained will 
become distorted into helices, and since the twisting is uniform, these 
helices may be expected to havea common axis. The filament lying along 
this axis will be the only one to remain straight, and its position in any 
cross-section will be called the distortion centre of the cross-section, and 
the axis itself will be called the distortion axis of the prism. ; 

Fig. 1 shows a prism being twisted. The centroidal axis has become 
distorted into the helix G’GG”, which has a radius R and makes an angle ¢ 
with the generators of the cylinder containing it. The torque necessary to 
maintain the twisting is T, and since its effect at all cross-sections is the 
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same, it must act about an axis parallel to the axis AA’ of the cylinder 
containing G’GG”—that is, about the distortion axis of the prism. GL is 
the generator passing through G of the cylinder containing G’GG", GD is 
a radius of the cylinder, and GJ is a line perpendicular to GD and to the 
tangent GK to G’GG", —_ 

The torque T may be resolved into components T cos ¢ twisting the 
prism about its centroidal axis, and T sin ¢ bending the prism about GJ. 

It can be shown from the geometry of the helix that the twist 7 per unit 
length of the prism about its centroidal axis and its curvature I/p about 


GJ are given by 


=F sin ¢ cos 
and 1 oi ten lity Biiye. 2 Oke weekowe (1) 
Ss rie 
pom d. 


Consider an elastic prism constrained to twist in such a way that its 
centroidal axis remains straight, and let its twist per unit length be constant 
and equal to 7. Let r be the distance of a longitudinal filament from the 
axis. Then the radius of the helix into which the filament distorts will be 
r, and its angle ¢,=tan-1(rr). Further, if e, is the strain of the filament 
lying along the axis, the filament distant r from the axis will suffer an 
additional strain of 1—cos ¢,, due to its distortion into a helix, giving a 
total strain of e)-+-1—cos ¢,. ‘Here it must be mentioned that the warping 
of cross-sections of the prism out of their planes, due to the twisting, will 
cause axial displacements of points on the filament, but that since 7 is 
constant and warping unrestricted, this displacement will also be constant 
along the length of any filament, and hence will produce no additional 
strain. 

The strain of the filament considered, then, is e),+1—cos ¢,, and if it 
may be assumed that lateral stresses acting upon it may be neglected 
(see § 3, below), the stress in it will be E(e)+1—cos ¢,), where E is Young’s 
Modulus for the material of the prism. The component of this stress in 
the direction of the centrodial axis will be E(e)-+-1—cos ¢,.) cos ¢,, and since 
o,=tan (rr), this may be expanded in terms of rz to give 


Eleo+4(r7)?—deq(r7)?—B(r7)* . 2... 


If, now, ey and (rz)? may be limited to the usual small order of elastic 
strains, the axial stress at a distance r from the centroidal axis may be 
taken as Ele9+-}(r7)*], terms containing higher powers of ey and (rz)? being 
neglected. The implications of this limitation are discussed in the next 
section. 

Since there is no resultant axial force in the prism, the value of ¢, may be 
found by integrating these stresses over the cross-section and equating 
their sum to zero. This gives 


Bej=py>=—fRe,) 2 re) 
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where Po is the stress corresponding to e, and k is the polar radius 
of gyration of the cross-section about its centroid. 

Let Ga and Gy be the principal axes of the cross-section (see fig. 2). 
Then if (2, y) is the position of an element of area SA referred to these 
axes, and r is the distance of 5A from G, the axial stresses will have 
resultant moments M, and M,, about Ga and Gy given by 

M,=tE7fr2y dA 
and M,=—4Er? {rx dA , 


where the integrals are taken over the whole cross-section. 


Fig. 2. 


Fa = My 


For convenience, let [fdA be written (f), so that, for example, 
fr?ydA=(r?y). Then 


M,=427(r7y) | (3) 


and — M,=— 3E7*(rz). 


M,, and M,, have the sense shown in fig. 2. 

Thus if the prism is to twist about its centroidal axis, in addition to the 
couple required to produce the twist, moments M,, and M,, must be applied. 
In the absence of such moments the prism will tend to take up curvatures 
—M,/EL, and —M,/EL,, where I, and I, are the moments of inertia of the 
cross-section about Ga and Gy. These curvatures will combine to produce 
a curvature about some other axis through G, and this, together with the 
twist 7, will cause the centroidal axis of the prism to distort into a helix, 
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Consider again the prism constrained to twist with its centroidal axis 
straight. The moments about Ga and Gy will be given by (3). Let it 
now be caused to distort so that its centroidal axis becomes a helix of 
radius R and angle 4. Let D be the position of the distortion centre in 
fig. 2 and GJ be parallel to the axis about which the resulting curvature I/p 
of the centroidal axis takes place; DGJ will be a right angle. Let 
JGzx be «. Then the curvatures of the centroidal axis about Gz and Gy 
will be (1/p) cos « and (1/p)sin« and these will require moments about 
these axes of (EI, cos «)/p and (EI, sin «)/p. Thus the total moments to 
be provided about Gx and Gy are 


EL, cos a 


EI, sina 
p 

If this distortion is to be that actually suffered by the prism 
when subjected to a pure torque T applied about the distortion axis, the 
components of T about Gx and Gy must equal these required moments. 
Now it was shown above that the component of T tending to bend the 
prism is T sin ¢ acting about GJ, and so its components about Ga and Gy 
are T sin d cos wand T sin ésin«. Thus the equations of equilibrium for 
the bending of the prism about Ga and Gy are 


T sin ¢ cos <=(EI,/p) cos «+M,, 
and . T sin ¢ sin «=(EI,/p) sin «+M,,. 


and “ML 


(4) 


Furthermore, the component of T twisting the prism about its centroidal 
axis is T cos ¢, and thus, if CJ is the torsional rigidity of the prism, 


T cos 6=7CJ. (5) 


Finally, 7 and p are related to R and ¢ by (1) above, and M,, and M 
are given by (3). Eliminating T, 7, p, M,, and M,, from (1), (3), (4) and (5), 
it is found that j 


> 


_ E(r?y) cos? f 
Th eR aac) 

and : (6) 
at ame K(r°a) cos? 6 


2R(EL,—GJ) ’ 


Squaring and adding these equations, the following equation relating R 
and ¢ is obtained : 


aS ed) eee ee 
EF? costg (HI,-GE (EL_yp 


An alternative form of this equation relates R and the angle of twist of 
the prism per unit length measured along the distortion axis. If 6 is this 
angle of twist, from the geometry of the helix 


COs = (14 RPP The soy ol ae ee (8) 
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and thus 
£2R2(L-P R202)2 6. (72)? rye 
ert 00)2. (RIC sa |) 
where 
(ry) = fry dA 
and (p= |rerdA. | (10) 


the integrations being performed over the whole cross-section. 
If Gx is an axis of symmetry of the cross-section, (ry) will be zero, and 
the equation relating R and @ becomes 


2R(1-+R*e2)(EL,—CJ)=E(r'x). =. 2... (11) 


It is worth pointing out that although equations (9) and (11) are cubic in 
R, they each have only one real root. If 6 is given, and R26? is not 
immoderately large, R can be found quickly to any desired degree of 
accuracy by the usual iterative method for such equations. The position 
of D can then be found using equations (6) and (8). If Ga is an axis of 
symmetry, D is the point (R, 0). 


§ 3. ASSUMPTIONS INVOLVED. 


In deriving these expressions two assumptions have been made. 

The first is that the stresses in longitudinal filaments are E times their 
strains. Now since the sum of the axial stresses must be zero, and since 
4(r7)? is always positive, é), the strain of the filament lying along the 
centroidal axis, will be negative—see equation (2). Thus, in general, the 
assumption is not valid, for if the cross-section has a compact shape with 
no internal boundaries, the tension in filaments near the external boundary 
will cause lateral contractions tending to squeeze the centre of the prism, 
which experiences longitudinal compression ; and so the ratio between 
stress and strain will depend on Poisson’s ratio as well as on Young’s 
Modulus. However, prisms with such cross-sections may be expected to 
have distortion centres closely coinciding with their centroids, and their 
behaviour in this field is not of much interest. In other types of cross- 
sections, such as open and closed hollow sections, there will be little or no 
resistance to lateral contraction, and the ratio between stress and strain 
may safely be taken as Young’s Modulus. 

The second assumption is that e) and (r7)? may be restricted to the order 
of magnitude of small elastic strains. This also need not always be true, 
for if the material of the prism is at a fairly constant distance from the 
centroidal axis, and if rz could be at all large, e) and (rz)? might have large 
numerical magnitudes of opposite sign without overstraining the material. 
But this could happen only in the case of a circular, or nearly circular, 
thin-walled hollow tube with a longitudinal slit—a special case in which, 
again, the distortion centre would very nearly coincide with the centroid. 

In all cases which are likely to be of interest, therefore, these two 
assumptions may be taken as valid. 
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§ 4, TyprcaL Posrrions oF THE DistoRTION CENTRE—AN APPROXIMATE 
EXPRESSION. 

In the case of most thin-walled open sections CJ may be neglected in 
comparison with EI and Ré will be small compared with unity. Equation 
(11) then becomes 

R==(r2z)/27,, 3S Sea ae hae 

This approximation will be true enough for many practical cases, and 
typical results obtained from it for thin-walled symmetrical sections are 
given in the Table. It should be noticed that the 4 given in the Table is 
not equal to R since it is not measured to the centroids of the sections. 

Equation (12) can be obtained very quickly if two modifications are made 
in the above analysis. Firstly, the bending component Tsin¢ must be 
ignored ; this has the effect of omitting the term CJ from (9) and (11). 
Secondly, the expression for the longitudinal stresses obtained in §2 
must be referred directly to the distortion axis instead of to the centroidal 
axis ; the effect of this is to invalidate the argument of the second part 
of §3, and to cause the omission of the factor (1+ R*0?) from (9) and 
(11). These effects cannot be ignored in the case of flexible sections 
twisted considerably, or if it is desired to extend the analysis to 
include the effects of axial forces. A tentative solution due to Cullimore 
(1949) uses an energy method which yields results for channel and angle 
sections identical with those obtained from (12) above. The approximate 
nature of this solution seems to arise because an exact solution cannot be 
obtained by treating the problem as a small-deflection problem. 


§ 5. GENERAL COMMENTS ON THE DISTORTION AXIS. 


Previous writers on this subject have ignored the effect of the 
longitudinal stresses here considered, and the curvature they produce, on 
the ground that they are second-order effects (see Duncan, Ellis and Seruton 
1933). But they seem to have overlooked the fact that from (1) and (8) 
it can be shown that the twist 7+ of the distorted centroidal axis is 
6(1+ R26?) and its curvature 1/p is R6?(1+ R262). Thus the distortion 
of the centroidal axis into a helix does in fact require that its curvature 
should be of the second order in + even if R is large, and so in determining 
the position of the axis of distortion second-order effects cannot be ignored. 

It must be emphasized that all the analysis in this paper so far is true 
only if all the cross-sections of the prism are allowed to warp out of their 
planes in the manner required by the St. Venant stress distribution. If 
one section is restrained against warping, the torque will be resisted near 
that section partly by the twisting of the prism and partly by the 
differential bending of its various parts. In this case the second-order 
longitudinal stresses due to twisting will certainly be negligible compared 
with the first-order stresses due to the differential bending near the con- 
strained cross-section, and the usual methods for estimating the additional 
stiffness of the prism due to the constraint will be unaffected by the work 
contained in this paper, . 
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§ 6. Frnrre DIsPpLACEMENTS OF PRACTICAL CANTILEVERS. 
It remains to distinguish between the axes of twist and distortion by 
considering their effects on the behaviour of practical cantilevers. 


TABLE. 


A (measured to mid-line) 
defines the position of 
the distortion centre. 


Section. 
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Consider a prismatic cantilever as shown in fig. 3. At the root its cross- 
section is everywhere attached to a rigid plane support. Its length is /, and 
SE 


828 D. G. Ashwell on the 


( is the flexural or shear centre at its tip, while D : the AiSvOrnOe nay 
there. The corresponding points at its root are C and D’. Let u e 
twisted by a torque about an axis perpendicular to the supporting plane. 
The stress distribution along the greater part of its length will be ey 
predicted by St. Venant, and its cross-sections will tend.to warp out ° 

their original planes. This they will be able to do except near the root, 
where an additional stress-system must be applied to produce continuity 
between the prism and the plane support. These continuity stresses 
applied at the root will have no resultant force or moment, and their effect 
may be expected to die out fairly quickly along the prism, 80 that at no 
creat distance from the root they will have become negligible. Neverthe- 
less, although they are confined to a short length of the cantilever, they may 


Fig. 3. 


Root PLANE 


affect its displacement significantly by producing a rotation of it about 
some axis in the fixed plane. Their precise nature is not known in general, 
and so it is not possible to determine exactly their effect except in special 
cases. Duncan, Ellis and Scruton, however, have shown that if the canti- 
lever is twisted slightly, the tip ABC will rotate about C, the centre of 
flexure. In other words, the continuity stresses are such that the flexural 
axis of the cantilever remains perpendicular to the root plane. 

[f the rotation of the tip becomes large, CC’ will not remain straight but, 
according to the analysis given earlier in this paper, will distort into a 
helix with DD’ as axis. This will be true for the greater length of the 
cantilever, except near the root, where the continuity stresses will 
determine the distortion. It will still be true, however, that near the root 
the filament lying along CC’ will be perpendicular to the root plane. 
Thus, since the length affected by the continuity stresses is small compared 
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with J, it can be said that the distortion of the cantilever will be that dis- 
cussed in previous sections, while its displacement will be such that the 
tangent to the flexural axis at the root will remain perpendicular to the 
root plane. This statement is not strictly accurate owing to the short 
length affected by the continuity stresses, but if / is sufficiently large for 
the rotation of the tip to be appreciable, the error involved may be 
expected to be small. 

It is now possible to estimate the displacement of the tip for large 
deflections. If the twist of the cantilever per unit length is 6, and the 
distance C’D’ is 6, it can be shown from the geometry of the helix that the 
rotation of DD’ relative to the tangent to CC’ at ©’ takes place about 
C’D’ and is of magnitude tan-1(65). In most cases likely to arise in 
practice this rotation will be small and may be safely approximated to 68. 
This means that since DD’ remains straight, D suffers a linear displacement 
at right angles to C’D’ of magnitude 615. Now the rotation of ABC, ;, say, 
will equal 6, and so the linear displacement of D will be #s. 

If DD’ becomes inclined to the root plane, the applied torque, T, 
assumed above to be about an axis perpendicular to that plane, will have 
a component tending to bend the cantilever as a beam, but this component, 
T66, will be small, and as prisms of types likely to be met with are stiff in 
bending, the distortion caused by this component will be ignored. It 
would not be a difficult matter to take account of it in any particular case. 

Thus if the quantity 06 is small, and / is large compared with the length 
near the root affected by the continuity stresses, the displacement of the 
tip ABC consists of (i) a rotation about D together with (ii) a linear 
displacement of D equal to 46 perpendicular to C’D’. The usually accepted 
centre of displacement, C, is now seen to be rather the instantaneous centre 
of rotation of the tip for the position J=0. 

To determine the finite displacement of the tip of a twisted cantilever it 
is thus necessary first to determine the positions of the flexural and distor- 
tion centres for the value of @ concerned, and then to superimpose the linear 
and rotational displacements stated above. 


§ 7. EXPERIMENTAL RESULTS. 


A series of experiments was carried out to test the analysis given in this 
paper. They consisted of applying a pure torque to cantilevers with 
different cross-sections and determining the displacements of their tips for 
different angles of rotation. 

Three specimens were used. Specimens I. and II. were made from 
beryllium—copper strip of thickness 0-0125 in. Specimen I was formed by 
bending this strip into a symmetrical angle of side 0-576 in. and length 
10 in., the two legs of the angle being at right angles to each other. 
Specimen II. was a symmetrical channel of dimensions 0-633 in. web, 
0-242 in. flanges and length 16-5in. Specimen IIT. was a commercial 
steel pocket tape-measure of the type having a transverse curvature for 
stiffness. It was especially chosen from a number of such tapes for its 
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initial straightness, and its dimensions were : thickness 0-0065 in., width 
0-596 in., radius of transverse curvature 0-615 in., and length 17 in. Each 
specimen was provided with a specially made steel clamp which firmly held 
it at its root along the whole of its cross-section. 


Fig. 4. 


The testing apparatus is shown in fig. 4. The specimen was held 
vertically by its clamp, which was supported from the movable base-plate 
by means of the verticalrod. This base-plate could be moved about ae the 
horizontal table, and carried on its upper surface a piece of graph paper 
Attached to the bottom, free, end of the cantiliver was a horiecnt Laan 
piece that had a vertical pointer at some distance on each side of its centre 
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and a horizontal one projecting from each end. By sliding the base- 
plate about on the table the tips of these two horizontal pointers could 
be made to correspond with the tips of two other horizontal pointers 
that were fixed relative to the table. The apparatus was so arranged 
that when this was done the forces applied to the cross-piece—and 


Fig. 5. 


hence to the specimen—by the pair of horizontal strings, resulted in a 
pure couple about a vertical axis, while the tensions in the vertical strings 
' exactly balanced the weight of the cross-piece. In this way, by varying 
the equal weights on the two hangers, the displacements of the two 
vertical pointers relative to the graph paper, and hence to the support 
of the cantilever, could be measured for various angles of twist. It was 
hoped that these displacements could be read directly from the lower 
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ends of the vertical pointers, but the cross-piece tilted slightly as the 
specimen distorted, and this introduced an appreciable error due to the 
length of the pointers. So, instead, a small set square was used, with 
which, by sighting along the lines on the graph paper the coordinates 
of the centres of the pointers at the level of the end of the specimen 
could be determined. This method was found to give consistent and 
satisfactory results. 

In figs. 5, 6, and 7 are shown typical results for each specimen. As the 
tip of the cantilever rotated in an anti-clockwise direction the two vertical 
pointers traced out paths on the graph paper. In each case the cross- 
section is shown in its initial position. The full lines show the paths 
predicted by the analysis of this paper, while the broken lines show the 
paths that would be traced out if the tips of the cantilevers experienced 
only rotations about their centres of twist. The crosses are experimental 
points. 


Fig. 7. 


It Is seen that in each case the agreement between experiment and the 
analysis given in this paper is satisfactory, while the distinction between 
this theory and the assumption that the sections rotate about their centres 
of displacement is most marked. 
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SUMMARY. 


The paper discusses the theory of the reflection of very low frequency 
electromagnetic waves at the boundary of an ionized medium. The 
medium is assumed to be sharply bounded and homogeneous, and to 
have a steady magnetic field superimposed. The quasi-longitudinal 
approximation to the magneto-ionic theory is used (Booker 1935), and 
the reflection coefficients are then independent of the horizontal 
direction of the transmission path. The results of the theory are compared 
with the experimental observations of the propagation of very low 
oh EPRI TS SSS SETAE 0 Se eS Ee 
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frequency radio waves (Bracewell et al. 1951). It is concluded that the 
theory will account for the general form of some of the experimental 
results for frequencies of 16 ke./s. and above. It is likely to be most 
useful for extending the theory of the propagation of very low frequencies 
(10 ke./s. and lower) such as are used in the study of the propagation of 


§ 1. INTRODUCTION, 


1.1. Previous Work on the Theory of the Reflection of Very Low 
Frequency Waves. - 


A CONSIDERABLE amount of information now exists on the nature of the 
reflection of very low frequency radio waves from the ionosphere (Budden, 
Ratcliffe and Wilkes 1936, Bracewell et al. 1951). By “very low” 
we mean frequencies from about 10 ke./s. to 100 ke./s. It is of interest 
to consider what structure of the ionosphere would account for the 
observed phenomena. Several attempts at this have been made. 
A number of authors assumed that the ionosphere was a sharply bounded 
homogeneous medium containing electrons subject to collision damping, 
and neglected the effect of the earth’s magnetic field. Thus Yokoyama 
and Namba (1932) showed that for such a medium there is an angle 
of incidence, the quasi-Brewster angle, for which the reflection coefficient 
is a minimum. They attempted to explain the minimum which is often 
observed near sunset in the signal from a distant long wave station, 
by a change in the quasi-Brewster angle during the passage from day 
to night conditions. Bremmer (1949) and Rydbeck (1949), using similar 
assumptions about the ionosphere, were able to account for many of the 
observed facts relating to the propagation of very low frequency waves 
to great distances. However, it is known that when a plane polarized 
wave is incident nearly vertically on the ionosphere, the polarization 
of the reflected wave is different from that of the incident wave. This 
cannot be accounted for if the reflecting medium is isotropic, that is—if 
the earth’s magnetic field is neglected. 

Another group of authors has treated the ionosphere as a horizontally 
stratified medium, in which the constants vary with height in a distance 
small compared with one wavelength. The differential equations for the 
case where the earth’s magnetic field is not neglected were first given by 
Hartree (1929). Solutions for particular models of the ionosphere have 
been given by Wilkes (1940 and 1947) and Stanley (1950). 

It is the purpose of the present paper to give some numerical solutions 
for the case where the ionosphere is a sharply bounded homogeneous 
medium with a superimposed magnetic field, and to discuss their 
applicability to the experimental results. A problem very similar to 
this was discussed by Lorentz (1909) and Voigt (1915), who were concerned 
with the theory of the magneto-optical Kerr effect. This effect concerns 
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the change of polarization which occurs when light is reflected from the 
metal pole piece of a magnet. Lorentz and Voigt both derived a set 
of four reflection coefficients similar to equations (22) to (25) of this 
paper, but they discussed a reflecting medium with properties very 
different from those of the ionosphere. Bremmer (1949) has given 
expressions for the various reflection coefficients, but points out that 
their use in numerical computations is extremely complex. The few 
numerical results which he gives are discussed in § 4.4. They are based 
on values of the parameters which are probably different from those 
which would be expected to apply at the level of reflection of very low 
frequency waves. 


1.2. Scope of the Present Paper. 


The experimental results on the propagation of very low frequency 
waves refer mainly to frequencies of 16 ke./s. and higher. We shall show 
that for these frequencies some—but not all—of the experimental results 
can be explained by assuming the ionosphere to be homogeneous and 
sharply bounded. We conclude that at the level of reflection of these 
waves the variation of the characteristics of the ionosphere, in a vertical 
distance of one wavelength, is too small for the assumption of a sharp 
boundary to be valid. It would be necessary to adopt a full wave solution, 
of the type given by Wilkes (1940, 1947) and Stanley (1950), allowing 
for the finite rate of variation of the electron density and collision frequency 
with height. At still lower frequencies, however, such as those used in 
the study of radio atmospherics (Bureau 1930, 1932, Gardner 1950, 
Bowe 1951), the wavelength is very long (30 km. or more), and we may 
expect the theory of the sharply bounded ionosphere to be of some value. 

In § 2 of this paper we discuss the application of the Appleton-Hartree 
magneto-ionic formula to the region of the upper atmosphere, at about 
70 to 90 kms. where very low frequency waves are reflected. The 
quasi-longitudinal approximation (Booker 1935) is used. In § 3 we use 
the boundary conditions at the surface of the medium to deduce 
expressions for the various reflection and transmission coefficients. It 
is convenient to describe the incident and reflected waves in terms of their 
component electric fields parallel to, and perpendicular to, the plane of 
propagation, and to consider four coefficients denoted by 
1By> &,> Ry, ,R,, to indicate the complex ratio of a specified electric 
field in the wave after reflection to a specified electric field in the wave 
before reflection. The first subscript denotes whether the electric field 
specified in the incident wave is parallel (Il) or perpendicular (1) to the 
plane of incidence and the second subscript refers in the same way to 
the electric field in the reflected wave. In §4 the results of some 
numerical calculations are given and their relation to the experimental 
results are discussed. In §5 we discuss the validity of the various 
approximations used. 
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§2. THe ForMUL4 OF THE Magneto-lonic THEORY. 
The complete Appleton-Hartree magneto-ionic formula may be written 
(Appleton 1932) :— 
| (u2— 1) =1—jz—4y}|(1—@—ja) +4 (yr + 4yi(l —2@— je) 28) 
(1) 
where »=complex refractive index of the medium. 
x—=po|p*. 
Y=PulP- 
2=v/p. 
p=2n x frequency of the wave. 
pre=4rNe*/eym. 
Py=HpHe/m=gyro frequency of electrons in the earth’s magnetic 
field 
Yy—=component of y perpendicular to direction of wave propagation. 
y;,=component of y in direction of wave propagation. 
Pr=Yx/P- 
Pr=Y1/P- 
e, m=charge and mass of the electron. 
N=number of electrons per c.c. 
v=collision frequency. 
H=strength of earth’s magnetic field. 
€, 4o—electric and magnetic permittivities of free space. (This is 
the only place in this paper where py is used for the magnetic 
permittivity of free space. The same symbol is later used 


for the refractive index of the ordinary wave, but there 
should be no confusion.) 


In this paper we consider the ionosphere to be a homogeneous medium 
governed by equation (1), with a sharp lower boundary. For this medium 
there are two characteristic waves, the ordinary and extraordinary 
waves, corresponding to the two values of » given by (1). For each of 
these waves the value of » and the state of polarization depend on the 
direction of propagation. This would make the calculation of reflection 
coefficients very involved for the most general case. 

Kquation (1) becomes much simpler if it is permissible to make the 
* Quasi-longitudinal ” approximation (Booker 1935). This approximation 
is valid if 

| ef 40% | << Lia? ee ee eats | 
Then (1) reduces to 


pe=1—al(l—jety,) . . . . 2 2.” (8) 
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and it can be shown that the ordinary and extraordinary waves are 
circularly polarized with opposite senses (Booker 1935). Both waves 
also have a component of the electric vector in the direction of propagation, 
but it will be shown in § 5.3 that this component is usually small. The 
quasi-longitudinal approximation is used in obtaining all the numerical 
results given in this paper. Its applicability is discussed in § 5.1. 

The value of » given by (3) still depends upon the direction of 
propagation, on account of the term y,;. We shall at first, however, 
assume that y; is a constant, and the validity of this assumption will 
be discussed in § 5.2. Moreover, at the low frequencies considered here, 
Y;,>1 and we may write (3) 


p= 1 —9ju/(2 Jy.) =1—9( p,/p) exp (+7) see ears) 


If y;, is purely real, then 


p,,p=a V(F+yi)=PpVw +p) . . - - - ~ (5) 
tan r=y,,/2 =p,/v Shy SiR RS a esate ia Ca) 


Later, in § 5.2, we shall consider the possibility that y, may be complex, 
and the expressions for p, and 7 are then more complicated. 

We now take coordinate axes, ¥1, %j, 3, where x; is measured vertically 
upwards. We consider a wave with its wave normal in the (#,, 7) plane, 
inclined at an angle, 6, to the x, axis and travelling in the direction of 
both a, and x; increasing. Let the variation with time, #, be given by a 
factor e/”". We denote by E, the component of the electric vector of 
the wave in the (a1, #3) plane, and by E, the component in the direction 
of x. For such a wave in the northern hemisphere it can be shown that 
E,/E,=—j, or +j, for the ordinary and extraordinary waves, 
respectively. 

When a plane wave is incident on the ionosphere from below, the 
disturbance in the ionosphere consists of two waves, the ordinary and the 
extraordinary, for which the values of @ are different. We shall therefore 
write 

. For the ordinary wave : 


Complex refractive index=p)=1 —j(p,| per” SLRs tare ee (7) 
Direction of wave normal—6) 


Bolriztioneling Wyo ee ae fai eiere = <i bo eg? hens) (8) 


For the extraordinary wave : 
Complex refractive index=p,—1—j(p,/p)e""”  . . . . . + (9) 
Direction of wave normal—6, 


Polarization E,,/K,, =-+) Bip cee teeta s 6 Free ( 10) 
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§ 3. THE REFLECTION AND TRANSMISSION COEFFICIENTS 
FOR THE LONOSPHERE. 


Now consider a plane wave incident from below on the boundary of 
the ionosphere, with its wave normal at an angle 7 to the vertical, and 
with its electric field components given by : 


aa exp [jp{t—(a, sini+a, cosi)/c}]. . . . . (11) 


Le 


The reflected wave may similarly be denoted by the expression 


ee exp [jp{t—(x, sin i—a, cost)/c}]. . . . . (12) 
Abe 

In the ionosphere, the two transmitted waves are described by the 
expressions 


Ej exp [Jp tt— (a, sin 0) +23 COS O))Ho/¢ }] 1. at uh 3) 
and E,,, exp [jp{t—(x, sin 0,+-x3 cos 8,)u,/c}] ivy hy. eee) 


together with equations (8) and (10). 
The variation of all these waves in the x, direction at the boundary 
must be the same, and hence we have Snell’s law 


SiN t==p5 SiN Og==p, SUN 0, 7 ee a ieue cen OD 


It should be noted that because pz» and p, are, in general, complex, the 
angles 6), and @, are complex. Consequently y;, is complex. We shall 
investigate the consequences of this in § 5.2, where we show that, for 
practical purposes, it is sufficient to ignore the imaginary part of y;. 

Let E,, H, and E,, H,, be the components of the electric and magnetic 
fields K and 1 parallel to the xz, and a2, axes respectively, in either 
medium. Then the boundary conditions that must be satisfied are :— 


EK, and E, continuous . ..:. ..°. (16) 


H,; and H, continuous © “25. a 1) ere 
Application of these to the wave expressions gives the following four 
equations :—~ 


Ky; cos i— E,,. cos 4 = Ey cos 6) + E,,cos@, . . (18) 
K,; a La =—j Ky tf Bilas oc tommet rome L a) 
Ki; +E, = LoL mi MoH, AF indie Lea ae) 


E,,; cos i— EF, , cos ¢ =—jygHyp Cos 0) +ju,Ey, cos 0, . . (21) 


From these we can evaluate the reflection and transmission coefficients. 
For example, to find ,R,, we put E,,—0, and eliminate E,,, Ey) and 
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K,,, from the equations. The result is an equation for E,,/E,;, which is 
,8,. The results are as follows :— 


By = {(49 + He)(COS® i— Cos By COS A.) + (t44e—1)(COs +008 6,) cos 7}/D 


(22) 
R= 2) Cos (49 cos 97—p, cos 0,)/D Rae euteee acne We smut (23) 
R= 2) 60s i(9 Cos O.—PM,e COS A9)/D Me eRe ee Veet peak od) 
1Ry= {((Ho+H,)(Cos? i— cos 9 COs 8.) — (Hof4-—1)(cos O)+-cos 8,) cos +}/D 

(25) 


where 


D= (Uo +#¢)(COs? 1+ COs Oy Cos O,)+ (HoH, 1)(Cos O)-+ cos 8,) cosz. . (26) 


We may similarly evaluate transmission coefficients at the surface of the 
ionosphere. These will be denoted by To, ,T,, ,T, ,T,. Here the first 
subscript indicates whether the incident wave has its electric vector 
parallel (Il) or perpendicular (1) to the plane of incidence, and the second: 
subscript indicates whether the transmitted wave is the ordinary (0) or 
extraordinary (€) component. The transmitted wave is assumed always 
to be described by that component of its electric,vector which is parallel 


to the plane of incidence, so that 


ilo = Hyo/Ky3= JE yo/By; iTy= E/E = jE /E,; 
il,=E,,/E,,= —j#,,/Ky,; il. =E,,/E,;= =o) Bi 
The values of the transmission coefficients are given below :— 
nlop= 2/9 C08 4( Cost--p,cosé@,)/D . . .- . (27) 
he ie COS t\ mE COS ¢-1-2,/COB Op)/ 1) 0.5 4 + (28) 
70 = 2)H9 COS Us, cos 2--— cos 0,)/D. . . -. . (29) 
Ve=— Zit, COS t(4y COS*+ cCosG)/D . . . . (30) 


where D is given by (26). 
As a partial check of the formule (22) to (30), it is readily verified 
that they reduce to the Fresnel formule for a homogeneous medium, 


when we put po=p,, and O)=9,. 


§ 4. NUMERICAL RESULTS. 


4.1. The parameters p,/p and r. 

The preceding formule have been used to calculate values of Ri, and 
iR, as functions of the angle of incidence, for various values of the 
parameters p,/p and 7. The values of these parameters are given 
approximately by equations (5) and (6), although, as we have already 
mentioned, some slight modification may be necessary when account 
is taken of the fact that y, may be complex. This modification will be 
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discussed in § 5.2, and will be shown to be small. Equation (6) shows that 


tan 7 is proportional to the magnetic gyro frequency for electrons, and 
inversely proportional to the collision frequency, v. Equation (5) shows 
that p, is equal to p2/+/(v®+-pz) where py is the critical frequency that the 
medium would have in the absence of collisions and magnetic field 
p, is therefore directly proportional to the electron density N. 

If the earth’s magnetic field is zero, 7 is zero, Some results for this 
case are shown in figs. (1) and (9). 

If the collision frequency v is zero, then t=90°, and Pr=p2|py. Some 
curves for this case, for two different values of p,/p are shown in 
figs. 7, 8, 15 and 16. 


TABLE [. 
Frequency Collision Electron 
Figure No. P,|P T p27 frequency v Density N 
ke./s. se0.7) em. 

1,9 1 0° 16 (v/N=3-17 x 104) 
1,9 1 0° 100 (v/N=5-08 x 103) 
2, 10 3 15° 16 2-7 107 0-44 x 108 
2,10 4 15° 100 2:1 GLE 2-75 x 108 
Sil 2 15° 16 ark AOE 1-76 x 108 
ord 2 15° 100 2-7 X10" 11-00 x 108 
4,12 2 60° 16 4-2 x 108 0-53 x 108 
4,12 2 60° 100 4-2 x 10° 3°30 x 103 
5,13 4 60° 16 4-2 x 108 1-06 x 108 
5,18 4 60° 100 4-2 x 10° 6-62 « 108 
6, 14 8 715° 16 1-9 10° 1-85 x 108 
6, 14 8 75° 100 1-9 108 11-60 x 108 
Tal 2 90° 16 0 0-46 x 108 
7,15 2 90° 100 i) 2-88 x 103° 
8, 16 8 90° 16 0 1-83 x 108 
8, 16 8 90° 100 0 11-40 x 108 


We should expect the ionosphere to have characteristics somewhere 
between the extreme cases just mentioned. To consider the behaviour 
of a particular medium, we first decide what values to assume for 
v and p,, and derive the corresponding value of 7 from equation (6). We 
can then fix p, if we know the critical frequency of the medium p,. Any 
pair of values of p, and 7 therefore refers to fixed values of vy and N. 
The complex refractive indices can then be found from equations 
(7) and (9), using the known values of p,, 7, and the wave frequency 
p27. p, and p enter these expressions only as a ratio. Since it is 
convenient to discuss the results in terms of particular wave frequencies, 
the values of N and v appropriate to the curves of figs. 1 to 16 for selected 
wave frequencies are listed in Table I. The frequency 16 ke./s. has been 
included because many of the experimental results apply to it. The 
yalues of vy and N in Table I refer to vertical incidence when the angle of 
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magnetic dip is 60°. It will be shown in §5.2 that only small 
modifications are needed at other angles of incidence, and that the form 
of the results is independent of the relative directions of the transmission 
path and the magnetic meridian. 


4.2. The Reflection Coefficients R,, and ,R c 

Figs. 1 to 8 show some typical numerical values for the amplitudes of 
uR, and ,R,, plotted as functions of the angle of incidence i. Figs. 9 
to 16 show the corresponding curves for the phases of ,R, and ,R Ry. 
This last quantity gives the relative phases of the normal and abnormal 
components of the reflected wave, and thus determines its state of 
polarization. 

Figs. 1 to 8 show clearly that there is an angle of incidence for which 
|,R,| is a minimum. This is the quasi-Brewster angle, which exists 
also when the superimposed magnetic field is zero. Figs. 9 to 16 show 
that there is a change in the phase of ,R, as the angle of incidence 
varies. The change is most rapid near the quasi-Brewster angle. 


4.3. Comparison with Experimental Results. 

The experimental results for a frequency of 16 ke./s. may be summarized 
as follows :—(Bracewell e¢ al. 1951). For an angle of incidence of about 
30° the wave reflected from the ionosphere is roughly circularly polarized 
with a left-handed sense. The reflection coefficient | ,R,| at night is 
0-4to 0-5. Inthe daytime it is about 0-25 for winter and 0-1 for summer, 
and increases as the angle of incidence is increased. 

The theoretical curves of figs. 10 to 16 show that the polarization 
of the reflected waves would be elliptical with a left-handed sense. The 
theoretical curves for which the polarization, at 1=30°, would be most 
nearly circular are those of figs. 4 and 12. Here the phase difference 
between ,,R, and ,R, is about 240°, whereas for circular polarization 
it should be 270°. The corresponding value of | ,R,| is about 0-34, so 
that for this particular angle of incidence the curves of figs. 4 and 
12 represent the experimental results fairly well for a winter day. The 
appropriate values of the parameters are: p,/p=2, and 7=60°. 
Table I. shows that these correspond to the values : v=4:2 x 108 sec.—t and 
N—530 cm.-3. There are no values of the parameters which will fit the 
experimental results much better than these. 

The curves of figs. 10 to 16 show that as the angle of incidence increases, 
the polarization of the reflected wave changes, and for very oblique 
incidence becomes nearly linear. An effect of this sort has been reported 
by Weekes (1950), whose experiments show that for very oblique 
incidence the reflected wave is roughly linearly polarized with its plane 
of polarization at about 45° to the vertical. Weekes’ experiments do 
not give any information about the direction of tilt of the plane of 
polarization, but it would appear from the results of some experiments 
made by Mr. Bain in Aberdeen, that the observed direction of tilt is 


opposite to that given by figs. 10 to 16. 
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Curves showing reflection coefficients |,,R,| and |,R,| and phases of R, and 
1X,/,Ry as functions of angle of incidence. Quasi-longitudinal approxi- 
mation used. Results are independent of horizontal direction of trans- 
mission path. 
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In fig. 17 we show the values of the reflection coefficients |R,| and 
| ,R, | plotted as functions of frequency for an angle of incidence 30°. It 
is seen that | R,| decreases as the frequency increases, and this is in 
qualitative agreement with the observations. The theoretical decrease 
of |,,R,| with frequency is very slight, and the observed decrease is 
much more rapid. 

We conclude that, if we assume the ionosphere to be a sharply bounded 
homogeneous medium to which the quasi-longitudinal approximation 


Fig. 17. 


COEFFICIENT 


REFLECTION 


Ove 0-4 0-6 O-8 ee) 


P/b. (PROPORTIONAL TO FREQUENCY ) ——> 


Variation of reflection coefficients |,,R,| and |,,R,| with frequency. 
Angle of incidence i=30°. 760°. 


is applicable, we cannot account for the observed phase difference of 
90° between the normal and abnormal components of the reflected wave. 
For an exact theory a more complete treatment is necessary. For 
example, we might investigate the properties of an ionosphere with values 
of N and v for which the quasi-longitudinal approximation is not valid 
but Bremmer (1949) has shown that the analysis of this case is extremely 
complex. Alternative treatments taking account of the finite rate of 
variation of N and v with height have been made (Wilkes 1940 and 1947 
Stanley 1950), and, for reasons given in § 4.6, these are probably preferable 
to any extension of the theory of the sharply bounded ionosphere, 
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However, the sharply bounded model discussed in this paper gives a 


good rough approximation to the behaviour of the ionosphere for very 
low frequency waves*. 


4.4. Comparison with Bremmer’s Results. 


Bremmer (1949: Terrestrial Radio Waves, pp. 286 to .301) has 
discussed the reflection of waves at a sharply bounded, homogeneous 
non-isotropic ionosphere using the full magneto-ionic formula, without 
making the quasi-longitudinal approximation. He gives numerical results 
for two cases only. He assumes the earth’s magnetic field to be vertical, 
with py=9-4x< 108 sec.'. His results are given in Table II., together 
with the corresponding values of 7, p,/p and N. (Bremmer assumes the 
waves to have a time factor ¢/”. The signs of arg (,R,) and arg 
(2, /,R,) have been adjusted to agree with the time factor e/” and the 
system of axes used in the present paper.) 


TasB_eE II. 
Angle of Incidence 76° Ave 
Frequency (ke./s.) 300 30 
a 84° 84° 
v (sec.—1) 108 10° 
P,P 13 29 
N (cm.-) ao LOS 1-6 x 104 
Bul 0-231 0-427 
arg (Ry) —138° —55° 
Ry | 0-214 0-166 
arg (yRy/yRy) 312° 256° 


It should be noted that the values assumed for the electron density 
N are very large, and correspond to critical frequencies of 2-4 Mc./s. and 
1-1 Mc./s. respectively. These might apply somewhere near the maximum 
of the E-region, but are much too large for the levels, usually between 
65 and 95 km., where very low frequency waves are reflected. The 
values of |y?/(l—a«—jz)?| are of the order of -0004 and -007 respectively 
for the two cases discussed, so that the quasi-longitudinal approximation 
might have been applied with fair accuracy. 


4.5. The Transmitted Waves. 

It may be of interest to know the ratio of the amplitudes of the 
ordinary and extraordinary waves transmitted into the ionosphere. 
The values of this ratio for various values of the parameters have been 
calculated from (27) and (28) and are shown in the following Table IIT. 


* It is hoped to use this model to discuss the effects of the earth’s magnetic 
field on the propagation of a radio atmospheric, and it was originally developed 
for this purpose. For a discussion of the propagation of atmospherics which 
neglects the earth’s magnetic field, see Budden 1951. 
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It is found that the ratio is always fairly close to unity. It is also of 
interest to know the direction of propagation of the two waves in the 
ionosphere. Let ¢) and 4, be the inclination to the vertical of the 
normals to the planes of constant phase for the ordinary and extra- 
ordinary wave, respectively. Then 


tan dy= &(sin 0))/A(cos 9) ; tan d= A(sin 6,)/ B(cos 8.) 
The values of ¢, and ¢, are also shown in Table ITI. 


TaBxe III. 
Angle of incidence 10° 80° 10° 80° 10° 80° 
Pri D 2 2 2 2 8 8 
| T 15° 15° 60° 60° 75° 75° 

Ratio of amplitudes, ordin- 

ary wave to extraordinary 

wave . 1:04 | 1:05 | 1:09 | 1-40 | 0-84 | 1-06 
Inclination )) ordinary 8° 20° | 56° 12’ | 4° 50’ | 39° 20’| 4° 4’ | 19° 10’ 

of wave — 

normal extraordinary | 9° 46’ | 59° 18’ | 3° 40’ | 17° 52'| 1° 22’| 3° 4’ 


4.6. The Thickness of the Ionosphere. 


It is of interest to ask how thick the ionosphere must be for the 
assumption of homogeneity to be justified. The two transmitted waves 
are both attenuated as the height increases in the ionosphere, and it is 
found that the ordinary wave is the least attenuated component. The 
height at which the amplitude of the ordinary wave is reduced to one-half 
of its value at the boundary has been calculated for some selected values 
of the parameters. The results are shown in Table IV. 


TABLE IV. 


Angle of incidence 0° 


P,|_P 


a 


Thickness of medium in 
wavelengths for re- 
duction to half 
amplitude 
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The above results perhaps provide the strongest reason for rejecting 
the sharply bounded ionosphere when an accurate theory of the reflection 
of very long waves is required. Thus, for a frequency of 16 ke. /s. (wave- 
length 18-75 km.), the last row of Table IV. shows thicknesses of the 
order of 6 km., and there is little doubt that there are large variations 
of N and » in distances of this order at the level of reflection of very low 
frequency waves. 


§ 5. VALIDITY OF THE APPROXIMATIONS. 
5.1. The Quasi-longitudinal Approximation. 
We see from equation (2) that the quasi-longitudinal approximation 
is valid if 
lyf fy —e—je} ed. 1) 
The values of this quantity for some selected values of p,/p, 7 and i, are 
shown in Table V. 


TABLE V. 
P,P T Angle of incidence 7 | |y4/{4y2(1—a—jz)?}| 

4 15° S04 0-21 

2 60° 0° 0-0044 
2 60° 30° 0-031 
2 60° 60° 0-12 

2 60° 90° 0-19 

4 60° 30° 0-0055 
1 ghi 30° 0-79 

8 75° 90° 0-0039 
2 90° 90° 0-24 

8 90° 90° 0-0031 


These values all refer to propagation from south to north in the 
magnetic meridian, since it is for this case that the quasi-longitudinal 
appréximation would be least accurate. It is found that the quantity (31) 
is greatest for large angles of incidence, for large values of 7, and for 
small values of p,/p. It is clear that the quasi-longitudinal approximation 
may be treated as valid for most of the range covered by the curves of 
figs. 1 to 16. It begins to be inaccurate at high angles of incidence in 
figs. 2, 10, 4, 12, 7 and 15. 

It might be suggested that further curves of the type of figs. 1 to 
16 could be calculated using the quasi-transverse approximation, which 
is valid if 

| ya {4yi(1 —e—Jjz)?} [> 1. 
However, this condition requires values of v rather lower than 
2 108 sec.—, and this is less than would be expected to apply at the level 
of reflection of very low frequency waves. Moreover, for propagation 
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in the direction of the magnetic meridian the quasi-transverse approxi- 
mation would show that the reflected wave should always be linearly 
polarized in the plane of incidence, whereas the experiments (Bracewell 
et al. 1951) show that it is roughly circularly polarized for small or 
moderate angles of incidence. 

We conclude that as long as p,/p is large enough, the quasi-longitudinal 
approximation represents fairly well the behaviour of a sharply bounded 
ionosphere. It would therefore be. the best approximation to use in 
discussing the propagation of waves of very low frequency (10 ke./s. or 
less) such as are used in the study of the propagation of atmospherics 
(Gardner 1950, Bowe 1951). 


5.2. Variation of Parameters with Angle of Incidence. 

If d is the angle of magnetic dip, and we consider transmission in a 
horizontal direction ¢, measured east from magnetic north, then 
Y;, is given by 

y,=y(sin 0 cosdcos¢-+cos@sind), . . . . (32) 
where @ may be either 0) or 6... We shall assume in what follows that 
d=60°. 


TaBLE VI. 
Ordinary wave Extraordinary wave 
Angle of incidence 0° 30° 90° 0° 30° 90° 
T 61°4' | 60° | 56° 23’| 58° 18’ 60° 64° 40’ 
P,|P 3°95 4 4-18 4-05 4 3°89 


It was mentioned in § 3 that when Snell’s Law (equation (15) ) was 
used to derive the directions 6) and 6, of the two waves in the medium, 
these angles proved to be complex, because the refractive indices py) and 
uw, are complex. As a result y, is complex. We now investigate the 
consequences of this. 

If y;, is a complex number, the denominator of equation (4) must be 
written 


ZF I (y)tjAyy) - . - . . . « (83) 


(# means “real part of” and % means “ imaginary part of”). Both 
Aly, and 4(y,,) vary with angle of incidence and the result is a change 
in the values of p,/p and 7+. The modification is most serious at very 
oblique incidence, but even then the associated changes in the parameters 
p,|p and 7 are small. If we assume that r—60° and p,/p=4 are correct 
for an angle of incidence 30°, then for west to east transmission in the 
northern hemisphere, and for a frequency of 16 ke./s. the correct values | 
of 7 and p,/p for other angles of incidence are as shown in Table VI. 


‘ 
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The above variations of 7 and p,/p with angle of incidence are not 
sufficient to affect the general form of the results. It is sufficiently 
accurate to regard fixed values of p,/p and 7 as referring to fixed values of 
N and » for all angles of incidence. This is equivalent to treating y,, as 
a constant. An important consequence of this is that the reflection 
ee aponts are independent of the horizontal direction of the transmission 
path. 


5.3. The Longitudinal Component of the Electric Field. 


The whole of the preceding theory has been worked out on the 
assumption that the waves in the ionosphere had only transverse electric 
fields. It is known, however, that the electric fields also have components 
in the direction of the wave normal. We shall call these the longitudinal 
fields and denote them by E,, and E,,, respectively, for the ordinary 
and extraordinary waves, or simply by E,, when the formule apply to 
either wave. It can be shown that if the quasi-longitudinal approximation 
is valid, then for either ordinary or extraordinary wave 


| Hy /Ey [=| 4?—1].| ye/(1—a—jz)|. . . . (84) 
TaBLe VII. 
p,|p | 7 | Angle of incidence | |E,/E,,| | H,/E,|.|tan 6] | Component 
4 15° 90° 0-15 0-21 ext. 
2 60° Ox 0-48 0 ext. 
2 60° 30> 0:74 0-14 ext. 
2 60° 60° 0-99 0-50 ext. 
2 60° 90° 1-08 0-66 ext. 
+ 60° 30° 0:67 0-036 ord. 
I 75° 30° 1-32 1-1] ext. 
8 75° 90° 4:2 0-53 ext. 
2 90° 90° 2-2 2:7 ord. 
8 | 90° |. 90° 1-2 0-18 ord. 


The longitudinal electric fields will in general have components 
parallel to the boundary of the ionosphere. This would mean that 
additional terms should appear in the boundary condition equation (18). 
The ratio of the magnitude of the additional terms to those already on 
the right hand side of (18) is proportional to the expression (34) and to 
tan 6, where @ is either 0, or @,. 

The ratio‘is therefore small for nearly vertical incidence, for then 
6 is small and the longitudinal electric field is nearly perpendicular to the 
boundary. The values of | E,,/E,,| and of | E, [E,|.| tan @| for some other 
cases are shown in Table VII. These all refer to propagation from south 
to north in the magnetic meridian, and the figures apply to that wave, 
ordinary or extraordinary, for which the value is the larger. 
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A comparison of Tables V. and VII. shows that the errors which arise 
by neglecting the longitudinal electric field become serious for about the 
same values of the parameters as those for which the quasi-longitudinal 
approximation becomes inaccurate. It is found that for very large 
values of p,/p the errors are small, which provides further encouragement 
for using the theory in discussing the behaviour of the lowest frequencies 
(10 ke./s. and lower). 
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LXXXVI. Note on the Analytic Continuation of the S-matriz. 


By N. G. van Kampen, 
Institute for Theoretical Physics, Copenhagen, Denmark*. 


[Received April 23, 1951.] 


SUMMARY. 


It is shown that methods which try to avoid the redundant poles are not 
likely to be successful. On the other hand, the analytic continuation of the 
S-matrix is an ill-defined procedure, which cannot be applied until more 
properties of the S-matrix are known. The redundant poles are caused by 
the use of the asymptotic expression of the wave function instead of its 
actual value and have no physical meaning. 


§1. PRELIMINARIES. 


THE scattering of a non-relativistic particle by a fixed central field of 
force may be described by the equation 


u(r) +{E—V(r)u(r)=0 2 ww. CL) 


with the boundary condition u(0)=0. Let u(E, 7) be the solution, the 
arbitrary constant factor being determined by w’(H, 0)=1. For each r 
the function u(E, r) is an integral function in the complex E-plane. For 
real positive E we shall write E=k? with k>0. Then the asymptotic 
expression of u(E, 7) for r+ 00 is 


u(E, r)~O(k) sin {kr—n(k)} for E=R>0. . . . (2) 


The problem at issue is to find the discrete eigenvalues E,, of (1) from the 
phase function 7(é) only. 

Let I’ denote the complex E-plane with the exception of the positive 
real axis. For each value of E in I’ there is (apart from a constant factor) 
one solution of the equation (1) which goes to zero for r>oo. It can be 
defined uniquely as the solution f(E, 7) of (1) which is asymptotically 
exp [i(4/E)r], where 1/E is made unique in I" by the condition Img ./E>0. 
For each 7, f(E, 7) is a one-valued function of Ein I’. When E approaches 
the cut along the positive real axis from above, f(E, 7) tends to that solution 
for E=k2>0 which is asymptotically f(E+,7)~exp [tkr]; when E 
approaches from below f(E, r) becomes f(E—, r)~exp [—*kr]. 

As shown below (see also Jost 1947), the function f(E, 0)=f(E) is 
analytic in I” without singularities. Its zeros are the energies of the 
bound states ; they are necessarily on the negative real axis : 


f(Ez)=0;,E,, <0; Waeeameriaerr (3) 
* Communicated by Prof. C. Meller. 
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The discontinuity of f(E) on the cut is connected with the S-matrix by the 
relation (Jost 1947) 


S(E) =exp [—2in(k)]=f(E—)If(E+). . - . + @) 


It should be noted that as yet no appeal has been made to analytic 
continuation, with the consequence that S(E) has a meaning only for real 
positive E. Now consider a second sheet I’, of the complex E-plane, 
connected with I along the positive real axis in such a way that one 
enters J", on crossing this axis upwards. If f(E) can be continued analytie- 
ally in I’,, taking the values f,(E) say, then S(E)=/,(E)/f(E) is an analytic 
continuation of (4)in I’. From (3) it follows that S(E) has poles (in general) 
for E=E,,, which suggests that the discrete spectrum can be found from 
the analytic continuation of the S-matrix (4). 

It is clear, however, that S(E) may have “ redundant poles ” (Ma 1946), 
namely, if f,(E) has poles. Also, if f,(E,)=0 for any E,, then this 
eigenvalue is not exhibited as a pole of S(E). Because of these difficulties 
it is tempting to seek a method of determining the discrete spectrum 
from (4) without an analytic continuation of f(E). This was tried by 
Wildermuth (1950) without conclusive result. We shall now show why it 
is not possible. 


if 

L * 
° io) ° ! ct 
EP oe, O | 

le 


The complex E-plane. 


§ 2. No Anatytic ConTINUATION OF S8(E). 


Levinson (1949) in a different connection introduced the function 
F(E)=1+ | exp [i(4/E)r]V(r)u(E, r)dr, . 2... (5) 
0 


and proved the following properties. (i) F(E) is analytic in I’ and tends to 
1 at infinity. (ii) F(E) tends to a limit when E tends to a point E=k? on 
the positive real axis, namely 


F(E+)=hO(k) exp [in(k)], - F(E—) = C(x) exp [—#n(h)] 


(k chosen positive). (iii) The zeros of F(E) are simple and lie on the 
negative real axis; they determine the discrete spectrum. This last 
property suggests that F(E) is identical with f(E), which can be shown as 
follows. : 

From the constancy of the Wronskian determinant follows 


w(E, r) f(B, r)—u(B, r) f’(E, r)=f(E, 0)=f(E). . . . (6) 
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For large r the left-hand side becomes (if E is in I’) 
(u'(E, r)—i(4/E)u(E, 1) } exp [i(1/E)r] 


On the other hand, it can be seen that the latter expression is equal to 


1+] exp fil VE)rIVir)u(B, 1) dr, ee ak ee ee 7) 


‘ because it has the same value for r—0 and the same derivative for r>0. 


- 


Hence, by taking r= ©, one finds that (6) is equal to the limit of (7), which 
is (5). Thus it is proved that F(E)=f(E) and, consequently, f(E) is 
analytic in and tends to 1 at infinity. 

Now Levinson showed that F(E) can be expressed in terms of the 
function 7(k), when F(E) has no zeros, 7. e. when there are no bound states. 
We shall now derive a similar relation for the case that there is a finite 
number of bound states. In that case, the function 


E, 
1B) =loe yA =log f(E)—2 log (1- =) a3) 


is one-valued analytic in I" (because the zeros of f(E) are simple) and 
vanishes at infinity (when the appropriate branch of the logarithm is 
chosen). It has a discontinuity along the cut : 


g(E+)—g(E—)=log f(E+)—log f(E—)=log S(E). . . (9) 


By means of the Cauchy integral g(E) can be expressed as an integral 
along the path L in the figure, with the result (for E in I’) 


| awe Lf a 


Inserting (8) and (9) one obtains 


E 
te ee ie Bet an’ +2 log(1— Ft). Ee 6) 


g(E)= Qi Qt 


on i’ 


where all logarithms are chosen in such a way that they vanish when 
E (or E’) tends to infinity. 

From (10) it is clear that the values of S on = positive real axis 
together with those of E,,, determine the function f(E).* On the other hand, 
when only S is given, the values of E,, may still be chosen arbitrarily, 
because for each choice an f(E) results from (10), which has the zeros E,, 
and the required a eae (4) along the cut. Hence it appears aa 
the discrete spectrum cannot be determined by using ei function f(E 
inside T° only, so that the analytic continuation of f(E) beyond the ee 
becomes inevitable. 


* The same conclusion was reached by Jost (1947). 
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It may be that f(E) must satisfy some further requirements, which 
reduce the arbitrariness of the choice of the E, in (10). For instance, it 
follows from (10) that for small E 


log f(E)=— ‘oni log 8(0)+N} log E-+finite function, 


where N is the number of zeros E,,. Now, according to its definition, 
f(0) is finite when rV(r) is finite for r=0. In that case, therefore, the 
choice of E,, is restricted by the condition 


—(1/2zi) log S(0)=n(0)/rr. 


(If N<n(0)/7, then f(0)=0 and there is a bound state with energy 0 
(Bargmann 1949).) However, it is unlikely that the choice of the E,, can be 
made unique by such general conditions. Consequently, our conclusion 
can still be maintained. 

Wildermuth’s method (1950) consists of defining an auxiliary function 
with a parameter «>0 


a)—=f(E—1«)/f(E+7«). 
This function is ate in the E-plane with the two cuts 7+ A and 


—ia+aA(0<A< oo). Hence it is determined by its values on the positive 
real axis. For E on the negative real axis one has ¢ 


lim S(E, «)=1 for «a0, 


except if f(E)=0. In that case one finds, by taking the derivatives of 
numerator and denominator, lim 8(E,,, «)=—1 (since f’(E,,)40). Conse- 
quently, the discrete spectrum ean be determined if S(E, ae is known for all 
E>0 and for « in a certain interval 0<a<a,. All the same this is not a 
solution of the problem, for S(E, «) itself is not determined by the phase 
function 7(k). It seems, Ronee premature to conclude from this 
(Wildermuth 1950, Heisenberg 1950) that the knowledge of (k) alone is 
insufficient for the determination of the E,, because there is still the 
alternative possibility of analytic continuation of S(E) (as in Kramers’ 
original idea). 


§ 3. Tor Usz or ANALytTIc CONTINUATION. 


The problem of finding the analytic continuation of a function given on 
the real axis is equivalent to the problem of solving the potential equation 
in two dimensions, with given values of the potential and the charge 
density on an open boundary. It is the initial-value problem for a partial 
differential equation of elliptic type, and can be solved only if the boundary 
function is given by means of some analytic expression (see for example 
Courant-Hilbert 1937). That is indeed the case when the function S(E) 
is obtained by solving equation (1) rigorously. However, in case the 
S-matrix is found experimentally, or from some approximate calculation, it 
is known only with a certain error. Then the analytic continuation is 
undetermined and the poles E,, cannot be found, 
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In order that the concept of analytic continuation be applicable to 
that case, it is necessary that S(E) is known on the whole boundary of the 
domain I’. Nowit appears from Ning Hu’s work (1948) that the behaviour 
of S(EK) at infinity can indeed be found from general considerations. 
However, the values on the lower side of the positive real axis (0. e.» 
S(E—)=f,(E--)/f(E—) for E>0) do not seem to follow from a priori 
arguments. Usually it is postulated (Heisenberg 1946, Ning Hu 1948) 
that S(E—)—S(E)*, but there is no guarantee that this is the value which 
follows from analytic continuation through I. Until this difficulty has 
been cleared it is impossible to find the bound states from the S-matrix in 
practical calculations. 

In the same way, when S(E) is given for E>0 with a finite precision, it 
is impossible to find the poles that do not correspond to bound states. 
Indeed, we know that they can always be introduced or suppressed by an 
arbitrarily small alteration of S(E) obtained by changing the potential at 
some large distance (Jost 1947). This change in the potential has no 
physical meaning, since any actual measurement is performed at a finite 
distance. Hence it is seen that the redundant poles are caused by using 
the asymptotic expression of the wave function instead of the actual value 
at the point of observation. One might ask whether this means that the 
S-matrix theory is based on an unjustified separation of physical systems 
or only that the equation (1) is too much simplified. At any rate, an 
uncertainty is involved of the same order as the error made in cutting off 
the potential (ter Haar 1946) at a distance comparable to the distance of 
the measuring instruments. Therefore it is legitimate to discard the 
redundant poles by this device. 


I am indebted to Prof. C. Moller for his interest and advice, to 
Prof. Niels Bohr for the opportunity to work at the Copenhagen Institute, 
and to the Rask-@rsted Foundation for a grant which enabled me to stay 
in Denmark. 


Note added in proof.—Recently P. Jauho (Ann. Acad. Scient. Fennicae, 
Al, nr. 80, 1950) has extended Levinson’s work so as to include a Coulomb 
potential. In the present article that case has been excluded because of 
eq. (2), but it can be checked that our conclusions remain valid. In 
particular, it is not possible to build an S-matrix theory on Jauho’s 
function x (the counterpart of Levinson’s F), as he suggested (p. 29). 
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LXXXVII. The Attenuation of Nucleon Cascades in Lead. 


By W. G. V. Rosser and M. W. Swirt, 
The University, Manchester*. 


[Received May 1, 1951.] 


ABSTRACT. 


The frequency variation of nuclear interactions under lead has been 
studied by comparing the rates of stars in Ilford G5 emulsions exposed 
above and below 30 ems. of lead. The size frequency distribution is the 
same for stars in both sets of plates. . 

The attenuation lengths in lead are :-— 


(305+7) gm./cm.? for all stars, 

(380+ 65) gm./cm.? for penetrating showers (n,>2), 

(405+31) gm./em.? for events with charged relativistic primaries, 

(260-34) gm./em.? for events with uncharged primaries (excluding 

simple evaporation stars). 
The ratio of the attenuation length to the interaction meanfree path for 
penetrating showers is (1-8--0-3). 

The percentage of all stars produced by charged relativistic primaries 
increases from (13-1+-9) per cent for the unshielded plates to (17-0-41-0) 
per cent for the plates under lead due to the formation of stars by 
a mesons generated in lead. The ratio of charged to uncharged 
penetrating-shower primaries increases from (0-86-+-0-21) for the air 
plates to (1:94--0-32) for the plates exposed under 30 ems. of lead. It 
is suggested tentatively that the increase is due to the production of 
penetrating showers by fast + mesons created in the lead. 


§1. INTRODUCTION, 


THE main properties of cosmic rays in the atmosphere can be accounted 
for in terms of nucleon cascades initiated by the primary cosmic rays in 
the upper atmosphere. The primary protons, «-particles and heavier 
nuclei collide with air nuclei and fast nucleons and mesons are emitted. 
Following Rossi (1948) the interacting particles will be called N-rays. 

It is important to define the interaction mean-free-path ), and the 
attenuation length A, of the interacting particles. The interaction 
mean-free-path is defined as the average distance an N-ray goes before 
giving rise to a nuclear interaction. The attenuation length is defined 
as the thickness of absorber in which the intensity of N-rays is reduced 
by a factor 1: e. It is found experimentally that in the lower atmosphere 
the variation of N-rays with altitude is roughly exponential so that 


I=I, (exp —a/A,) 
* Communicated by Professor P, M. 8, Blackett, F.R.S, 
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where x is the thickness of the absorber. In each nuclear interaction 
the primary particle loses energy so that in general the next nuclear 
interaction will be a less energetic event. The attenuation length is 
defined for the frequency variation of a particular type of event. It is 
not confined to the absorption of particles in a given energy range which 
produce these events, because after the first collision these particles will 
have left the energy range. The attenuation length therefore represents 
the decrease in intensity of a spectrum of particles which undergo 
collisions with energy loss, the attenuation being due mainly to the 
fact that there are fewer high energy particles to take the place of lower 
energy particles as the N-rays traverse matter. 

Following Fujimito and Hayakawa (1949), we shall call the primaries 
of penetrating showers’ the A-rays and the primaries of evaporation 
stars the B-rays. The division is not complete but is useful in discussing 
the nature of nucleon cascades. The average energy of B-rays is about 
1 BeV. and these rays are, in the main, secondary to the A-rays. 

The present experiment had the following main objectives :— 


1. To investigate the effect of lead on the A and B components 
separately. 

2. To compare the value of A, obtained for penetrating shower 
primaries under lead with the value of A; determined by counter 
experiments. 

3. To ascertain what fraction of stars under lead can be ascribed to 
fast 7 mesons whose mean-free-path for star production in lead is 
shorter than the mean-free-path for decay into » mesons. 


§2. EXPERIMENTAL PROCEDURE. 


Ilford G5 plates with emulsion 300 in thickness were exposed for 
21 days at the Pic du Midi, altitude 2860 m., one group of the 
plates (the ‘‘ lead ’’ plates) being placed under 30 cm. of lead built in the 
form of a hemisphere, and the other plates (the “air” plates) 
1-5 m. from the lead. Both groups of plates stood on the concrete floor 
of a hut with a thin aluminium roof, which was covered with about 
1 inch of snow during the exposure, making a total thickness of less 
than 2 gm./cm.”. It is estimated that this reduced the star rates by 
about 0-1 stars/c.c./day and, since both sets of plates were equally 
affected, no correction was applied. The plates were collected as soon 
as dry at Messrs. Ilfords, taken at sea level to Bagnéres-de-Bigorre and 
then transported as quickly as possible to the summit of the Pic du Midi. 
The plates were brought back to sea level immediately after exposure. 
Since the times at various altitudes were known accurately, the 
background rate could be calculated. It amounted to 3-3 per cent of 
the total number of stars observed on the air plates. 

The plates were developed in acidified Amidol, using the temperature- 
cycle method of Dilworth, Occhialini and Payne (1948). They were 
shown to be sensitive to minimum ionization by observation of ue decays 
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in the emulsion. The efficiency of detection of minimum ionization 
tracks was probably only about 90 per cent (Cosyns et al. 1949) because 
of the loss of tracks dipping steeply into the emulsion. Since, however, 
the main conclusions deduced from the present experiment depend on 
the relative rates of events above and below lead, if the same percentage 
loss of minimum ionization tracks occurs in both sets of plates the effect 
is not important. There is also a possibility of non-associated tracks 
passing so close to a star that they are accepted as being associated 
with the star. A strict criterion was adopted, namely, two minimum 
ionization tracks were not accepted as belonging to a star unless there 
was an observable deflection between the tracks at the star. In this way 
spurious events were cut to a minimum. 

The choice of thickness of emulsion was dictated by the need for fast 
and accurate scanning. The rate of accumulation of data increases with 
increasing thickness of emulsion, but the percentage of events missed 
also increases. After some test scanning on plates of different thickness 
it was decided that a thickness of about 300 was the optimum valu 
for our experiment. 

The’ results obtained, when the same area of plate was scanned using 

first a 10x objective, then a 20 objective, then a 45x objective, 
were compared. It was found that with the 10 x objective 20-30 per cent 
of all stars were lost, but with the 20x objective less than 2 per cent 
of all stars were missed. This was less than the statistical accuracy 
aimed at and, furthermore, the losses affected both sets of plates equally. 
The scanning was therefore carried out using 20 Cooke objectives and 
10 eyepieces. The stars were examined in detail with 45 x oil immersion 
objectives and grain counts were made on the tracks. To be accepted 
in the analysis a star had to have at least three prongs or two prongs 
and a short nuclear fragment. No two-pronged stars were accepted, 
as it was difficult in a large number of cases to distinguish between a 
two-prong star and a scattering of a heavy particle. Stars which had 
all tracks less than 65 » long were excluded on the grounds that they 
were probably of radioactive origin. No separation of o meson stars 
was made, and thus stars of two or more prongs produced by 7 mesons 
at rest were included in the analysis. 

The stars were classified in the same way as was used by Brown 
et al. (1949), and the results are shown in Tables I. and Il. The number of 
tracks with a grain density less than 1-5 times minimum is denoted by 
n,, and the number with grain density greater than 1-5 times minimum 
by Nx. 

§3. EXPERIMENTAL RESULTS. 

The results obtained with the air and lead plates are shown in 
Tables I. and II. respectively. The total rate of stars observed in the air 
plates was 9-0+-0:3 stars/c.c./day, after correcting for the background 
rate. This rate can be compared with rates obtained by Brown et al. 
(1949) and Bernardini et al. (1950) at 3500 m, Assuming the attenuation 
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length for all stars in the atmosphere to be 135 gm./com.?, the rate 
observed at the Pic du Midi is equivalent to a rate of 14:1-L0-4 
stars/c.c./day at 3500 m. The rate obtained by Brown et al. was 
10-8 stars/c.c./day. This value is significantly lower than the present 
result, but the difference is confined to 3- and 4-prong stars. The present 
result is, however, consistent with the rate of 14:2--0-5 stars/c.c./day 
observed by Bernardini et al. 

The rate observed under 30 cms. of lead was 2:96-0:09 stars/c.c./day 
after correcting for the background rate. 

Assuming that the absorption is exponential, the attenuation length 
of all stars in lead is 305+7 gm./em.2. This value compares favourably 
with the values of 310+20 gm./em.2 and 30020 gm./cm.? obtained 
by George and Jason (1949) and Bernardini et al. (1949 a) respectively. 

The attenuation lengths for the various types of events are shown in 
Table IIT. 

Tase IIT. 


Attenuation Lengths in Lead. 


Type of event Attenuation length 


gm./em.? 

CaN stars 305-+7 - 
Op+lp+2p+... 405 +31 
In+2n+3n+... me |? 260 +34 
2n4+-2p+3n+3p+... 380 +65 


(Penetrating showers) 


The significance of these results will be discussed in the next section. 


§4. Discussion oF EXPERIMENTAL RESULTS. 
4.1 Constancy of the Size Frequency Distribution. 

The experiments of Bernardini, Cortini and Mantredini (1949 a) 
showed. that the size frequency distribution of the stars observed in 
Ilford C2 emulsions varied little between sea level and 20,000 m., which 
suggests that the energy spectrum of N-rays has the same shape at all 
depths in the lower atmosphere. To compare these results with those 
from the present experiment, the numbers of stars with different numbers 
of heavy tracks are collected in Table IV. for both sets of plates. The 
numbers are normalized to 100 for Nyj=5. 

Comparing row (i.) with row (ii.), it will be seen that the frequency 
distribution for different values of Nj, is the same above and below 
30 ems. of lead. The results of Bernardini et al. (1950) show, however, 
that there is a much larger proportion of 3- and 4-prong stars under 
2 ems. of lead; under 30 cms. of lead this transition effect vanishes, 
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“The results of Malaspina et al. (1950) showed a higher rate of stars in 
plates placed on top of a lead absorber than in plates not placed on lead, 
suggesting that there is a larger backward emission of secondary neutrons 
capable of producing 3- and 4-prong stars than in an equivalent thickness 
of air. The difference is probably connected with the larger size of stars 
in lead than in air. Under large thicknesses of lead the excess 
backward flux of star primaries will combine with the secondary N-rays 
emitted in the forward direction in the lead above the plates, giving 
the same form for the spectrum of N-rays as in air. In the present 
experiment the plates were placed on concrete, so that there is a deficiency : 
of star primaries incident on the plates from the concrete. Comparison 
with the results of Bernardini et al. (1950) under 2 cms. of lead shows that 
the decrease should be equivalent to about 6 per cent of all stars. There 
is, however, a compensating effect which increases the proportion of 
3- and 4-prong stars, namely, an increase in o meson stars due to the 


TABLE LV. 


Size Frequency Distribution of Stars. 


Numbers of stars (normalized to 100 for Ny=5) 


ey 3 4 5 6 tf 8 9 210 26 : 
Absorber = 


(i.) light roof 251410 1684+8 10046 5244 3144 3044 1845 3213 162+9 


(ii.) 30 em. lead 245410 16847 10046 5345 3144 1743 1643 3543 15248 || 


(iii.) 2cm. lead. (After 
Bernardini eé al. 1950)} 3604-20 244+20 100+6 150-5 | 


fact that in lead 7 mesons have a greater probability of stopping before 
decaying into » mesons than in air. Two-prong stars produced by 
o mesons are included in the present analysis. The number of such 
stars in the air plates depends on the proximity of dense absorbers. 
Various estimates place the proportion of o meson stars with two or 
more evaporation prongs to be about 3 per cent of all stars. According 
to Bernardini e¢ al. (1949 b), the proportion of such stars increases by 
a factor (2:6+-0-7), under large thicknesses of lead, making an overall 
increase of about 5 per cent in all stars. The increase occurs mainly 
in stars of 3 and 4 prongs. This effect tends to compensate the 
transition effect and accounts for the constancy of the size frequency 
distribution. 
4.2 The Attenuation Length for all Stars. 
The attenuation. length in lead for all stars is : 


A,=305-+7 gm./em.2, 


It is interesting to compare this value with that obtained under other 
absorbers by various workers. 
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It has been suggested that the values are proportional to Al, suggesting 
that the cross-section for star production is proportional to the 
geometrical cross-section. As pointed out by Ritson (1950), however, 
the values given in Table V. are not very consistent with an A law. 

The picture of nucleon cascades allows us to interpret the results in 
a different way. Since the B-component is to a large extent secondary 
to the A-component the attenuation length of the primaries of evaporation 
stars depends on the attenuation length of the parent penetrating-shower 
primaries. This dependence is not complete under 30 cm. of lead since 
a large number of the N-rays present are due to the B-component 
incident from the atmosphere. Moreover, since the nuclear interactions 
of the A-component are included in the analysis, A, for all stars depends 
to a large extent on A, for penetrating showers (n,>2). The interaction 


TABLE V. 


Attenuation Length in Different Materials (in gm./cm.?). 


Material 


Authors 


Carbon Ice Aluminium Lead _ Air 
Barton, George and Jason 
(1951) 166+7 
George and Jason (1949) 310+20 150 
Bernardini et al. (1949 a) 220 300 +20 139 
Harding, Lattimore, Liand 
Perkins (1949) 200 +10 


mean-free-paths of the primaries of penetrating showers have been 
determined for lead and carbon by Cocconi (1949 a, b) and by Walker 
et al. (1950) using counter arrangements. The value for lead is 160 gm./cm.? 
and is consistent with a cross-section equal to the geometrical cross-section. 
For carbon Walker et al. find (81--5) gm./cm.? which is longer than the 
interaction mean-free-path corresponding to the geometrical cross-section, 
suggesting that there is some transparency of the carbon nucleus. 
We see that 
mo Ed 
ealeo. 2 
Aso 
Aa(Po) 
for penetrating showers. This is the same ratio as the attenuation length 
for all stars. Thus the ratio of the values of A, for all stars under lead 
and carbon can be understood in terms of the absorption of penetrating- 
shower primaries, with the B-rays keeping in equilibrium with the A-rays. 
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This is not completely true since some B-rays are incident from the 
atmosphere but due to the transparency for penetrating showers a need. 
not necessarily be proportional to At’. Furthermore, approximately 
one-quarter of all stars are due to B-rays emitted in the backward direction. 
The number of these depends on the nature of the absorber beneath the 
plate so that the value for the attenuation length of all stars varies with 
experimental conditions. 

Another interesting fact, pointed out by Barton, George and Jason 
(1951) is that for all stars A, (carbon) is greater than A, (air) even though 
the atomic weight of carbon is less than the atomic weight of air. Barton 
et al. suggested that the difference was due to the production of stars 
by 7 mesons since for carbon the mean-free-path for star production is 
smaller than the mean-free-path for 7+ decay whereas in air the mean- 
free-path for 7p decay is shorter. Therefore the attenuation length for 
all stars in dense absorbers depends on the local production of 7 mesons. 


4.3. Comparison of A, and x, for the Primaries of 
Penetrating Showers. 


The interaction mean-free-paths of the primaries of penetrating showers 
in lead has been determined by Cocconi (1949 a, b), Walker (1950) and 
Sitte (1950) from experiments with counters. These workers found 
\;~160 gm./em.?, which is consistent with a cross-section equal to the 
geometrical cross-section of the lead nucleus. Walker (1950) found that 
A; increased with decreasing number of relativistic secondaries required 
to trigger his penetrating shower set, for example, A; was (150-++8) gm./cem.? 
for n,>5 and 208 gm./cm.2 for n,=3. . 

Penetrating showers in the photographic emulsion will be defined as 
events with n,>2, that is events of the types 2n+2p+3n+3p+.... 
The attenuation length is determined by observations made with the 
photographic emulsion and the value obtained is Ay=(380-+-65) gm./cm.?. 
The comparable value from Walker’s experiment for A, is probably 
~208 gm./cm.?. It is seen that the attenuation length is longer than the 
interaction mean-free-path. The probability that A, is less than or 
equal to 208 gm./cem.? is less than 0-5 per cent so that the increase can 
be taken as significant, which means that on the average penetrating-shower 
primaries give rise to more than one penetrating shower in lead. The 
experimental value of A,/A; is (1:85-L0-3) for events with n,>2. 

The value of the interaction mean-free-path in carbon is rather uncertain, 
because of the inconsistency in the experimental results, so that A, for 
carbon cannot be compared with the attenuation length for penetrating- 
shower primaries in air. It must be borne in mind that the value of 
A,/A; for lead may be larger than the corresponding ratio for air because 
of the contribution from 7 mesons produced locally in the lead. On the 
other hand the ratio may be reduced because of the larger probability 
for successive nucleon collisions within the same nucleus, with a 
consequent larger energy loss per interaction. 
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4.4, Stars Produced by Fast 7 Mesons under Lead. 

In the present experiment it is possible to observe stars produced 
directly by fast + mesons as they should lead to an increase under lead in 
the proportion of stars with charged relativistic primaries. The 
attenuation length in lead for events of the type Op+1p+2p-+... is 
(405+31) gm./em.?, and for events of the type 1lv+2n+3n+... is 
(260+-34) gm./em.?. In the lower atmosphere these lengths should be 
approximately equal. The value of A, for events of the type 0p-++-1p+2p 
is larger than the value for events of the type 1n+2n-+... by almost 
exactly three times the standard deviation of their difference so that it 
is safe to conclude that some stars under lead are due to fast 7 mesons 
generated locally in lead. From Tables I. and II. it can be seen that the 
ratio of Op+-lp+2p-+. .. events to all stars increases from (13-1-L0-9) per 
cent for the air plates to (17-0-L1-0) per cent for the plates under 30 cm. 
of lead, giving an increase equivalent to (3-9--1-3) per cent of all stars. 

Twenty-six of the fast charged primaries in the lead plates were suitable 
for identification from measurements of their multiple Coulomb scattering 
and grain density (cf. Fowler 1950), and of these eight were shown to be 
less massive than protons. This result is consistent with the estimated 
proportion of such events which are due to fast 7 mesons. 

The number of stars due to fast + mesons under 30 cm. lead is 
approximately equal to the number of o meson stars produced by 
7m mesons at rest. The two are comparable because the probability of 
a 7m meson stopping in the emulsion is of the same order as the probability 
of a fast 7 meson passing through the emulsion giving rise to a nuclear 
interaction (cf. Bernardini, Cortini and Manfredini (1949 b). 

It is possible that some of the events of the type lp, 2p, 3p, ete., in 
the lead plates are the same as the events observed by George and Evans 
(1950) at 60m. water equivalent underground. According to these 
authors the rate of such events at the altitude of the Pic du Midi should 
be about 0-03 stars/c.c./day, that is, about 1 per cent of all stars in the 
lead plates or 12 per cent of 1p+2p+3p-+... events. It is not possible 
to confirm the presence of such events due to » mesons because of poor 
statistics and because such events may be produced by 7 mesons generated 
locally. 

Since the observed proportion of stars due to 7 mesons in the lead 
plates is (3-9-.1-3) per cent of allstars probably less than 6 per cent of all 
stars are produced directly by 7 mesons generated locally. This proportion 
in itself is not sufficient to account for the difference between A, (air) 
and A, (carbon) observed. by Barton et al. (1951). However, the 7 mesons 
generated locally may give rise to energetic secondary nucleons in nuclear 
interactions in the absorber above the plates thus increasing the 
proportion of stars associated with the 7 mesons generated locally. 

From Tables I. and II. the ratio of charged to uncharged primaries of 
penetrating showers (7. e. events of the types n,>2) is (0-86+-0-21) for 
the air plates and (1:67--0-37) for the plates exposed under 30 cm. of lead. 


866 W. G. V. Rosser and M. W. Swift on the 


For nucleon cascades in air the ratio should remain constant and close 
to unity at all depths in the lower atmosphere. However, under lead 
the ratio increases to 1-67-L0-37). Our statistics alone are not good 
enough to establish the increase with certainty, but a similar experiment 
carried out at Bristol leads to a similar result. During a preliminary 
analysis the Bristol workers observed 84 penetrating showers under 30 cm. 
of lead and of these showers 58 appeared to have charged primaries and 
26 neutral. Combining the two sets of resuits the ratio of charged to 
uncharged penetrating shower primaries under 30 cm. of lead is (1-94-- 0-32). 
This figure is greater than one with a probability of more than 99 per cent. 

It is possible that part of the increase may be due to events of the type 
observed by George and Evans (1950) who found that four of the ten 
stars produced at 60m. water equivalent underground by relativistic 
primaries had more than one relativistic secondary. This would correspond 
to an increase of about 15 per cent in the ratio of charged to uncharged 
penetrating showers but is not sufficient to double the ratio. ; 

A more likely explanation is the production of events of the type 
2p+3p+... by mmesons generated locally in the lead. One possible 
scheme for the production of events of the type 2p is the elastic collision 
of a 7 meson with a proton in a nucleus, both the 7 meson and proton 
emerging from the nucleus with relativistic velocities. ‘Such a process 
is energetically possible for a primary 7 meson of energy about 2 x 10° eV. 
The probability for the process depends on the cross-section for the elastic 
collision of 7 mesons with protons, and since this cross-section is not yet 
known the contribution cannot be estimated. Two further possible 
theoretical schemes for the production of penetrating showers by 7 mesons 
have been suggested by Heitler and Janossy (1950). In the first scheme 
it is suggested that az meson might give a large fraction of its energy to a 
nucleon inside a nucleus and the nucleon then give rise to a penetrating 
shower. According to Camerini et al. (1950 b) a nucleon needs on the 
average an energy of 4 10° eV. to give rise to an event of the type n,>2, 
so that the 7 meson would need more than this energy. Now it has already 
been shown that approximately 4 per cent of the stars under 30 em. of lead 
are due to fast 7 mesons. Less than 10 per cent of these 7 mesons would, 
however, have energies greater than 410%eV. These 7 mesons could 
contribute to a small increase in the ratio of the charged to uncharged 
penetrating-shower primaries; the maximum increase due to this effect is 
estimated to be 20 per cent, which is hardly sufficient to account for the 
observed ratio of (1-94-+0-32). In the second scheme suggested by Heitler 
and. J. anossy 7 Mesons May give rise to penetrating showers by a multiple 
process in which the primary 7 meson splits up into several mesons. 
These authors suggest that if the energy of a 7 meson is about 10° eV. 
it might split into two or three secondary mesons which would be observed 
as events of the type 2p and 3p in electron sensitive emulsions. Approxi- 
mately 30 per cent of the 7 mesons under 30cm. of lead have energies 
greater than 10° eV. so that if the probability for such a process were 
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high it could lead to an appreciable increase in the ratio of charged to 
uncharged penetrating-shower primaries. 

At present no definite conclusions can be drawn about the nature of 
these events, but it is suggested tentatively that they are mainly 
produced by fast 7 mesons generated locally in the lead. 
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SUMMARY 


The yield point of a plastic-rigid body is defined as the moment when 
deformation first becomes possible as the load is increased. The practical 
significance of the yield-point load for an actual plastic-elastic body is 
discussed. It is then shown that part of the plastic zone at the yield 
point, namely the part where deformation occurs, depends only on the 
current surface tractions and not on the previous loading programme. By 
means of the principle of maximum plastic work and its complimentary 
minimum principle (Hill 1950) approximations from above and below to 
the yield-point load can be obtained ; some examples are given. Recent 
American work on the plastic limit design of structures is critically 
reviewed and shown, when properly regarded, to be a particular application 
of the principle of maximum plastic work. 


§ 1. THe DEFINITION AND SIGNIFICANCE OF THE YIELD Pornt. 


THE following considerations relate to an ideal work-hardening plastic- 
rigid material in which Young’s modulus is indefinitely large. Consider a 
mass of such material subjected to some loading sequence. The moment 
when it first becomes locally plastic is usually known as the elastic limit, 
while the moment when it begins to deform will here be called the yield 
point (the term collapse, which has also been suggested, is suitable only 
when the subsequent deformation would render the body useless). When 
the elastic limit and the yield point are not coincident, the body remains 
entirely rigid in the corresponding loading interval, even though partly 
plastic. The calculation of the stress, state and development of the 
plastic zone while the whole body is rigid have already been discussed 
by the writer (1950 a, pp. 237-245). Here, it is proposed to examine 
certain properties of the stress state at the yield point itself which depend 
only on the current surface tractions and not on the previous loading 
history. 

In the absence of passive constraints the rate at which the loads need to 
be increased to enforce further distortion following the yield point is 
governed mainly by the rate of strain-hardening—though also by the 
changes of shape (so that even if the hardening is zero the load does not 
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generally remain constant during a finite distortiont). When the rate 
of hardening is of the order of the yield stress, only a slight increase in 
load is sufficient to produce strains of a magnitude which would, in an 
actual structure, be regarded as amounting to callapse. 

Correspondingly, in a slowly hardening plastic-elastic body such strains 
generally occur under loads differing little from the plastic-rigid yield- 
point loads. The previous permanent deformation, occurring since the 
elastic limit, is of an elastic order of magnitude, due to the constraint of 
the non-plastic material ; an overall load-displacement curve has therefore 
a pronounced bend in a small loading interval in the neighbourhood of the 
plastic-rigid collapse load. (An obvious exception is when the plastic 
zone is entirely contained within an elastic zone which accommodates the 
severe local distortions, as when pressure is applied to a small cavity). 
In this sense, the plastic-rigid yield-point loads have a significance for an 
actual (plastic-elastic) body in which the definition of a precise moment of 
collapse, or yielding, must be somewhat arbitrary. More exactly, if a 
non-hardening plastic-elastic body is monotonically strained, the requisite 
loads approach asymptotically (though not necessarily steadily) to the 
plastic-rigid yield-point loads when changes of shape are ignored in the 
analysis. The asymptotic values are generally not reached in any 
finite strain (the contrary is often wrongly stated), but the actual loads 
differ from them only by an insignificant amount after a permanent strain 
which is still of an elastic order of magnitude. 

The main purpose of the present paper is to examine the factors 
governing the distribution of stress in a plastic-rigid body at the 
yield point, and to draw attention to some rather surprising properties 
which have not been observed before. 

These properties are intimately related to certain extremum theorems for 
a plastic-rigid body, namely the principle of maximum plastic work and its 
complementary minimum principle established by the writer (1948 and 
1950 b). These principles refer to the rate of working of forces on the 
surface of a plastic zone, at all points of which deformation is occurring 
(the zone does not necessarily extend over the whole body). Since any 
equilibrium system of forces acting on a rigid zone does no work, the 
extremum principles apply also to the work of the forces on a body which is 
only partly plastic. These theorems supply a means of obtaining 
approximations, from above and below, to the actual rate of work, and 
also, in certain situations, to the yield-point loads themselves. However, 
by contrast to the analogous methods for a perfectly elastic body, they do 
not furnish approximations to the whole internal state of stress in a partly 
plastic body at the yield point since this depends, in part, on the previous 
loading programme. It is only for particular surface tractions, such that 
deformation can occur under them at all points, that the whole internal 


a ei Jf re eae eta pee. ta te 

+ For this reason the term “ flow limit ”’, suggested in a recent paper by 
Pell and Prager, is unsatisfactory, as “flow” signifies deformation under 
constant stress, 
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state of stress is unique (this was first shown by the writer in discussing the 
problem of a uniformly bent and twisted beam ; 1948). 

Recently, approximation methods, identical in principle, have been 
applied by American writers on limit design of structures (e. q. Greenberg 
and Prager 1951). Owing, presumably, to the simplifications involved in 
the conventional idealization of a structure, it has escaped notice that the 
relevant theorems are special cases of the maximum work principle and its 
complement. This is even unremarked by Pell and Prager (1951) who 
treat the continuum problem of the yield point of a completely plastic plate 
under uniform loading, and also by Drucker, Greenberg, and Prager (1951)f. 
To the latter is due the first illustration of the usefulness of the methods 
when the body is partly plastic, viz. that of the collapse of a pressure vessel. 
A further illustration has been given by the writer in the problem of the 
indentation of a thin strip by a pair of opposed dies (described by A. P. 
Green, 1951, who has obtained the exact solution). Other applications 
are outlined below. 


§ 2. THe PLASTIC STRESS-STATE AT THE YIELD POINT. 


‘Let the criterion of yielding of an element be 
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in some given state of hardening. The rate of strain under continued 
loading is taken as 
of 


PDO . . . . . . . . . (2) 


where A is a positive scalar factor of proportionality. This relation 
between the ratios of the components of the stress and strain-rate tensors is 
suggested by experiment, and has also been derived theoretically for an 
aggregate of crystals which deform by slip according to the Schmid law 
(Bishop and Hill 1951). 

Ifthe criterion (1) be regarded as a surface in stress hyperspace, equation 
(2) states that the vector representation of €,,; is parallel to the outward 
normal at the stress point o;;. If the yield surface is everywhere concave 
to the origin (so that the relation between the directions of the stress and 
strain-rate vectors is one-one), it is evident that 
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where oj; is any stress within or on the surface (the summation convention 
being adopted.) This is the maximum work principle for an element ; it 
states that the actual stress does work at a rate not less than any other 
stress compatible with the yield condition. 
——$<—<—— ee 

} The work of all these writers is marred by erroneous statements that the 
theorems apply to a plastic-elastic material ; furthermore, collapse is defined 
as deformation under constant load—a situation which generally does not 
occur, as mentioned previously. 
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Consider, now, the yield point of a body under given forces F per unit 
area over a part Sy of the surface 8, and given velocities u; over the 
remainder Sy (more generally, certain components of F, and the comple- 
mentary components of uw; may be given over the same part of the surface). 
Consider the set of equilibrium states of stress, and associated velocity 
distributions, corresponding to the same distribution of hardening and 
_ the same given surface conditions. The stress must satisfy (1) and (2) 
wherever the respective strain-rates do not vanish, and must not violate 
the yield condition elsewhere. Included among this set are, of course, 
actual yield-point stress states and their associated distortion modes ; 
such actual states may be considered to be reached from different initial 
states of residual stress or by different loading programmes. Let (o 4j, W;) 
and (o;;, u;) be any two members of the set. Across certain internal 
surfaces some of the components of stress may be discontinuous ; for 
equilibrium, only the stress resultants need be continuous. Across other 
internal surfaces the velocity may be tangentially discontinuous ; such 
a surface of discontinuity has to be regarded as the limit of a narrow layer 
through which the velocity component parallel to the surface changes 
rapidly. The state of stress at each point of a discontinuity is therefore 
that corresponding to a simple shear strain-rate in the direction of the 
relative velocity vector. The material being assumed isotropic, the stress 
is thus a pure shear, k say, along the discontinuity, together with some 
hydrostatic component ; since fis concave to the origin and is assumed 
not to involve the hydrostatic component of stress, no shearing stress 
component can exceed k& in a given state of hardening. Let 7* be the 
shear component of oj, along a discontinuity of u,;, and 7 be the shear 
component of o,; along a discontinuity of u;. Then, by virtual work, 

§(F FF) (u — uz) AS= J(o 9%) (€ gs —) EV + [(h—7")[u] ASp 
+J(k—7)[u*]dSp, . - - (4) 
where the square brackets denote the absolute amount of the 
discontinuities across the respective surfaces S, and 85. The left-hand 
side vanishes identically since one or other of the factors is zero at each 
point of the surface of the body. According to (3) and the fact that 7 and 
7* cannot exceed k numerically, all terms on the right-hand side are non- 
negative. Hence they must each be zero. This requires that the 
deviatoric parts of o;; and oj, are equal on Sp and Sj, and wherever ¢,; 

or €;; is not zero. 

In other words, stress states of the set considered, and in particular 
actual yield-point states, are identical wherever deformation is occurring 
and can differ only in the common rigid zone. In particular, in plane 
strain, the slip-line field in the plastically-deforming zone is uniquely 
determined by the given surface conditions and internal hardening at the 
yield point. On the other hand, the present analysis says nothing about 
the possible modes of distortion (except in so far as they must be compatible 
with the stress). This theorem justifies the procedure adopted in the past 
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for calculating yield-point loads : namely, to ignore the development of the 
plastic zone and to find some stress and. velocity distribution satisfying 
boundary conditions at the yield point only. It is necessary also to show 
that there is an associated stress distribution in the rigid zone, not violating 
the yield condition ; in so far as this is left unexamined, previous solutions 
are not completely rigorous, as the writer has frequently emphasized. If 
the yield condition is in fact locally exceeded in the rigid zone by a certain 
amount, such a solution is valid for a material which is at least that much 
harder there than the one under consideration. It is intuitively obvious, 
and an immediate consequence of the maximum work principle (5) below, 
that the yield-point load will then tend to be overestimated. 

An example of the use of this theorem may be instanced. Consider the 
slip-line field and associated velocity solution proposed by Prandtl to 
satisfy the surface conditions at the yield point of a semi-infinite block 
indented by a flat die. Ifit could be proved (as seems probable) that there 
exists an associated stress distribution in the rigid zone, not violating the 
yield condition, then it would follow that the actual plastic zone certainly 
extends at least as far as Prandtl’s slip-line field. Similar conclusions hold for 
the notched bar problem (Hill 1950 a, pp. 249-251). It would not necessarily 
follow, however, that the mode of distortion suggested by Prandtl is the 
actual one for an initially stress-free and uniformly hardened block. To 
prove that an assumed mode is correct, it is essential, as remarked by the 
writer (1950 c¢), to examine the stress increment occurring during an 
increment of distortion and to show that all elements assumed to deform. 
have in fact remained plastically stressed. (There is no possibility of 
overlooking this point when dealing with a plastic-elastic material since the 
stress increment in any case enters the relations between the ratios of the 
components of stress and total strain-rate). This applies equally to the 
alternative mode suggested by the writer for this problem (op. cit.). To 
what extent a yield-point mode is independent of loading history is 
a matter for future research. 


§ 3. ExTREMUM PRINCIPLES. 


Let oj, now refer to any equilibrium state that does not violate the 
yield condition at any point where ¢ ;;40 (oj; is not necessarily any actual 
yield-point state). Then, as in (4), 

J(F ,—Fj)u,; dSy=J(o 5;—ofj)e ; V+ J(k—7*) [u] dSp>0 . . . (5) 
if KF; =F; on Sp. This is the maximum work principle for a finite body 
(Hill 1948). 

In particular, if F7=F; all over the surface the left-hand side of (5) 
vanishes. Since each term on the right-hand side is non-negative, both 
must be zero and oj, must equal o,; wherever deformation is occurring 
(including, of course, velocity discontinuities.) Thus, if given forces act on 
the surface of a body in a given state of hardening and if it is possible to 
find some corresponding internal equilibrium distribution of stress oi; and 
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velocity w;, satisfying (1) and (2) where €,;40 (with positive A), then, 
wherever ¢,;40, the body must certainly be plastic and this o;; must be 
the actual stress. 

For the complementary minimum principle consider any velocity distri- 
bution uv; having zero divergence and such that ut=w,on Sy. Then, by 
virtual work, ° 

fF ,u; dS= fo ,, |, dV-+Jr[u*] d8*. 


A] 


Let oj; be the stress that would be needed to produce a strain-rate ei 
in the given state of hardening ; the distribution oj, is not necessarily in 
equilibrium and is undefined where <7,=0. Then, by the maximum work 


inequality (3) and the fact that |r|<k, 
Fw; dSy <Joj,e;, dV+ fh[u*] dSi—JF;ufdSp, . . . (6) 


i) 


which is complementary to (5). This result is more general than 
one obtained previously by the writer (1950 b), which was valid only for 
the yield criterion of von Mises. For the latter, | oj;|= /(o{0/;)=const., 
where oj is the stress deviator, and oj¢;;—|o;||¢;|. In plane 


strain, oje;,—ky* whatever the function f, where y* is the maximum 


shear strain-rate (twice the tensor component). 


§ 4. APPROXIMATE ESTIMATION OF YIELD Pornt Loaps. 


According to (5) a lower bound to the rate of work by the forces over Sy 
is to be obtained by selecting any equilibrium distribution of stress not 
violating the yield condition. An upper bound, on the other hand, is to 
be obtained from (6) by means of any distribution of velocity. With a 
well-judged choice these bounds can be brought close enough to give a 
useful approximation to the actual rate of work on Sy. 

When the prescribed velocity on Sy is uniform in magnitude and 
direction, a knowledge of the rate of work is equivalent to a knowledge of 
the load on Sy in that direction (and similarly when a component of velo- 
city and the complementary component of force are given). This situation 
arises, for example, when a single load is applied to the body through a 
rigid plane die (as in indentation) or through a rigid part of the body 
itself (as in a notched-bar test). Similarly, when the velocity distribution 
on Sy amounts to a rigid-body rotation or twist, the bending or twisting 
couple is immediately calculable from the rate of work. 

A slightly different problem presents itself when several loads are 
applied, their ratios at the yield point being given ; this includes pro- 
portional loading, to which writers on limit design have restricted their 
attention. The single parameter specifying the absolute magnitudes of 
the loads is derivable from the extremum principles. To see this, let Sp 
now denote the part of the surface where the load-ratios 7; are given, 
S,, standing as before for the part where absolute values of certain loads 
are prescribed. In addition there may be fixed frictionless constraints, 


é 


but specified non-zero velocities can no longer be admitted since the 
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corresponding forces would appear as further unknowns. If F,=Fr; 
on Sp at the yield point, then, as in (5), 
BSS. 5 le Ree 


yn cancelling the factor fr;w;dSp, where oj, is such that F7=F*r; on 
S, and F*=F,onS,. It follows from (7), as remarked by Greenberg and 
Prager (1951), that the value of F is unique. When the deformation 
zone develops from a single point a limiting argument appears sufficent to 
establish the same result. Also, as in (6), 


Fr; ut dSq<Jotet, dV+fk{u"] dS§—F;ufdSp. . . . (8) 


§ 5. Some APPLICATIONS. 


It is to be noted that the entire discussion and the approximation 
methods are equally relevant for those problems of steady motion where 
there is no hardening and where the external contour of the body is known 
a priori. Steady motion can, in fact, be regarded as the integrated result 
of incremental distortions due to a succession of yield points. So, for that 
matter, can any plastic problem, but since the external configuration is 
continually changing, or is not known beforehand, the present methods 
cannot be applied directly. 

As an example of the estimation of a steady-state load, consider the 
drawing of a sheet in plane strain through a smooth wedge-shaped die 
(without bulging). An intuitively obvious lower bound is provided by the 
stress distribution corresponding to radial flow towards the virtual apex 
of the die: r being the fractional reduction, (5) gives 2k In{1/(1—r)} for 
a lower bound to the drawing stress, as in elementary theories of drawing. 
An upper bound is provided by a velocity distribution in which the 
reduction is effected by block-sliding over a criss-cross of lines, of which the 
extreme pairs pass respectively through the points of entry and exit. If 
the lines make equal angles @ with the axis, there is a constant velocity 
discontinuity across each pair. This is sin «/sin (@—«) at the first, where 
2 is the die angle, and is multiplied at each successive pair by a factor 
sin (0-+-«)/sin (@Q—a). Since the dimensions of successive pairs are multi- 
plied by the inverse factor the work of shearing for each line is the same. 
If there are n pairs (6) gives an upper bound of 


2kn sin « cos a 
sin (@—«) sin (@-+-«) 
where 
sin (@—«) 


(1—rjin 
sin (0-+-«) 


The best approximation, with the corresponding value of 6, is to be obtained 
from the integral value of » making this upper bound least. In the 
practical range of reductions only the values n—1 and 2 need be considered. 
For those reductions not so far covered by an exact theory (see Hill 1950 a, 
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p. 167), the greatest difference between these upper and lower bounds is 
of order 5 per cent. In this range, then, the elementary theory (which 
assumes homogeneous distortion) is a fairly good approximation. 
However, contrary to general belief, there are reasons for thinking that 
the relative amount of redundant work does not tend to zero steadily as the 
reduction increases, but oscillates with a steadily decreasing amplitude. 

Consider, finally, the yield point of a symmetrically-notched bar 
extended in plane strain. Let the width of the bar be w times the smallest 
section. An upper bound is obtained from a deformation consisting of 
block-sliding across a pair of lines from the notch roots, intersecting the 
longitudinal axis in an angle 6, say, and a pair passing from this point to the 
sides of the bar and making 47—@ with the axis. For the stress at the 
ends of the bar equation (6) gives an upper bound k(cot 0-+-w tan 6). The 
smallest value of this is 2k\/w, when cot 6=1/w. This is valid for any 
shape of notch geometrically compatible with this fictitious mode of 
yielding, but is obviously likely to be best for a parallel-sided notch with 
zero root radius (since this should have the highest constraint factor). 
An exact solution (Hill 1950, unpublished) gives a yield-point stress of 
2-05 x 2k for w=4:-7, compared with the upper bound of 2-17 2k. For 
w>8-6 exact calculation gives 2:57 x 2k (Hill 1950 c), while the upper bound 
exceeds this value when w>6-6. We may infer that the constraint factor 
is not more than a few per cent less than 1/w when w is less than about 3. 
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SUMMARY. 


In the binary alloys of Al-Zn there exists, above a certain temperature, 
a mixed crystal region in which the two phases « and «’ both have the 
game face-centred cubic structure but with different compositions and 
lattice spacings. It is shown that an explanation for the formation of 
these phase boundaries is provided by the electronic theory of the elastic 
properties of single crystals of aluminium. A maximum in one of the 
shear constants at an electron concentration corresponding to the 
ata’ region gives rise to a subsidiary maximum in the temperature- 
dependent part of the free energy, which leads, according to the usual 
construction, to a mixed phase region. 


§ 1. INTRODUCTION. 


In a recent paper (Leigh 1951) it has been shown that in aluminium 
solid solutions a maximum is to be expected in the shear constant 
3(C,,—C,,) at an electron concentration in the neighbourhood of 2-67 per 
atom. The elastic constants }(C,,—C,,) and C,, are shown in fig. 1 as 
functions of the concentration of a divalent solute. In the theory of the 
elastic properties of pure aluminium it was found necessary, in order to 
explain the observed values of the shear constants, to assume that there 
are a small number of electrons overlapping the square faces of the 
Brillouin zone. The sharp peak in 4(C,,—C,,) occurs at the point at 
which this overlap disappears. 

It is now shown that the maximum in 4(C,,—C,,) provides an 
explanation for the existence of a mixed crystal region which, below 
355° C., interrupts the range of stability of the primary solid solution 
of zine in aluminium. The aluminium-rich end of the phase diagram 
of the Al-Zn system (Raynor 1944) is shown in fig. 2, with zine 
concentration given in atomic per cent. The mixed crystal region in 
question is shown as «--«’. The phases « and «’ both have the face-centred 
cubic structure, and only differ in composition and lattice spacing, the 
gap in the solubility disappearing above 355°C. The f phase is the 
primary solid solution of aluminium in zinc, and has the hexagonal close- 
packed structure. 
oe a es 

* Communicated by Professor H. Jones, 
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Fig. 1. 
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The existence of the solubility gap must be due to a subsidiary maximum 
in the free energy of the solution regarded as a function of zinc 
concentration. The positions of the phase boundaries are then given 
by the points of contact of the common tangent to the parts of the free 
energy curve on either side of the maximum. The present theory shows 
that the calculated maximum in }(C,,—C,,) leads to a subsidiary 
maximum of the required form in the contribution to the free energy 
from the thermal vibrations of the lattice. This provides some 
confirmation of the hypothesis of overlap over the square faces in pure 
aluminium. 


§2. THe FREE ENERGY. 


The free energy of a solid solution can be expressed as a sum of 
contributions from the electronic energy, the mixing entropy and the 
thermal vibrations. 

The elastic constants of the solution, as shown in fig. 1, have been 
obtained by neglecting changes in the density of states curve as zine 
is taken into solution. To this approximation the first derivative of the 
Fermi energy with respect to electron concentration cannot have a 
discontinuity anywhere ; and the second derivative, given by 1/N(Q), 
where N(¢) is the density of states at the Fermi surface, cannot be 
anywhere negative. In the present theory, therefore, no subsidiary 
maximum in the free energy can arise directly from the Fermi energy. 
(Compare the discussion by Hume—Rothery (1944), who also postulated 
overlap over the square faces in aluminium.) It will be assumed that the 
whole electronic energy varies smoothly and approximately linearly 
with zine concentration. It can then be neglected, as it will not affect 
the positions of the phase boundaries. 

The contribution to the free energy per atom from the mixing entropy 
is given, assuming complete disordering, by the expression 


kT [x log x+(1—za) log (1—2) J, 


where T is the absolute temperature and x is the solute concentration. 
The contribution to the free energy per atom from the thermal vibrations 
is given by the expression 


ir {Rhv-EKT log [1— exp (—hv/KT)]}f(0) dy, 


where f(v) dv is the number per atom of normal modes of vibration of the 
lattice having frequencies lying in the interval v to v-+dv. We are here 
concerned primarily with the effect on the frequency spectrum f(v) of 
changes in the elastic constants, rather than with the actual form of the 
spectrum. It is therefore reasonable to neglect the effect on the spectrum 
of the dependence of the wave-velocities on wavelength, since this effect 
may be assumed to remain constant or change smoothly with zinc 
concentration. As a first approximation we use the Debye spectrum 
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for an Snag gle crystal, taking for the value of the rigidity modulus the 
mean of 3(Cy,—Cy2) and C,,. The contribution to the free energy per 
atom is then given by the expression 


9 T \8 (9p/T 
3hOn+94T( 5) fs fog Gl -=er ede) eae ke 41) 
where @,, is the Debye characteristic temperature, given by the relation 
hf 3 \% 
p= k (a) rs 
Q is the atomic volume, and c is a mean wave velocity given in terms of 
the density p, the bulk modulus K, and the rigidity modulus p, by the 


equation 
5 p 3/2 p\3/2 
one (es) a e) 
Fig. 3. 
40 
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Zinc concentration 
The Debye characteristic temperature as a function of zine concentration. 
The thin curve shows how @p would vary if the elastic constants remained 
constant. 


In calculating @) as a function of zine concentration from the above 
formulz, we can neglect the variation of 2 and K. The change in 2 is 
small, and it is easy to see that the effect on ¢ of any change in K is small. 
Moreover, there is no reason to suppose that they do not vary smoothly 
with zinc concentration, and the effect on the positions of the phase 
boundaries would therefore be very small. In calculating c, K has been 
taken constant and equal to its value for pure aluminium, which is 
7-78 1011 dynes cm.~? (Lazarus 1949). The large variation of p with 
zine concentration, arising from the differing atomic weights of zine 
and aluminium, has been taken into account. 

Figs. 3 and 4 show the variation with zinc concentration of the Debye 
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temperature ©, and the sum of the mixing entropy and Debye 
contributions to the free energy. The zero of F has been fixed arbitrarily 
at the value for pure aluminium. In each case the thin curves show the 
variation due to p alone, with p constant at its value for pure aluminium. 
In calculating the free energy, T has been taken equal to 548° K., the 
lowest temperature at which the «-+«’ region exists, though, in fact, the 
choice of the value of T has practically no effect on the calculated phase 
boundaries. The common tangent to the two parts of the curve is shown, 
the points of contact being at approximately 20 and 70 atomic per cent 
of zinc. The observed phase boundaries at 275°C. are at 16-0 and 
59-4 atomic per cent. The positions of the calculated phase boundaries 
obviously depend almost entirely on the position of the peak in 3(C,,;—C9), 
and, in view of the approximations involved in the calculation of the 
elastic constants, the agreement is surprisingly good. 


Fig. 4. 


_ -O8 


x—> 
Zinc concentration 
The free energy (Debye and mixing entropy contributions) as a function of 
. > . ne = r 4 5 
zinc concentration. The thin curve shows how the free energy would 
vary if the elastic constants remained constant. 


It is perhaps worth noticing that the qualitative result of this caleu- 
lation, namely that the maximum in }(C,,—C,,) leads to a subsidiary 
maximum in the free energy, does not depend on the use of the Debye 
approximation to the frequency spectrum. Apart from neglecting the 
dependence of the wave velocities on wavelength, this approximation 
involves (i.) neglecting the dependence of the velocities on direction, and 
(ii.) allowing for the fact that the number of degrees of freedom of the 
crystal is limited, by the simple device of cutting off the total spectrum 
of the continuum (i.e. the sum of the contributions from the three 
frequency branches) at a single maximum frequency. A more. general 
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form for the contribution to the free energy per atom, in which only the 


dependence of the wave velocities on wavelength is neglected, is given 
by the expression 


1 (47 Q\ 3 9 T\3 p@rxT 

3 (=) ie | (ons E kOt-+9KT( Ge) i : log (1 —e*)é ae |e 3 
: (2) 
in which the free energy is expressed as a sum of contributions from the 
three frequency branches, each contribution involving an integration 
over direction. dw is an element of solid angle, and o* is the magnitude 
of the maximum wave vector in the direction of dw. In general the wave 
vector is restricted to lie within a polyhedron. (If the polyhedron is 
replaced by a sphere of equal volume, o* is constant and equal to 
(3/477.Q)!), @; is a temperature corresponding to the maximum frequency 
in the direction of dw in the ith branch, and is given by the relation 

= z, o'C;, 

where c, is the corresponding continuum wave velocity in that direction. 
The dependence of the ¢,; on direction is very complicated, but the extreme 
values of the velocity in each branch are simple functions of the elastic 
constants. The following table gives the values of the c; in the principal 
sets of crystallographic directions (see, for example, Blackman 1938). 


(100) (110) (111) 
pc? Ci Cu pos sac: 3 Ki een 
pcs Cag (Cy, —Cy2) $ (Cy —Cy2+ Cay) 
pce CW 4(C,, +0 ,.+2C,,) $ (Cy, +2C0,.+4C44) 


These values include the extreme values in each branch. In fig. 1 these 
values of pc? are shown as functions of zine concentration. The bulk 
modulus K, which is equal to 4(C,,+2C,,), has been taken constant 
at its value for pure aluminium as before. The effect of an increase in 
1(C,—C,,) on a value of c; corresponding to a given direction is to increase 
it or leave it unchanged, but never to decrease it. The peak in 3(C,,—Cjp), 
regarded as a function of zine concentration, therefore gives rise to 
similar peaks in the c; in every direction, the size of the peak being very 
large in the neighbourhood of some directions, e.g. the (100) directions 
in the third branch, and vanishing in the neighbourhood of some others, 
e.g. the (100) directions in the second branch. The same is therefore 
true of the @;, and, since the expression in square brackets in (2) is an 
increasing function of @;, it follows that there must be a subsidiary 
maximum in the free energy. 
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“The present theory does not account for the disappearance of the 
solubility gap at higher temperatures. This is not surprising, as the 
theory of the elastic constants only refers to absolute zero. It might be 
thought that the explanation lies simply in the direct effect on the 
Fermi energy contributions to the elastic constants of the finite tem- 
perature form of the Fermi distribution, this leading to a smoothing off 
of the sharp peak in }(C,,—C,,). Calculation shows, however, that this 
is not so, the effect on $(C,,;—C,,.) being only to shift the peak a little 
in the direction of lower electron concentration. 
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SUMMARY. 


The concept of a structural time-constant is introduced by considering 
the time auto-correlation function of any waveform, and a structural 
frequency-constant is similarly defined in terms of the spectral auto- 
correlation function. By setting up a combined time and frequency 
auto-correlation function x, an absolute time-frequency constant is 
obtained. Its value is always unity, and it is suggested that this expresses 
time and frequency uncertainty more precisely than the relation 
dt. dfS1/47. To illustrate the function xy, two special waveforms are 
considered, namely, the complex Gaussian pulse and the real Gaussian 
‘pulse-train. 


§1. INTRODUCTION. 


No waveform ;/(¢) can be confined within a very short time interval and 
yet have a frequency spectrum s(f) which occupies only a very narrow 
band offrequencies. A limit is set by the well-known relation 


5) AY See ee 


which in its wave mechanical form is of course the Heisenberg 
inequality. This is a property of any pair of Fourier transforms ¥(¢) and 
s(f); the quantities 5¢ and 8f are the standard deviations of ¢ and f in 
||? and |s|* respectively. Unfortunately, the relation (1) is an equality 
only when # and s are Gaussian in form. It is therefore tempting to seek 
some stronger expression which, while reducing to the equality (1) for 
Gaussian waveforms, would yet remain as an equality for waveforms 
in general. It might be, for instance, that the weakness of inequality 
arises from using the standard deviation as a measure of spread. This 
line of attack, however, is quickly ruled out by considering frequency 
modulated waveforms. 

It is suggested in this paper that the elusive feature in terms of which 
the uncertainty principle can be exactly formulated, is contained in the 
idea of joint time and frequency auto-correlation. The statement of 
uncertainty which results from this approach becomes a statement, not 
Cel ee ie Sa aoa SD a ec aE EEE 

; * Original version received April 12, 1951. 
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about spreads, but about the lack of orthogonality between the waveform 
itself and the same waveform displaced both in time and frequency. 
Such a time and frequency domain was first discussed by Gabor (1946) 
and has also been treated by Ville (1948), though not perhaps in precisely 
the same way as is done here. 


§ 2. GENERALIZED TIME AND FREQUENCY CONSTANTS. 


The notion of a “time constant ” is familiar to every physicist and 
engineer, but it is not usually defined precisely except for a simple 
exponential transient. One of the simplest linear electric circuits, for 
example, is the one whose impulsive response is given by 


0, t<0, (2) 
y= herein. bOe 
mf ie exp (—t/t)) #t>0, 


and t, is known as the time constant of the circuit, or of the waveform ¢. 
When the input is not a simple impulse the output can be deduced from 
the impulsive response by the principle of superposition : the output 
is simply the convolution or resultant of the input and the impulsive 
response. Rather loosely, it is a degraded form of the input in the sense 
that any time-structure which is fine compared with t) is smoothed out. 
To formulate without explicit recourse to a frequency description what is 
meant by the “‘ time-structure ’’ of any waveform is a problem which 
has received considerable attention in the field of random noise 
fluctuations. Here time-structure may for many purposes be adequately 
described in terms of a temporal auto-correlation function, though the 
formalism of auto-correlation is by no means restricted to the study of 
random functions. For a non-random (and in general complex) waveform 
b(t), the temporal auto-correlation function is given by 


ane fo _ woyrtets a a 


If 4(t) be normalized as in wave mechanics, the function g has the value 
unity at the origin, and it can be shown that |g| never exceeds unity. 
It therefore seems reasonable to define the amount of temporal auto- 
correlation, whether it be positive, negative or complex, by the integral 


el \a(z) |? dr, 1, ey 


since the vanishing of g implies no correlation, or more simply, the 
orthogonality of %(t) and y(t-+-7). The constant 7 may be called the 
structural time-constant of the waveform % without conflicting with 
established use of the unqualified term, since it is readily shown that 
when # is given by equation (2), 7’ is equal to tj. Using the proposed 
definition, the structural time-constant of a Gaussian waveform is given by 
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where d¢ is the measure of spread which occurs in the Heisenberg inequality. 
In a similar manner, we may define a structural frequency-constant for 
any spectrum s(f) by the equation 


F= | |G(¢)P ag, iia teomE es a (6) 


where G(¢) is a frequency auto-correlation function given by 


ad)=| sMPsfted. om 


The relation between b(t) and s(f) being the Fourier transformation 


7 oO 


WO=| _s(f) exp @mift) df, Pe ac yaar: (3) 


f= [WW exp (—2nift) dt par Tris (0) 


it is easy to show that g(7) and G(¢) may be expressed in forms alternative 
to (3) and (7), thus 


gin)= [-_ s()P exp (—2nifr) df eer (10) 


G(d)= { [p(t)P exp (—2midt)dt.  . . . . (12) 
By applying Parseval’s theorem to equations (10) and (11), we may now 
express the structural constants in the forms 


fel ae df eee ede Nees ec £12) 
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By comparing equation (5) with the Heisenberg equality, we see that for 
real Gaussian waveforms 7'F is equal to unity, but this product is certainly 
not invariant for all waveforms, nor should we expect it to be. Since 
|s|? encloses unit area—its normalization resulting from that of | |?—it 
can be seen from equation (12) that 1/7’ measures, in a peculiar way, the 
frequency spread of |s|?. Similarly, 1/F is a measure of the extent to 
which ||? is spread out in time. We might therefore expect that when 
we have 8t. 5/1/47, we could also have TF <1, and it will be shown 
in $4 that certain very simple frequency modulated waveforms do 
satisfy both these conditions. 


§3. Jomnt Trmze-FREQUENCY AUTO-CORRELATION. 


So far, we have considered auto-correlation independently in time and 
frequency, but although there may be no correlation between +(?) and. 
w(t-+7), nor between s(f) and s( f+¢) for particular values of 7 and 4, 
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it can happen that a combined time and frequency displacement (r, ) will 
restore correlation. Let us therefore investigate the joint time-frequency | 
auto-correlation function 


x(7 @)= | MOdAC+t) exp (—2nigt) dt, . . (1) 


suggested by comparison of equations (3) and (11). A little manipulation 
reveals the alternative form 
} 


x(t, =| s*(f)s(f-+4) exp (—2mifr) df, (15) 


which is suggested by comparing equations (7) and (10). This function 
is not entirely a new one. Ville (1948) has considered it in a very similar 
form, and Gabor (1946) has considered it in discrete form, except that 
in his theory, the two functions within the integral are different: one 
of them is a Gaussian analysis function, the other is the waveform under 
consideration. His treatment corresponds to cross-correlation rather 
than auto-correlation. 

The main purpose of the present paper is to point out that the amount 
of auto-correlation in time and frequency combined, may be measured 
by the value of 


c={_ | [xPdrag. PARR Ts. EE 


This is merely a generalization into two dimensions of the one-dimensional 
structural constants 7’ and F. It is an absolute structural constant, because 
it turns out to be invariant for all waveforms. Thus, applying Parseval’s 
theorem to equation (14), we have 


[-_IxPag=[- [HF |or(i+r)P dt 2 2... (17) 


and integration with respect to 7 gives 


fr oO co foe) 
C= | 24h xP dp dr | WOPRd=1, . . . (8) 
since ||? is assumed normalized. This expresses the “ uncertainty 
principle ” of waveform analysis in an exact form. One may verify that 
in wave mechanics, C would equal Planck’s constant. Whether or not it 
proves to be an interesting or useful expression of the uncertainty 
principle depends solely on whether the concept of time-frequency auto- 
correlation is found to have any genuine physical significance. It has 
already arisen in the theory of radar; upon its relevance, if any, to wave 
mechanics, the writer is not qualified to judge. 
It seems worth remarking that, while the functions 


g(7)=x(r, 0), jeg eta ere ye ee Reel) 
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do not provide a complete description of the waveform 7s which generates 
them, there is no loss of information in passing from y to x, a fact which 


may be proved by considering the two-dimensional Fourier transform 
of (7, ¢). 


§4. THe CompLEX GAUSSIAN WAVEFORM. 


It has been remarked that certain “frequency modulated’? waveforms 
are such that 7 <1. Thus, apparently, there is a lack of auto-correlation 
in time and frequency when considered separately, but not in a joint 
_ time-frequency domain, in view of equation (18). The simplest waveform 
which illustrates this is given by } 


%(t)=A exp (—af?), Cprivramseaise nici aateses (21) 


where « is complex with positive real part. This waveform has a Gaussian 
modulus, and its instantaneous frequency, or rate of change of phase, 
is a linear function of time. If a large constant frequency shift were 
included, the real part would repressent a high-frequency pulse enveloping 
a frequency modulated carrier. Since, however, neither a frequency 
shift nor a time shift affects the function y, it will be sufficient to consider 
the simplest form of ~%—as given by equation (21)—in which the 
instantaneous frequency is zero at =0, where the modulus is a maximum. 

When « is purely real, there is of course no frequency modulation, 
and figs. 1 and 2 illustrate typical forms of ||?, |s|?, and |y[?. In 
representing | y|?, a shaded area has been given, its boundary representing 
a single contour. Since all the functions considered here have a Gaussian 
modulus, this representation by means of an ellipse should be adequate. 
In these two examples, it will be understood that 


SEO Litre ame en hee es (22) 
eee Rea Te. (23) 


and as always |x|? is unity at the origin and encloses unit volume. 
Suppose now, that « is complex. It is then possible to obtain both 
a wide ||? and also a wide | s|* as in fig. 3, such that 


Neola ee aes ese (24) 


The ‘‘ bandwidth ”’ 8f is large because of the spread in frequency occurring 
in h(t), which does not show up in |%|?, and therefore does not influence 
St. The structural time and frequency constants, however, are each 


small, such that 
TPM ON Ore Tea OL) 


because of the phase-structure which is masked both in ||? and in |s/?. 

The form of |x|? shows that auto-correlation is distributed in a manner 

which makes mere time and frequency cross-sections quite deceptive. 
Tt can be shown that if we write 


g=a+ib, ° Sheree! Sik wh Son Petes MMR Nasal) (26) 
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the inclination 6 of the major axis of | x|? to the 7-axis is given by 


2nb 
ee tera eae Wren an 2) 
The major axis is thus at 45° when |«|=7. Pure frequency modulation 
is obtained by allowing the real part of « to tend to zero. The 
|x|? function then becomes drawn out into a “knife-edge” of unit 
amplitude and infinite length along a line through the origin which has 
a slope given by 


Cy OaeO meee acs, ees) 


Apart from sign, this slope is simply the rate of change of the instantaneous 
frequency, given by 


d 
2 et eg) ee eee 8, (29) 


The reason for the difference of sign is that |y(7, ¢)|? is unity when a 
time-advance of +, combined with a frequency advance of 4, leaves 
#% unchanged. If ty and dy have the same sign, then the instantaneous 
frequency must decrease with time. 

Knife-edge distributions of |x|? can be obtained mathematically only 
as the result of a limiting process. A knife-edge along the ¢-axis arises 
from a waveform consisting of a single impulse occurring at any time, 
a knife-edge along the 7-axis arises from a purely sinusoidal waveform 
or monochromatic spectrum. Pure linear frequency modulation gives 
a knife-edge at some other angle and is associated with infinite values 
of both the 5f and the 6¢ in the Heisenberg relation. 


§5. THe ReaL Gaussian PuLsE-TRAIN. 


As a final illustration, we may consider a train of real positive Gaussian 
pulses, varying in amplitude according to a wider Gaussian envelope, 
as indicated in fig. 4. It is assumed that the pulses do not overlap 
appreciably and that there are a large number having comparable 
amplitudes. All the following results apply to this approximation. 

The forms of ||? and |s|? are the same: they differ only in their 
parameters. Denoting as previously the standard deviations of ¢ and 
f in |%|? and |s|? by 5¢ and df respectively, it is easily shown that 


Sot ee ee tae. (30) 


Further, let us denote the standard deviation of each individual pulse 
by At and that of each broadened spectral line by Af. Then it can also 
be shown, as is well known, that 

dt. Af= At . sf=1/4n, ee: career (ak) 
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so that the Heisenberg relation is an equality for the wrong measurements ! 
One may further show that the structural time constant 1s 


T= 51. At ee 


where a is the time separation of the pulses, and that the frequency 
constant is 


F=4na6f . Sf, ey i es ee BEY 
where 1/a is the separation of the spectral “lines”. Consequently, as 
will be seen by comparing equations (31), (32) and (33), we have 

TR=1, 5 Ae ee ae 
Fig. 4. 
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Gaussian pulse-train. 


The reason for this is that the time-frequency auto-correlation function 
x is separable as a product 


x(t, ')=9(7)G(d)) oa ee eee oos 
as also in the case of a single real Gaussian pulse. The special characteristic 
of frequency modulation is that y is not separable without a rotation 
of axes. 

Shaded areas for |x|? are shown in fig. 4: the function is made up of 
a lattice of elliptical Gaussian peaks having amplitudes which vary 
according to a Gaussian envelope indicated by the dotted contour. The 
enveloping distribution has the same eccentricity as the cellular peaks, 
since from equation (31) we have 
bt Ls At 
she af’ Cee A ee 


and the area of each unit lattice cell is unity, 


(36) 
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It should be clear from these examples that it is an easy matter to make 
the product of 8¢ and df, which occur in the Heisenberg relation, greatly 
exceed 1/47. Provided, however, that the time-frequency auto-correlation 
function remains separable, the product of the structural time and 
frequency constants remains equal to unity. Whether 7'F can ever 
exceed unity is not yet known, but the absolute structural constant 
defined by equation (16) remains equal to unity invariably. 
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ABSTRACT. 


The beta-spectrum of tritium, H®, has been carefully examined by the 
proportional counter method down to energies ~1 KeV. In the present 
work it has been more thoroughly explored using improved techniques 
designed to permit accurate measurements of the spectrum shape down 
to a lower energy limit ~200 eV. The peo vemcns consisted chiefly 

in employing a counter which was made “ end-free” by means of 
' field-adjusting tubes and the use of an electronic method of measuring 
the pulse distribution in addition to the moving film technique. The 
intensity of the spectrum below 1 KeV. was found to be definitely less 
than theoretically predicted. 


$1. INTRODUCTION. 


THE beta-decay of tritium is an allowed transition. The recent work of 
Curran et al. (1949 a) has shown that between ~1 KeV. and the upper 
energy limit, the spectrum is in good agreement with Fermi theory. 
A tentative suggestion of some deficiency in the number of beta-particles 
in the low energy region was made by these authors. Theoretically, the 
intensity of the spectrum (particles per unit energy interval) at zero 
energy should be about 86 per cent of the intensity at the maximum of 
the distribution. Cooper and Rogers (1950) have constructed a composite 
of the various published spectra and have concluded that all work on the 
spectrum indicates an excellent agreement with Fermi theory from 
5 KeV. to 18 KeV. The introduction of the improvements to the 
technique makes it possible to examine the spectrum shape to lower 
energy limits than in the past without appreciable uncertainties arising 
from spuriously small pulses. These small pulses are those which have 
their origin in the part of the counter near the ends of the wire, since 
in this region the gas multiplication is reduced below the normal value 
unless field-adjusting tubes (Cockroft and Curran 1951) are employed. 


§ 2. APPARATUS. 


A diagram of the counter normally used in this work is seen in 
fig. 1 (a). The field-adjusting tubes were held at the potential appropriate 
to their diameter, thus ensuring a minimum of “ end-effect ’. The value 
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of this potential was checked by ensuring that pulses from a homogeneous 
beam of X-radiation, fired through both windows W, and W,, were of 
equal size at the centre and end of the counter. The eounton was lined 
with a thin sheet of aluminium to minimize secondary effects at the 
cathode. The lining was perforated at the window positions to prevent 
any ionization in the window “ neck ” from being recorded at the central 
wire. 

The pulses from the counter after amplification, were recorded on 
a 35 mm. film. Using this technique, it was found difficult to examine 
the spectrum shape accurately below ~300 eV. The thick stems of the 
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higher energy pulses and the brightness of the baseline itself caused a 
fogging of the film in the very low energy region. To overcome this 
limitation, the pulses were measured using a single-stage pulse analyser. 
This consisted essentially of two scaling units type 1009A. The 
discriminator of one unit was used to control the first stage of both units. 
The insertion of a steady source of potential difference (V volts) between 
the grids of the first stages in each scaling unit supplied a channel width 
of V volts-_the small difference in the triggering potentials of each unit. 
To measure the shape and width of the channel, homogeneous pulses 
were fed into the analyser. Graphs of the difference in the counting 
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rate between the two scalers against the discriminator setting were plotted 
(fig. 2) at various points throughout the range of the discriminator, and 
it was found that the channel width was constant and that the channel 
had “‘ rectangular ” shape. 

The tritium source had argon associated with it to act as a carrier. 
Before the gas was introduced into the counter its temperature was 
‘educed to about —40° C. to remove tritiated water vapour. Measure- 

ents proved that no tritium stuck to the walls of the counter after this 
treatment. 
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§3. Putse Distripution at Low ENERGIES. 


lt will be noted that all the above precautions, both in counter design 
and handling of the source, were introduced to eliminate the possibility 
of affecting the true spectrum shape, particularly in the direction of 
moving higher energy particles towards the lower energy region. 

Three separate measurements were made on the spectrum with the 
counter already described. In the first a small quantity of H3, giving 
a counting rate of ~17,000 ¢./m. (the background counting rate was 
~300 ¢./min.), together with a mixture of methane at a pressure of 
16 em. Hg. and argon at slightly more than 3 atmospheres, was introduced 
into the counter, fig. (a). The spectrum shape between ~200 eV. and 
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2 KeV. was examined, the pulses being recorded on 35 mm. film. The 
energy scale was calibrated by superimposing on the spectrum the 
fluorescence Ka X-rays of Gallium (energy 9-2 KeV.) arising from the 
K-capture of Ge7!. The amplifier gain was reduced by a factor of 4 to 
achieve this calibration. The pulse distribution from this measurement 
on the spectrum is seen in fig. 3, and the intensity below about 
1 KeV. appears to be much less than expected on theoretical grounds. 
We have already mentioned the uncertainty of the curve below ~300 eV. 

The experiment was repeated under almost identical conditions, 
using the pulse analyser to record the spectrum shape. Two measure- 
ments were made, one with a counting rate ~6000 c./min. and the other 
with a counting rate ~20,000 c./min. (fig. 3). It was found that the various 
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points repeated satisfactorily after several hours. In the case of the 
low counting rate, it was not possible to obtain readings below about 
300 eV. because of poor statistics. The channel width was kept at 
100 eV. throughout the experiments, and the energy scale was calibrated 
in the same manner as before. Since the results obtained from these 
determinations with different counting rates agree, it seems unlikely 
that the apparent lack of low energy pulses is due to any limitation of 
the amplifier. The particular limitation in mind was the possibility of 
affecting pulse amplitudes when the pulses occurred within too sort 
an average interval. 

In all the above measurements, the field-adjusting tubes were held 
at the potential appropriate to their diameters, thus ensuring a counting 
tube free from “‘end-effect ”’. Curran et al. (1949 a) have published 
a result obtained using a counter which was not “ end-free ”’ and their 
experimental curve did not demonstrate any tendency to fall steeply 


896 G. M. Insch and 8S. C. Curran on the 


at low energies. To achieve a result comparable with theirs, the 


field-adjusting tubes were earthed and the spectrum shape examined 
(fig. 4) in a manner similar to the two previous measurements. The 
counting mixture giving a counting rate of 20,000 c./min. was used. The 
result obtained is in fairly good agreement with that of Curran et al. 

A“ beaded ” counter was built for further examination of the spectrum. 
This additional step was deemed advisable in view of the somewhat 
surprising result obtained. This counter consisted of a tube 24 in. in 


diameter and 24 in. in length (fig. 1b). The central wire was divided — 


into two sections, one 6in. long and the other 14 in. long, insulated 


electrically from each other by a glass bead, 1 em. long and 1 mm. in © 


diameter. This provided two independent counting volumes having 
similar ‘“ end-effects ’, and the difference between the spectrum in the 
long and short ends represents the distribution which would be obtained 


Fig. 4. 
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from a counter of length 8in. free from ‘‘ end-effect ”. Tritium was 
introduced into the counter, together with a mixture of methane at a 
pressure of 16 cm. Hg. and argon at 2 atmospheres. The counting rates 
were about 20,000 ¢./min. in the long end and 8000 ¢./min. in the short 
end, compared with background counting rates of 600 c¢./min. and 
200 c¢./min. respectively. The spectrum shape obtained from this 
examination is also seen in fig. 3. The pulse analyser was again used, 
and the calibration was carried out in the same manner as previously. 


$4. Discusston. 
The very good agreement between the results achieved by the differing 


methods employed is striking. We believe that the curves shown are . 


accurate to within 10 per cent down to an energy value of 400 eV. and 
that even below this energy the shape of the curve is significant. The 
results obtained are in disagreement with Fermi theory in the low energy 
region, since the theory predicts that the curve should meet the ordinate 
at an intensity value about 86 per cent that of the maximum. 
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There are two possible explanations of the experimental results which 
remove disagreement with Fermi theory. (1) If the counter were not 
100 per cent efficient for detection of very low energy electrons or if the 
energy lost by the particles per ion pair produced varied with the energy 
of the ionizing particles, the true shape of the spectrum could possibly 
be distorted into a shape similar to that obtained. (2) If there were 
present, in coincidence with the fB-particles, quanta or particles of low 
energy, the results could again be explained. 

Fortunately, sufficient data is at hand to discuss the first possibility. 
The work of Kirkwood et al. (1948) together with that of Curran et al. 
(1949 b) seems to establish that in argon, the energy expended on the 
production of an ion pair is constant down to very low energies. 
Additional work carried out in these laboratories with the counter shown 
in fig. 1 (a) has confirmed this result. The K and L X-rays of chlorine 
(energy 2-6 KeV. and 200 eV., respectively) resulting from K- and 
L-capture in A®’ are largely converted by Auger effect. It is found, 
by using a source of A®’ that the energy expended per ion pair by slow 
electrons in argon is constant down to 200 eV. It has also been found 
that the relative intensities of the two peaks are in agreement with 
the results of Pontecorvo et al. (1949). The low energy peak consists 
of two unresolved peaks, one due to the K-quantum escaping from the 
counter, involving the L,, or Ly, shell (200 eV.) and due to nuclear 
L capture involving L, electrons (280 eV.). From these results, it seems 
most unlikely that the counter used is failing to act proportionally down 
to the lowest energy values. 

In connection with the second possibility, a homogeneous quantum 
radiation of energy ~400 eV. in coincidence with the f-particles would 
explain the spectrum shape. Such a radiation would give a peak (Curran, 
Cockcroft and Angus 1949) of shape approximately represented by 
«9 exp (—a), and this radiation, together with the B-rays, would give 
a distribution similar to that observed. The origin of such radiation, 
however, is difficult to explain theoretically. We know that quanta 
of total energy about 60 eV. will be released in the de-excitation of the 
He? ion following B-emission from H’. These quanta could give rise to 
about 3 ion-pairs per disintegration, but we require about 13 ion-pairs 
on the supposition that radiation is responsible for the shape of the 
spectrum observed. It is extremely improbable that the radiation has 
its origin in nuclear excitation of He* in the normal sense of the term. 
It is also possible that a number of low energy quanta (energy ~50 eV.) 
released in coincidence with each other and the B-particles could account 


for the results. 


§5, COINCIDENCE MEASUREMENTS. 


The presence of quanta of energy <400 eV. is extremely difficult to 
establish in a coincidence measurement due to the extreme absorbability 
of the radiations. However, the following method was adopted in an 
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endeavour to ascertain whether or not such radiations were present. 
A counter of 54 in. diameter and 16 in. in effective length (fig. 1 (¢) ) 
was divided into two similar but independently operating volumes by 
a very thin glass bead (~} mm. in diameter and 1:5 em. long). The 
counting tube was lined with aluminium sheet. Tritium was introduced 
into the counter with a trace of He to act as carrier and methane was 
added to give a total pressure of 1 cm. Hg. The counting rate was 
~10,000 ¢./min. in each volume, and the counter operated satisfactorily 
at this low pressure with about 1-5 KV. on the case, giving a high gas 
gain sufficient to detect pulses due to radiation of energy as low as 
40 eV. The value for the mass absorption coefficient of methane for 
low energy quantum radiations does not appear to be accurately known, 
but from the data available it seemed likely that if quanta of energy 
400 eV. were present, coincidences would be obtained. On the other hand, 


Fig. 5. 
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if the quanta had an energy value much less than 400 eV. no coincidences 
would be recorded. The counter was placed with its axis perpendicular 
to a magnetic field of sufficient strength to ensure that there were no 
B-B coincidences due to an electron travelling from one counting volume 
to the other. The results obtained from this investigation were indecisive. 
As the counting tube volts were increased the two counting volumes lost 
their independence and coincidences due to ‘ mutual firing’? were 
recorded. This phenomenon increased in frequency as the gas gain was 
increased. A curve giving the ratio R of the coincidence counting rate 
to the total counting rate in terms of the case voltage was plotted (fig. 5). 
This curve rose steeply at high gas gain. Ka fluorescence X-rays of 
copper were fired through the central window, resulting in an 
approximate doubling of the counting rate in each volume. An exactly 
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similar coincidence curve was again plotted (also fig. 5), and it was found 
that the points on this curve overlapped the points of the previous curve 
to within the accuracy of the experiment. The X-rays of copper gave 
pulses of average size roughly the same as those produced by the tritium 
B-rays. The close similarity of the two sets of observations shows that 
very few, if any, of the coincidences, observed with tritium alone, can be 
ascribed to the hypothetical B-quantum decay process. It is unfortunate 
that this coincidence technique fails to give definite results in the region 
of 400 eV. Clearly, it makes the technique quite inadequate for testing 
the hypothesis of much softer quanta (<400 eV.). 

We are, thus, forced to conclude that the spectrum of H®, as observed 
in a proportional counter, is not in accordance with Fermi theory at 
‘energies below about 1 KeV. We are unable to establish experimentally 
the reason for the discrepancy, although it seems possible that it might 
have its origin in the simultaneous emission of a quantum or quanta 
of total energy about 400 eV. This value is considerably higher than the 
expected energy release (~60 eV.) in the de-excitation of the He? ion. 
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SUMMARY. 

A theory is suggested to describe the distribution of stress and velocity 
in a block of ductile material symmetrically indented on opposite sides 
by two smooth flat dies. The theory is two-dimensional and an ideal 
plastic-rigid material is assumed. Slip-line fields are proposed for all 
ratios of width of dies to height of block greater than one. It is found 
that the velocity field transforms into its corresponding slip-line field. 
Other plane plastic problems where this correspondence occurs are cited. 
Since the velocities are known along only one slip-line boundary, the 
usual numerical procedure for calculating the field has to be modified 
and the use of influence coefficients is introduced. Solutions are computed. 
for width/height ratios between one and two, and it is observed that 
simple analytical relations hold between these ratios and the corresponding 
velocity discontinuities in the field. The calculated mean pressures on 
the dies are compared with upper and lower bound curves derived by 
Hill on the basis of extremum principles. The deformation after different 
finite indentations is illustrated by photographs of experiments using 
(a) plasticine models, and (6) aluminium alloy strip. 

Similar fields are proposed for extrusion or drawing through a wedge- 
_ Shaped die with large reductions, but no solutions have been computed. 


$1. IntRopuctTiIon. 


Tue plane stress-strain curve of a strip, before or after rolling, can 
conveniently be obtained by symmetrically indenting the strip on opposite 
sides by two flat dies completely spanning its width. The test was 
suggested by Orowan and is described by Ford (1948). Provided that 
the indented area is sufficiently long and narrow, the constraint of the 
non-plastic material on both sides of the dies inhibits lateral spread, 
and the deformation of the strip is essentially plane-strain, with the 
exception of narrow regions near the edges (Hill 1950 a). In this paper 
the theory of plane plastic strain due to Hencky and Geiringer is applied 
to analysing the stress and velocity distribution during the course of the 
indentation when the dies are smooth and their width, 2w, is greater 
than the thickness, 2h, of the material (fig. 2). 


* Communicated by the Author. 
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When the dies are smooth and w/h <1, a complete solution has already 
been found and the pressure on the dies for different values of w/h has 
been calculated (Hill 1950 b, p. 256). For perfectly rough dies and 
w/h>3-64, a complete solution was first obtained by Hill, Lee and 
Tupper (1945 or 1951; for a shorter account see Hill and Lee 1946). 
When 3-64>w/h> 1 a solution has been suggested by Hill (1950 b, p. 230) 
which is only valid when there is friction on the dies ; but it is uncertain 
for what range of frictional conditions the solution is appropriate. 

In all these theoretical solutions the material under test is supposed 
to be rigid up to the yield-point ; after which the maximum shear stress 
retains a constant value during plane plastic strain. Such a plastic-rigid 
body cannot deform until the plastic region attains a certain critical 
size (depending on the problem); the load-indentation curve has a 
discontinuity in slope at the yield-point and the load when distortion 
begins is called the yield-point load. When a real material is indented 
by rigid dies the bend in the load-indentation curve is well defined if 
the work-hardening is small. The bend marks the initiation ot large 
plastic strains ; at lower loads the amount of indentation is restricted 
by the elastic resistance of the non-plastic material and all plastic strains 
are of an elastic order of magnitude. The mean load in the interval 
corresponding to the bend should approximate closely to the yield-point 
load of the plastic-rigid material. If the material is annealed, it is not 
possible to define a load which may be compared with that calculated 
by the theory. The solution presented here can therefore be applied 
to the initial yielding of a material with a fairly well defined yield-point 
and to the continued indentation of a non-hardening material. 


§2. Bastc EQuaTION OF PLANE-STRAIN Puastic FLow. 


The theory of plane plastic strain and the slip-line field is set out most 
clearly by Hill (1950 b) who has considerably simplified and improved 
the original presentation by Hencky and Geiringer. Only a brief résumé 
of the basic theory will be given here. 

Since the volume of an element of plastic-rigid material does not alter, 
each incremental distortion in a state of plane strain consists of a pure 
shear. Hence, for the ideal isotropic material, the state of stress at each 
point is a pure shear stress /, assumed constant, together with a hydro- 
static pressure p. If Tresca’s yield criterion is used, k is equal to Y/2, 
while for von Mises’ criterion it is equal to Y/1/3, where Y is the yield 
stress in uniaxial tension. It is assumed that inertial stresses are negligible 
compared with k and that the problem can therefore be treated as quasi- 
static. The plastic stress equations in the plane of flow, consisting of 
the yield criterion and the two equations of equilibrium, prove to be 
hyperbolic. The velocity equations in the plastic region express the 
condition for zero volume change and the fact that for an isotropic 
material the principal axes of stress and plastic strain-increment must 
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coincide. These equations are also hyperbolic and their characteristics 
coincide with those of the stress equations. There are two orthogonal 
characteristic directions at any point in the plastic field; they are the 
directions of maximum shear stress or shear strain-rate. The two 
orthogonal families of curves whose directions at every point coincide 
with those of the maximum shear stress are known as slip-lines and will 
be labelled «- and f-lines. By convention, if the «- and f-lines are regarded 
as a pair of right-handed curvilinear axes, the line of action of the 
algebraically greatest principal stress falls in the first and third quadrants. 
The equations governing the variation of p along a slip-line, which were 
first derived by Hencky (1923), are 


p+2k¢=constant, along an «-line, a5 
1 


p—2kd=constant, along a f-line, 


where ¢ is the anticlockwise angular orientation of a slip-line with respect 
to some fixed direction. These relations are completely equivalent to the 
equilibrium equations. 

The equations of Geiringer (1930) for the components of 
velocity u, v, along the « and f slip-lines, are 


du—vd¢d=0 along an «-line, 
dv-+-udd=0 along a f-line. j o 
They state that the rate of extension along any slip-line is zero. 

It can be shown that the normal components of velocity may be 
arbitrarily prescribed without inconsistency along two intersecting 
slip-lines. It follows that it is possible to have a discontinuity in the 
tangential component of velocity across a slip-line. This implies a 
momentarily infinite rate of shear; but this is permissible since work- 
hardening is neglected and a slip-line is a direction of maximum shear 
strain-rate. The discontinuities correspond to what in a real material 
would be a narrow transition region, where the shear strain-rate is large. 
It is evident that such transition regions will tend to be less sharp in a 
metal with a high rate of work-hardening, which emphasizes the fact that 
this theory is not applicable to such a material. 

Hencky’s first theorem, which is a direct consequence of equations (1), is 
the geometrical expression of the fact that the pressure difference between 
two points, calculated by proceeding along either of the two slip-line 
paths between the points, is the same. If we consider a curvilinear 
quadrilateral APBQ (fig. 1) bounded by two «-lines, QA and BP, and 
two f-lines, QB and AP, we deduce from (1) that 


Pa Pe bq = Opi ar eae ae ee ea) 


Any two families of orthogonal curves satisfying this condition 
constitute a possible slip-line field for a plastically deforming mass under 
certain boundary conditions. In general the plastic region is divided 


Compression of a Ductile Material between Smooth Flat Dies 908 


into two parts by one or more slip-lines: a part where deformation is 
occurring, and a part which is rigid. The plastically deforming part 
consists of a number of subsidiary domains bounded by slip-lines across 
which the curvature of the slip-lines and/or the tangential components 
of velocity change discontinuously. The plastic material in the rigid 
part is prevented from deforming by the constraint of surrounding 
non-plastic material. In order to solve any problem the slip-line field 
must be found in the plastically deforming zone. The probable approximate 
location of this zone is first suggested by physical intuition, aided when 
possible by experiment. The slip-line field is then constructed to satisfy 
all the stress and velocity conditions that directly concern the zone, as 
well as the condition that the rate of plastic work should be everywhere 
positive. Finally, the solution is valid if the associated stress distribution 
in the rigid material is such that the yield limit is nowhere exceeded. 


Fig. 1. 


Slip-line field network showing (x, y) coordinates, 


The most convenient coordinates for defining a point in a slip-line 
field are the quantities (x,y). They are the coordinates of the point 
P: under consideration, referred to axes passing through some fixed 
origin 0, and parallel to the respective slip directions at P (fig. 1). It 
is easy to show that for any two orthogonal families of curves 

dy+azd¢=0 along an «-line, 
dx—yd¢=0 along a B-line. | 

In applying these conditions to a slip-line field the geometrical property 
(3) is introduced which determines ¢ at all points of a network of slip-lines. 

The stress and velocity boundary conditions are particularly easy to 
deal with when the material is in contact with either a smooth or a 
perfectly rough rigid surface. The slip-lines meet a smooth surface 


(4) 
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at 45°, whereas a perfectly rough surface is cut orthogonally by one 
family of slip-lines since the frictional stress is equal to k. These 
conditions are independent of the normal pressure and we can therefore 
proceed with the solution of the slip-line field without specifying 
beforehand the value of p at any point. The only other simple frictional 
condition is that the shear stress on the surface has some constant value 
less than k. Each slip-line then intersects the surface at the same angle. 
This condition implies, however, that the coefficient of friction decreases 
with increasing normal pressure on the surface, which is an artificial 
condition not likely to occur in practice. If the coefficient of friction is 
constant, the angle of intersection varies with normal pressure along the 
surface of contact. This is a difficulty inherent in the solution of problems 
involving intermediate frictional conditions. 


§3. Frypine THE Suip-LinE FIELDs. 


In this indentation problem geometric similarly is not maintained, 
and therefore we have to follow, from the very beginning, the progressive 
changes in the configuration. At each stage of the indentation there 
will be a different slip-line field, with its corresponding stress and velocity 
solution. The initial shape of the overhang does not affect the plastic 
problem, but for convenience we shall discuss the symmetrical indentation 
of a rectangular block. 

The development of the plastic zones and the momentary boundary 
between the rigid and the plastically deforming material has been discussed 
in general terms for indentation by rough dies (Hill and Lee 1946). The 
present problem is similar in many respects. It may be expected that 
plastic zones will be initiated at the sharp edges of the dies and will 
spread inwards until they meet in the geometric centre of the block. 
Since, however, in this problem there are no frictional stresses, it does 
not seem possible that there should be a non-plastic area left in contact 
with the dies at the centre. For the solution presented here, no 
deformation consistent with Geiringer’s equations is possible until the 
plastic zone has spread completely through the block. Plastic deformation 
is also restricted by the overhang, part of which must be non-plastic 
since its surfaces are stress free. It will be assumed that the overhang, 
even if partly plastic, moves out as a rigid whole. This was the 
assumption made by Hill and Lee, and it is confirmed by experiment 
(Nye 1947; see also the plasticine experiments described in this paper). 
It follows that the external plastic-rigid boundary must pass through 
the edges of the dies. 

For integral width-height ratios the solution, suggested by Hill 
(1950 b, p. 232), is very simple. The block deforms momentarily, as a 
number of independent rigid units, which slide along a criss-cross of 
straight slip-lines, which meet the smooth dies at 45°. If the dies move 
with unit speed, there is a constant tangential velocity discontinuity 
of 1/2 initiated at the block centre and propagated by successive reflection 
from the dies along the whole length of the block. The solution is only 
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valid for integral values of w/h because for all other dimensions the 
discontinuity terminates on the exist slip-lines, which is incompatible 
with the rigid body motion of the overhang. The solution for intermediate 
values of w/h must be such that it passes continuously into the simple 
straight line criss-cross as w/h approaches an integral value. 

Consider first the general form of the field in the upper left-hand 
quadrant of the block, as determined by the boundary conditions (fig. 2). 
Velocity conditions at the plastic-rigid boundary, OB, determine that 
this boundary must be a slip-line, and the sense of the relative straining 
on either side of OB. requires that it should be an «-line. Once the 
position of OB is chosen, the slip-line configuration is uniquely determined 
by equation (3) and the stress boundary conditions. These are that the 
slip-lines cut the axes of symmetry, DE and EB, at 45°, and that along 
OD they meet the smooth surface at 45°. 

Assume for the moment that there is a stress singularity at O, and that 
the radius of curvature of OB does not change sign along its length. 
OB must then be convex towards the overhang, and the «-lines within 
a certain angular span pass through O. The field defined by OB and the 


Fig. 2. 


Rectangular block during indentation. 


singularity at O can be continued round O until the «-line that meets the 
surface of the die at 45° is reached. To avoid ambiguity, the point 
singularity will in future be denoted by O, or A, according to whether 
the point is considered to lie on the plastic-rigid boundary or on the 
surface of the die. Let the difference in ¢ be equal to % between O and A. 
The geometry of the boundary conditions on AD and BE determines 
that the angle turned through by OB is also ¢. These boundary conditions 
also make it possible to build up successive domains of the field, the 
boundaries of which meet the surface of the die at D,, D,, etc., and the 
horizontal axis of symmetry at E,, E,, etc. (fig. 3). The angle turned 
through by each domain boundary is equal to 4. The vertical axis of 
symmetry must, clearly, be one of the domain diagonals D,E,, D,Ks, ete. 
The number of domains through which the field must be extended depends 
on the value of w/h. For reasons of symmetry, velocity discontinuities 
can only take place along either AE,D,E;.. . or OBD,E,D; soe 
This distinguishes two types of field, I. and II., with their velocity 
discontinuities, V, and V,, ending at A and O respectively. 


. SER. 7, VOL. 42, NO. 331.—AUGUST 1951 3Q 


906 A. P. Green on a Theoretical Investigation of the 


Fortunately, it is not necessary to obtain a complete solution of the 
velocity field in order to verify that the form of the slip-line field developed 
above is compatible with the velocity boundary conditions. The 
verification depends on the fact that, in this instance, the velocity field 
can be transformed into its corresponding slip-line field by means of the 
relations 

Uae, V=Y 5  a-sf. aoa teeth ee oe 


It can be seen, by comparing equations (2) and (4), that the variation 
of uw along an «-line corresponds to the variation of x along a f-line, 
provided that the direction of change of ¢ is the same ; v and y correspond 
in a similar manner. The normal boundary velocities w along a f-line 
transform into the x values along an «-line boundary. Now consider 
a tangential velocity discontinuity 4u across an «-line, which means 
that along any f-line as it crosses the discontinuity » and ¢ remain 
constant, and u increases by 4u. Upon transformation, this becomes 
a straight «-line region in the (x, y) field of length du. Similarly as we 
travel along any straight «-line through the region, y and ¢ remain 
constant and x increases by du. The velocities along a plastic-rigid 
boundary transform into a singularity in the slip-line field provided that 
the rigid part is not rotating; in addition, if there is a velocity 
discontinuity along the boundary, it becomes a circular are with the 
singularity as its centre of curvature. 

To make the correspondence exact in the present problem, the scales 
and origins of the two fields must first be adjusted. Let the outward 
speed of the rigid overhangs be U, corresponding to unit speed of the dies. 
The upper left hand quadrant is given a velocity U parallel to the dies 
towards the centre of the specimen so that the speed of the overhang 
is reduced to zero. The scale of the slip-line field is adjusted so that the 
half thickness, h, is always unity. The width of the dies will then be 
equal to 2U. We choose O as the (x, y) origin. The boundary conditions 
for the velocity and slip-line fields can then be written as follows :— 


Boundary condition Slip-line field | Velocity field 
OA, f-line, (x, y) origin. x=y=0 a 
OB, w-line, velocity of rigid overhang is zero. — 400) 
AD, slip-lines meet die at 45°; vertical xu=y U—v=/2 
velocity is unity. 
BE, axis of symmetry ; vertical velocity w—y=/2 U=v 
is zero. 


DE, axis of symmetry (AD=U) ; _hori- 
zontal velocity is U, “+y=U/2 ut+v=U/2 
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It is clear that if the relations (5 a) are applied, the velocity boundary 
conditions transform into those of the (x, ¥) field, except that certain 
boundaries are interchanged. The zero boundary velocities v and w on 
OB, which is an «-line, correspond to the (2, 7) origin OA, a fB-line. The 
PomtesOA, BLD Eines DE... in the velocity field correspond to 
the points O, B, A, E,, D,, . . . E,, D,, respectively, in the (2, y) field. 
The directions of change of ¢ are the same at all corresponding points. 
It is reasonable to expect, therefore, that the solutions of both fields 
are identical. This implies that there are straight line regions in the 
(x, y) field corresponding to the velocity discontinuity. The two types 
of slip-line obtained by transforming the two types of velocity field, 
I. and IL, are shown in fig. 3. Although it has not been proved that 
the rate of work is positive everywhere in the field, it is reasonable to 


Fig. 3. 


B E, E, Ea 
Two types of slip-line field (not to scale). 


assume that this is so since the directions of shearing across all the 
velocity discontinuities imply positive work. An examination of the 
already existing solution for w/h<1 shows that here also the two fields 
are equivalent. Discussion of the transformation of fields is continued 
later in the paper; for the moment, we confine our attention to the 


indentation problem. 


§4, NUMERICAL CALCULATION OF THE FTELDS. 


The basis of the usual numerical method of solution is the replacing 
of the differential form of equations (2) by linear difference relations 
across a network of slip-lines. In the present problem we are able to 


302 
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choose the network so that the increments of slope 4¢, between con- 
secutive points, is the same throughout the field. This is possible because 
the angles turned through by the boundaries of all domains are equal 
to %. The field is thus split into an assembly of equal unit cells of which 
a typical one is shown in fig. 1. In order to calculate the velocities at 
P in terms of those at A and B, equations (3) are replaced by 
Up—Up=zAAg(vp+ Up), 
Up—Va= —pwAd(Upt+Ua); 
where 4d denotes the numerical value of the change in angle, and A and 
u must be assigned the values of 1 or —1, according to whether ¢ increases 
or decreases in travelling towards the point P along the respective 


slip-lines from A and B. Solving for up and vp, and retaining powers 
of 4¢ only up to the third, there results 


up=(1— Pw Ag) y+ PAAG(1— FAAP") (04 +p) Pe gPug, 


vp=(1— Jw Ag), — WAG (1 Pw A$2)(u44 + 4_) Ps Adv. | 
If the normal components of velocity are known at net points along two 
intersecting slip-lines, the velocities at all net points within the domain 
defined by these slip-lines can be calculated by repeated application of 
equations (7). 

We shall now consider how these equations can be used to obtain 
numerical solutions of the velocity field for width/height ratios between 
one and two. The boundary velocities are known initially only on OB, 
the outer slip-line on the rigid side of which uw and v are zero. Neither 
y, nor the velocity discontinuities that correspond to different values 
of U, are known in advance ; and before solutions over the whole range 
of U have been completed, even the range of % can be guessed only 
approximately. There is, clearly, not sufficient information to calculate 
the velocities directly through the field, in the same way that the velocity 
solution was calculated, for example, in the rough die indentation 
problem. It is necessary first to try normal components of velocities 
on a slip-line adjacent to OB. -OA was chosen because the velocities 
along this slip-line correspond to the (x, y) coordinates along OB, which 
are needed when calculating the stress distribution. OA and OB were 
divided into five equal arcs so that 4¢ throughout the field was equal 
to p/5. 

The normal component of velocity on OA at O, within the plastic 
region, is either zero, or Vg, according to whether a solution of 
type I. or I. is being considered. At the five successive net points on 
OA after O, the normal components of velocity, u, were taken as 
a, b,c, d and e. Using the second of equations (6) .and the condition 
that v=0 at O, a step integration along OA gives the values of v at the 
net points on this boundary. Since the equations (7) are linear in w and 
v, the calculated velocities at all net points are linear functions of the 
norma) boundary velocities on the starting slip-lines OA and OB, 


(6) 


(7) 
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Therefore, groups of boundary velocities containing the different unknowns 
can be integrated through the field separately and the results superposed 
at each net point. We require a number of solutions in the selected 
range in order to be able to draw a graph of the mean pressuse on the 
dies against 1/U. Therefore, instead of repeating the computation 
several times for chosen values of A¢, it is more convenient to compute 
in terms of 44, using equations (7), and neglecting powers of 4d higher 
than the third. This can be done fairly rapidly provided an efficient 
method of tabulating the results is used. 

The velocities on AC, and BC, are first calculated by integrating through 
the domain OAE,B. The boundary conditions on the diagonals AD, and 
BE, enable the calculation to be continued through AC,D, and BC,E,. It 
should be noticed that the boundary condition on AD, introduces a 
numerical term which must be integrated separately from the terms 
involving the various unknown velocities. The field in the domain 
C,D,E,C, is defined by the calculated velocities on C,D, and C,E, and 


Fig. 4. 
8B a 
-£459 gene 
4 -22+ 4% 
] ri a 
(e) (b) «) 


is computed as far as the six diagonal points on D,E,. The velocities 
on this diagonal must satisfy the condition, w-+v=U,/2, and the 
velocities at A, in contact with the die, must satisfy the condition, 
u—v=1/2. We thus have seven equations for seven unknowns, 
a, b, c, d, e, U and either V, or V,, which are. solved numerically, 
substituting a definite value for 4¢. 

The process of integrating through the different domains of the field 
is considerably shortened by the use of “ influence coefficients ”’. ; These 
coefficients are the calculated velocities at net points of a domain, due 
to a unit normal velocity at a single net point on the initial boundary. 
The tangential velocities on the initial boundary are found in the usual 
manner by applying one of the equations (6). Three typical groups of 
initial boundary velocities are shown in fig. 4. The integrations of 
grounds (b) and (c) of this figure are clearly similar. The field of (c) is 
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the same as that of (6) shifted one step along the network. All the 
influence coefficients for the network can therefore be obtained by 
integrating groups (a) and (b). As we have already observed, the 
linearity of the velocity equations enables any group of velocity boundary 
values to be “ integrated ” through the domain by simple super-position 
of their individual effects. We can, if we wish, proceed directly from the 
initial to any final boundary chosen in the domain, without having to 
calculate the velocities at any intermediate net points ; all we need to 
know are the influence coefficient’ on the final boundary, corresponding 
to each net point on the initial boundary. 

For the solution of the present problem influence coefficients were 
worked out in terms of 4¢ (retaining powers of 4¢ only up to the third), 
for integration between the following boundaries :—OA to AC,, and 
BC,; BC, to BE,, and C,E,; AC, to €,D,; C,D, to D,E,; C,H, to 
D,E,. The reason for calculating the velocities on the diagonal BE,, 


Fig. 5. 


PRESSURE 
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r\|= 


SURFACE OF BLOCK CENTRE 
Pressure distributions. 


where u=v, is that this corresponds to the surface of the die AD, in 


the slip-line field, along which x=y. The coordinates are required ’ 


along AD, in order to calculate the mean pressure on that surface. 
The use of influence coefficients reduced the labour of computing the field 
by about a half. For a problem in which one fixed value of A¢ is 


sufficient, it would probably be more convenient to calculate the influence’ 


coefficients numerically, and the proportion of time saved would probably 
be greater than in this problem. 

When the slip-line configuration has been completed the values of the 
mean compressive stress p at all net points can be found from Hencky’s 
equations (1), provided that the value p, at one point of the field, say, 
Bis known. This value is provided by the condition that the horizontal 
component of the total force acting over OB is zero. Fig. 5 shows the 
distributions of pressure on the die at different stages of the indentation. 
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The relation between the mean pressure on the die and the ratio 
hyw=1/U, is shown by the full curve in fig. 6. 


§ 5. Discussion or REsuuts. 


The velocity fields, and hence the slip-line fields, of types I. and IT., have 
been determined for values of equal to 5°, 10°, 15° and 19°30’. In 
fig. 7, the slip-line fields of both types, for % equal to 15°, are shown to 
‘scale. In order to plot these fields in detail the shapes of the initial 
slip-lines, OB, were found by the methods described above. The fields 
were then re-calculated numerically, using an equiangular network 
with 4¢ equal to 5°. This was quicker than working out the coordinates 
of net points by substituting known values of 4¢ and the initial boundary 
values in the original calculations. 


Fig. 6. 
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Calculated mean pressures and the two bounds. 


The following simple relations are found empirically from an 
examination of the numerical results :— 


U24+4V2—U2_Vi=e ; 
UUs=2 


ee are oes amen se |. (8) 
and a : 


The suffixes 1 and 2 refer to the two types of solution, I. and IL., 
respectively. The accuracy with which these relations are satisfied by 
the calculated values is shown in Table I. It can be seen that the 
discrepancy is nowhere greater than 0-02 per cent. This is a useful 
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check on the accuracy of the calculations. It seems likely that fairly 
accurate results would have been obtained even if the cube powers of 
Ad had been neglected throughout the solution. The change over from 
the solution of type I. to that of type I. takes place smoothly when 


Fig. 7. 
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Slip-line fields to scale. 


U;=U,=~/2 and V- i Var: Since the last two equations of (8) express 
“ ae pes Aebee in terms of U, and P,, it appears that only one type 
of field need be solved to obtain the variation of P 

tea aes over the whole range 
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PLASTICINE EXPERIMENTS. 
Fig. 8. 


(b) w/h=12. 


wlh=/2->2. 


wlhh= \/6>1/12. 


Fig. 11. Fig. 12. 


wlh=1—1°34. whh=/2—>1-86. 


To fuce page 912 
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A theoretical explanation of equations (8) is not yet known. It seems 
probable that the general forms of these relations, valid for all U greater 
than unity, are 

sit 
Ui+ 3 Vi=m(m-+1), where m<U,<v[m(m-+1)], 


m-+1 
ui-( )¥i=m(m-+1), where /[m(m+1)] <U,<m+1; 


2 
U,U,=m(m-+1), - for the same m, if 4,=d, 
P,=P,, m=1,2,... (9) 


Tt should also be noted that when U is less than unity but greater 
than 1/8-74 
U?—iV?=0, 


which is in agreement with the second of equations (9) when m=0. 


TABLE I, 


¢ | U, | V, | U, | V, /U2+3V2]U3—Vv3] U,U, | Py/2e | P,/2e 


0 1 2 2 2 2 2 2 1 1 


5° | 1-1646 | 1-1348 | 1-7174 | 0-9744 | 2-0002 | 2-0000 | 2-0001 | 1-0106 | 1-0108 


10° | 1-2985 | 0-7922 | 1-5400 | 0-6098 | 1-9999 | 1-9993 | 1-9997 | 1.0277 | 1-0279 


15° | 1-3859 | 0-3977 | 1-4429 | 0-2858 | 1-9998 | 2:0003 | 1-9997 | 1-0376 | 1-0377 


19° 30’| 1-4140 | 0-0118 — — Pegg aaa — 1-0389 | — 


§ 6. Upper AnD Lowzr Bounps ror THE MEAN PRESSURE, 


It is interesting to compare the calculated mean pressure on the block 
with upper and lower bounds derived by Hill}; they were obtained by 
the application of two extremum principles for a plastic-rigid body 
(Hill 1950 c), extended to the case where part of the body is non-plastic 
by Drucker, Greenberg and Prager (1950). 


The lower bound is 
pe cae wera tees. I~ (10) 


which expresses the otherwise obvious fact that it is not possible for the 
mean pressure to be less than it is in uniform compression. The upper 
bound is obtained by applying the inequality 
: [FyutdS<k{fy*dV-+ f[u*]ds}. 2. we (LY) 
ee ee Se a ee os 
+ Private communication, 
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F, is the actual distribution of external stresses on the surface 8. The 
asterisks denote a fictitious velocity distribution, u;, which satisfies the 
boundary conditions and the conditions of no volume change ; y* denotes 
the corresponding maximum shear strain-rate in a volume element dV ; 
[w*] is the absolute value of the velocity discontinuity, and the 
corresponding integral is taken along it. The fictitious velocity solution 
used to obtain the upper bound of P is as follows :—the block is divided 
symmetrically into a number of independent rigid units, which slide 
along a criss-cross of lines inclined at the same angle 6 to the faces of the 
dies. For integral values of U this is, of course, the actual manner of 
deformation with @=45°. The integral involving y* is now zero, and 
[u*] is equal to 1/sin @ between the rigid units. Assuming that there are 
m units in contact with each die, it follows from (10) that 

r ds km 

) sind ~ sin’ 
Hence the least upper bound is defined by 


P line U 
<3 (Uta) o\ \ia al ey Deen E<ce (12) 


| Futas—PU <k where cot 6= = 


where [m(m—1) ]}<U<[m(m+1) ]; m=1, 2, .. . 

The two bounds (10) and (12) are shown as dashed lines in fig. 7. 
The closer the similarity between the fictitious and the real velocity 
fields, the nearer the upper bound approaches to the actual mean 
pressure. Peak values of both the upper bound and the calculated 
mean pressure occur when U= 1/2, and here the discrepancy between 
the two curves is greatest. This might be expected since the actual 
velocity field for this value of U involves no velocity discontinuity and 
is therefore very different from the fictitious field. The peaks of the 
upper bound are closely approximated by the expression 

3s 1 

ak T BUF 
which rapidly approaches unity as U increases. There is no need, therefore, 
to extend the accurate solution of this problem beyond U=2, since for 
all U greater than 2, P/2k differs from unity by less than 2 per cent. 


§7. INDENTATION EXPERIMENTS. 
(a) Plasticine. ; 
The instantaneous velocities can be found at net points of the slip-line 
field for any chosen value of w/h, but it would clearly be a very laborious 
process to calculate the distortion of a square grid after a finite indentation. 
It is comparatively easy, however, to obtain an approximate picture of 
the deformation by the use of plasticine models : for it has been shown, 
by comparing plasticine model experiments with calculations made on 
the basis of the theory for a non-hardening plastic-rigid material, that 


in plane strain plasticine deforms in a similar manner to such an ideal 
material (Green 1951). 
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Rectangular plasticine blocks were made in two halves, tin. thick. 
On the inner face of one-half a 1/10 in. square grid was lightly stamped 
in ink. The two halves were clamped between parallel glass plates 
lubricated with vaseline to ensure plane strain conditions, and French 
chalk was used to prevent the two halves sticking together. The blocks 
were indented from one side only, the horizontal plane of symmetry 
BE being replaced by a plane lubricated base. The photographs of the 
deformed grids after indentation between different initial and final 
values of w/h are shown in figs. 8, 9 and 10 (PI. XXIV.). For 1 <w/h <1/2 
the deformation is mainly confined to regions of intense shear radiating 
from the edges of the die (fig. 8a). As indentation continues, the 
mode of deformation clearly changes from type I. to type II., as shown 
by three relatively undeformed triangular regions in fig. 8b, where w/h 
has reached the value 2. This is also illustrated in fig. 9 where the 
formation has been of type II. throughout. Fig. 10 shows a block 
which has been compressed from w/h~4/6 to w/h~ /12; there are five 
triangular regions in which the deformation has been small, which is 
what we should expect from the theory. 

(b) Metal}. 

The deformation of metal strips symmetrically indented on opposite 
sides by two flat dies is illustrated in figs. 11 and 12 (Pl. XXIV.). Aluminium 
7 per cent magnesium alloy was used because the deformation can be 
clearly demonstrated with this metal. When it is rolled the segregated 
impurities in the metal are strung out in straight lines along the length 
of the strip. The distortion of such lines by indentation was shown up 
by polishing the central longitudinal section of the strip, and etching it 
with acid ferric chloride. The width of the strip was five times the 
width of the dies, which were made of tungsten carbide and lubricated 
with calcium oleate. Since the metal had been cold worked by rolling 
it to the required thickness, the rate of work-hardening during 
indentation was small. If the non-dimensional quantity H is defined 
as the mean slope of the stress-strain curve over the range of strain 
during indentation, divided by the mean stress in uniaxial compression 
over the same range of strain, then H was equal to approximately 
0-4 and 0-25 for the indentations shown in figs. 11 and 12 respectively. 
The flow patterns are in qualitative agreement with the theory and are 
similar to those shown in figs. 8a and 9, the corresponding plasticine 
experiments, though in the metal the regions of intense shear are more 
sharply defined. 

It will be noticed that experiment confirms the assumption made in 
the theory, that the overhang moves out as a rigid whole: the material 
outside the initial area of deformation remains undisturbed and the 
general shape of this region agrees with the theory. All the photographs 
demonstrate the inhomogeneous nature of the deformation. It is clear, 


a Aire cents oe ees a Pe 
+ Metal experiments carried out by the Metallurgy Section of the Metal 
Flow Research Laboratory, Sheffield, under the supervision of Mr. L. Bourne. 
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therefore, that severe work-hardening will modify the nature of the 
flow. For example, in similar experiments with annealed metal strip 
the regions of large shear were much more diffuse, though the type of 
deformation was similar. 


§ 8. TRANSFORMATION OF Fre~ps—DIscuUSSION AND APPLICATION. 


In addition to equation (5a) alternative relations for transforming 
the velocity field are : 


u=y, v= —2z', ara’, d=?’ (Ob), 
u=2z, v= —2z', «fp, ¢=—¢' (C); ese et) 
w=y', v=7', ae’, $= —#' (0). 


Starting with the velocity field of any quadrant in the present problem, 
these transformations all result in an (2, y) field of the same shape as that — 
of the original quadrant. However, either the «- and f-lines are inter- 
changed, in effect reversing the directions of the maximum shear stress 
throughout the field so that the work of distortion would be negative 
(6 and c); and/or the field is the mirror image in an axis of symmetry 
of the field in the original quadrant (c and d). 

The relations (5) can be used to transform any velocity or slip-line 
field. In general, the velocity field cannot be expected to transform 
into its own slip-line field ; but it may transform into a slip-line field 
which has already been computed, thus saving a considerable amount 
of labour. For example, the transformed velocity fields for symmetrical 
steady motion problems, such as drawing or extrusion, sometimes involve 
the field defined by equal circular ares ; this is a standard field that has 
already been computed and which can also be found analytically. In 
certain problems with particularly simple stress and velocity boundary 
conditions the two fields will be equivalent. This occurs, for example, 
in the following plane strain problems :—symmetrical indentation of 
a block on opposite sides by two rough flat dies for w/h greater than 
unity ; extrusion or drawing through a smooth square or wedge-shaped 
die (Hill 1948 a, Hill and Tupper 1948); indentation by a smooth wedge 
of the plane surface of a semi-infinite block (Hill, Lee and Tupper 1947) ; 
expansion of a semi-cylindrical cavity in the surface of a semi-infinite 
block (Hill 1948 b). The last two are rather trivial examples of this 
equivalence. 

An interesting example is the problem of sheet-drawing through a 
smooth wedge-shaped die. The solution which is valid for lower 
reductions fails for certain higher reductions, when the velocity 
discontinuity terminates on the exit slip-line (Hill and Tupper 1948). 
The situation is analagous to that encountered in the smooth die 
indentation problem when w/h has non-integral values greater than 
unity ; the slip-line configuration which overcomes the difficulty is 
similar in form (fig. 13). The (x,y) and velocity solutions are again 
identical if the origin of the (x, ¥) field, O, is taken on the axis of symmetry 
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at the point of intersection of the tangent to the die faces. No numerical 
calculations of the field or of the corresponding upper bound have yet 
been made. 

The importance of considering the velocity field associated with a 
slip-line field is emphasized by the solutions presented in this paper. 
The transformation of fields, and the existence of equivalent fields do 
not appear to have been discussed in the same manner previously. 
Geiringer (1937) speaks of “ conjugate ” solutions, but these involve the 


Fig. 13. 


' Slip-line fields for wedge-shaped die (not to scale). 


correspondence between one complete solution (velocity and slip-line 
field) and another. The transformation equations of the two kinds of 
conjugate solution which she discusses can be written 


x 5 Nhe mee Ist kind, 
y= —v ; w=, 
ear 2nd kind, 
S= uw: v= -—R, 


(she actually uses cartesian coordinates for both position and velocity). 
R and § are the radii of curvature of the «- and f-lines at a point. They 
satisfy the relations 

dS-+Rd¢=0 along an «-line, 

dR—Sdd=0 along a f-line, 
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The particular conjugate solutions which she derives are mainly of 
academic interest. Transformation relations exist, similar to those 
expressed by equations (5), but involving (R, 8) coordinates. These, 
however, are not likely to prove so useful in practice, since the stress 
boundary conditions can usually be expressed more easily in terms of 
(x, y) coordinates. 


§9. CONCLUSIONS. 


The value of the investigation presented in this paper, apart from the 
interest of the new solutions presented, lies in the formulation of new 
methods for solving plane plastic problems. The transformation of 
velocity fields and the use of influence coefficients will often save 
considerable computational labour. Provided that initially the form of 
the field can be tentatively built up, the possibility of integrating unknown 
boundary values from an initial to a final boundary overcomes the 
difficulty which arises when insufficient boundary values are known for 
a direct solution. The element of trial and error is thus reduced though 
not eliminated. The method due to Hill for obtainmg upper and lower 
bounds for mean applied loads is likely to prove most useful in some of 
the many problems of which the exact solution cannot be found or would 
be extremely laborious to calculate. 
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XCIII. The Development of Deformation Teatures in Metals.— 
Part II. Hexagonal Structures. 


By E. A. Cannan, B.Sc., and C. J. B. CiEws, Ph.D., 
National Physical Laboratory*. 


[Received May 11, 1951.] 


§ 1. INTRODUCTION. 


Hexaconat close-packed metals slip on a single crystallographic plane 
and not on a system of equivalent planes as is the case for cubic metals. 
Due in part to this simple slip mechanism, considerable attention was 
devoted to the deformation of the hexagonal metals, notably zinc and 
magnesium, by the German school during the decade 1925-1935 (Schmid 
and Boas 1950). There is, however, an additional factor important in 
the deformation of these metals, namely, the frequent occurrence of 
deformation twins. The crystallography of the slip and twinning 
processes has been well established, and the main features of the 
deformation textures observed have been explained in general terms, 
but there appears to have been no completely systematic treatment on 
the basis of slip rotations and twinning re-orientations. Mathewson and 
his co-workers have attempted to rationalize deformation textures as 
a whole, but their method is not applicable to hexagonal metals (see for 
example, Hibbard and Yen 1948). 

The purpose of the present paper is to show how the effect of twinning 
may be introduced into the treatment already developed for face-centred 
and body-centred cubic metals (Calnan and Clews 1950, 1951, 
subsequently referred to as Parts I. and II.), and how this treatment may 
be applied to the prediction of the deformation textures in hexagonal 


metals. 
§ 2. DEFORMATION MECHANISMS. 


Slip. All hexagonal close-packed metals deformed at room temperature 
slip on basal planes, (0001), in the <1120> directions. Approximate 
resolved shear stress contours are shown in fig. 1. It will be noted that 
the shear stress varies from the maximum possible value, 0°50 of the 
applied stress, at 45° from the pole of the slip plane and the slip direction, 
to zero at 0° and 90° from the pole of the slip plane. The calculation 
of accurate values of the cos y cos A function, where y and A are the 
angles which the slip plane normal and the slip direction make with the 
applied stress direction, was not necessary owing to the simplicity of the 
system. The rotations under stress are of the, usual kind leading under 
tension towards the slip direction and under compression towards the 
pole of the slip plane. Tension rotations are shown in fig. 4. Most of 


eer ere en ee ee es ee ee ee 
* Communicated by the Authors, 
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the rotations are those associated with single crystal slip, but at those 
orientations for which the effective stress, T,,, as defined in Parts I. and IL., 
may reach the boundary [0001] to [1010] duplex slip can occur. There 
is no position corresponding to multiple slip. 

Twinning. The orientation of a twin relative to the matrix may be 
considered as the result of rotation through 180° about a specific axis, 
the twinning plane normal. This does not, of course, imply a bodily 
rotation of the lattice, for the atomic movements correspond to a shear 
of each atom layer parallel to the twinning plane together with, in many 
cases, atomic movements perpendicular to this plane. The macroscopic 
effect of these movements is one of simple shear in the twinning direction. 
The twinning planes in hexagonal close-packed metals are of the form 
{1012}, and for the (1012) plane the twinning direction is the intersection 
of this plane with the (1210). 


Fig. 1. 


OOO! Ke){e) 


Unit triangle of standard hexagonal stereographic projection showing resolved 
shear stress contours and some lines of quickest descent for slip on (0001) 
planes in the [1120] direction. 


In establishing the conditions under which twinning occurs, it is 
necessary to take account of the fact that some metals have axial 
ratios, c/a, less than the ideal close-packed value of 1/3, e. g. magnesium, 
beryllium and titanium, while for others c/a exceeds this value, e. g. 
zinc and cadmium. Figs. 2 (a) and 2 (b) show the geometry of twinning 
for these two cases. For convenience, the positive sense of the twinning 
direction is taken as that lying on the opposite side of the c-axis from the 
twin plane normal. For c<,/3a, fig. 2 (a), it is obvious that twinning 
will occur under tension only when the stress is so oriented as to give 
a component in the positive sense of the twinning direction, OD, The 
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(0) 
Diagrams illustrating the geometrical conditions for twinning, (a) for c< 4/3a, 


and (b) for e>4/3a. The twinning direction, the twin plane normal, and 
the c-axis are all in the plane of the diagram. 
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resolved shear stress for twinning is T cos ycos6 where y and 6 are 
respectively the angles which the stress makes with the twin plane 
normal and the positive sense of the twinning direction. Thus for 
twinning T cos y cos 5 must be positive. A tension of sufficient magnitude 
in the direction T, will produce twinning, but one in the direction T, can 
never do so. Consideration of fig. 2b for c> +/3a shows that the condition 
for twinning under tension is that T cos ycosé shall be negative. 
Conversely, under compressive stresses these two conditions are 
interchanged. 

Boundaries between regions of slip and twinning. The boundaries 
between regions of slip and twinning are the loci of points at which the 
ratio of the resolved shear stress to the critical stress for slip is equal to 
the corresponding ratio for twinning, 7. é. 


T cos x cos A/C,=T cos y cos 8/C,, 


where T is the applied stress and C, and C, are the critical stresses for 
slip and twinning. The equation to the boundary is thus 


cos y cos.d6= COs y COB A, = =.= ae ee 


where n=C,/C,. Inserting now the sign conditions for twinning, 
(1) becomes . Y 
cosy cosd= -+ncosxycosA ..... . (la) 


for c <1/3a under tension, 
and c> 4/3a under compression, 


and 
cosy cosd= —ncosxycosA ... . . (16) 


for ¢<4/3a under compression, 


and c> 4/3a under tension. 


In the calculations n has been chosen to be unity, the reason for which 
will be considered later. This choice was not made merely to simplify 
calculation. The intersections of the boundaries with the sides of the 
30° unit triangle were determined for each of the six twinning systems 
with the most favourable slip system by appropriate substitution in 
equations (1 a) and (1 6). The angular distances of these intersections 
from the [0001] direction for the tension and compression cases of both 
zinc and magnesium are listed in Table I. The most favourable twinning 
system is that which has its boundary nearest to the maximum of the 
resolved shear stress contours for slip; these are marked with an asterisk 
in the Table and shown in fig. 3. 

Reorientation due to twinning. The new orientation resulting from 
twinning may be found easily on a stereographic projection. For 
magnesium under tension twinning takes place from the area My on 
the twinning system 3 (Table I. and fig. 3 (a)). The original and twinned 
orientations are equidistant from the operative twin plane normal on the 
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Biges: 


(c) (d) 


- Hexagonal standard stereographic projection, S is the most favourable slip 
direction for the reference unit triangle in heavy outline. 

(a2) Magnesium under tension. Showing the re-orientation from the 
twinning area M, to the twinned area My. The operative twin plane 
normal is the point 3. 

(6) Magnesium under compression. Twinning area, Mg, re-orients to 
the twinned area, Mj. Operative twin plane normal, /. 

(c) Zine under tension. Twinning area, Zp; twinned area, Z,; twin 
plane normal, 4. 

(d) Zine under compression. Twinning area, Z,; twinned area, Z; 
twin plane normal, 6. 

3R2 
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great circle, through these three points. Thus, to obtain the twinned 
orientation X’ corresponding to the point X a great circle is drawn through 
X and the twin plane normal 3; the angular distance between 
X and 3 along this great circle is measured, and an equal angular distance 


TABLE I, 
Intersections of Twin-Slip Boundaries. 
Magnesium : cla=1-624; n=1. 


Tension Compression 


Angles from [0001] along Angles from [0001] along 
system | (0001]-[1010] [0001 }-[1120] | [0001}-{1010] | [0001 ]-(1120] 


Twinning 


side side side side 
1 23-4° 22-4° 64-4° ies 
2 28-4° 22-4° 82-1° ys 
3 28:9°* 26-5°* 82-5° 90° 
4 25-4° 24-2° 65-5° 72-8° 
5 28-9° 24-2° 82-5° 72-8° 
6 28-4° 26-5° 82-1° 90° 
Zing: cja=1°858 ; n=1. 
Tension Compression 
the ae Angles from [0001] along Angles from [0001] along 
BOL: [0001 }-[1010] | [0001]-{1120] | [0001]-{1010] | [0001]-{1120}] 
side side side side 
1 66-5° 72-9° 25-7° 24-1° 
2 82-7° 72-9° 29-0° 24-1° 
3 82-1° 90° 27-4° 26-5° 
4 64-4°* thle ed 23-5° 22-2° 
5 82-1° 71-3° 27-4° 22-2° 
6 82-7° 90° 29-0°* 26-5°* 


* Most favourable twinning system. That the values for magnesium under 


tension and compression are equal respectively to those for zine under 
compression and tension is fortuitous. ; 


beyond 3 is marked off to give the point X’. Similarly the points Y’ and 
0001" may be obtained. The twinning areas and the regions into which 


they twin are shown in fig. 3 for magnesium and zinc under tension 
and compression. ; 


§ 3. TENSION AND CoMPRESSION TEXTURES. 


Magnesium under Tension. In the region AK[1010][1120] of fig. 4 the 
mly deformation mechanism is slip, which causes rotation leading to 
the [1120] direction. It is now necessary to consider what happens in 
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the neighbourhood of the twin-slip boundary XY. In accordance with 
the general concepts developed in Parts I. and II., it is visualized that if 
deformation occurs when the effective stress is at a boundary such as 
XY, then both slip and twinning will occur simultaneously in the small 
part of a grain under discussion. Since, however, the slip and twinning 
mechanisms are not necessarily initiated at the same point in this region 
both may begin. Whether one or both processes are halted by mutual 
interference is not important here, for the effect may be stated in general 
terms as producing some material which has slipped and some which has 
twinned. In a like manner, deformation occurring when the effective 
stress has reached a point such as Y, where two twin and two slip systems 
are equally favourable, will produce material of both twinned orientations 
and some of the orientation resulting from duplex slip rotation. For all 
the orientations in the area AXY the path of the effective stress leads to 
the slip-twin boundary XY and thence to Y. Thus, corresponding to 
the first part of the path of T,, single slip rotation towards [1120] results ; 
to the second part along XY, the same single slip rotation together with 
twinning on twin plane 3; and lastly, to deformation at Y, twinning 
on planes 3 and 5 and duplex slip rotation towards [1010]. For the 
adjacent region AKY also, T, may reach Y and the same twinning and 
slip may occur. The further the original orientation from the 
boundary XY, the greater the path length of T, during which there may 
be single or duplex slip, and consequently the more probable are these 
modes of deformation. As every element of slip rotation moves the 
original orientation away from the boundary the probability of twinning 
decreases fairly rapidly across the regions AXY and AKY until it is 
zero at the boundary AK. Similarly, there is some probability of slip 
for orientations within the twinning area [0001]XY. It will be seen, 
therefore, that in a polycrystalline aggregate, even assuming a unique 
value of the ratio of critical shear stresses for slip and twinning, there is 
no sharply defined boundary between orientations which twin and those 
which slip, but that the boundary XY defines the regions in which these 
processes predominate. 

In specifying the fibre textures it is sufficient to consider only a single 
unit triangle of the standard projection, and the orientation after 
twinning, M;, of fig. 3, is equivalent to the area My in fig. 4. That is to 
say, material from the area [0001]XY twins to My. Thus, after 
considerable deformation, all the orientations will have moved by slip 
and twinning to the edge of the unit triangle near to the [1120]-{1010] 
boundary. This corresponds to a texture with all directions in the 
basal planes parallel to the tension axis. As further deformation can 
occur only by slip, giving rotation towards the [1120] direction, in final 
texture with this direction parallel to the stress axis might be expected. 
However, on the [1120]-[1010] edge the resolved shear stress on the slip 
system is zero, and the result of T, reaching here can only be fracture 
of the grain or part of a grain under consideration. When this occurs 
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early in the deformation process, fracture of the whole specimen may _ 


or may not follow, but in the later stages when it has occurred at a 
number of places, failure of the whole specimen is practically certain. 
The probability of fracture is highest for orientations near the 
[1120]-[1010] edge, and so the development of a very strongly marked 
[1120] tension texture is unlikely. ; 

Magnesium under Compression. The magnesium compression texture 
is equally simple. The slip rotations lead towards the pole of the operative 


Fig. 4. 


Stereographic Unit triangle showing the slip rotations and twinning re-orienta- 
tions for magnesium under tension. I and [1 indicate directions of single and 
duplex slip rotations ; T is the general direction of twinning re-orientation. 


slip plane, namely the basal plane, and the twinning takes place from 
a region near the [1120]-[1010] edge to new orientations also near the 
pole of the basal plane, Mj in fig. 3. The texture is thus specified as 
[0001] parallel to the direction of compression. 

Zinc under Tension. In zinc the different geometry of the twinning, 
typical of axial ratios c/a>4/3, gives an interesting dynamical 
equilibrium of texture instead of the usual end points. Under tension 
the slip rotations lead towards the [1120] direction, but orientations 
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within the twinning area Z, of fig. 3 and all those which rotate by slip 
to this area twin back to a region near the pole of the basal plane, Z;. 
The ensuing cycle of slip and twinning is not, however, endless, for local 
fracture occurs whenever T, reaches the [0001] position where the resolved 
shear stress on the slip system is zero. Thus, there will be a concentration 
of orientations in the neighbourhood of the slip-twin boundary, more 
marked near its intersection with the [0001]-[1120] edge. In the case 
of »=1, this corresponds to [1120] directions at not less than 20° from the 
tension axis. For higher values of n this angle decreases. Due to the 
twinning of material back to Z, there will be a continuous spread of 
orientation from the [0001] pole towards the slip-twin boundary. 

Zinc under Compression. Here the twinning region is near the pole 
of the basal plane, which is also the end point of the slip rotations. 
Again there may be a cycle of slip and twinning, and the resulting 
texture, defined by the slip-twin boundary, is, with n=1, that of 
[0001] directions not less than 25° from the direction of compression. 
In addition, there is a spread of orientation across the unit triangle. 


§ 4. Rotting TextTurES FoR MAGNESIUM AND ZINC. 


The stresses associated with sheet rolling may be simulated by 
compression parallel to the rolling plane normal and subsidiary tension 
parallel to the rolling direction. The development of the rolling texture 
in terms of the basal plane pole figure will now be considered. 

For magnesium, the compression twins all grains having orientations 
with (0001) poles more than about 68° from the compression axis 
the centre of the pole figure in fig. 5. This twinning brings the (0001) poles 
to the centre of the figure to which point also the slip rotations lead. 
This corresponds to a texture in which the basal planes lie in or near the 
rolling plane. At the same time, the tension tends, somewhat weakly, 
to line up [1120] with the rolling direction. These orientations are 
mutually consistent, but there is some departure from circular symmetry 
in the pole figure due to further twinning effects. In the region less than 
28° from the rolling direction characterized by the point a, only tension 
twinning occurs to new orientations in the vicinity of the transverse 
direction. These then twin under compression to the centre of the pole 
figure. This twinning of the region around a@ leaves less material to 
rotate by slip to the centre and gives a corresponding contraction in the 
intermediate contours of the texture, fig. 5. 

In the rolling of zinc the compression gives slip rotations bringing 
basal plane poles to within 25° of the centre of the pole figure, fig. 6. 
Inside this circle twinning occurs to orientations near the circumference 
of the pole figure. The tension causes twinning of all material with basal 
plane poles more than 70° from the rolling direction, the new orientations 
being around the rolling direction. Since the tension removes the material 
from the central band the compression twinning is principally that from 
orientations such as C. These twin to positions near the transverse 
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direction and thence under tension to near the rolling direction. Thus, 
the twinning gives some concentration of orientations in the neighbourhood 
of the rolling direction rather than at other points around the 
circumference. From these considerations the semi-quantitative pole 
figure in fig. 6 was drawn. 


§ 5. Discussion. 


In magnesium and zinc it can be shown that the texture for cold-drawn 
wires predicted from a consideration of basal plane slip and twinning on 
the {1012} planes is the same as that for pure tension. This is because 


Fig. 5. 
R.D: 


(0001) pole figure for rolled magnesium sheet. 


the textures produced by longitudinal tension and radial compression 
are mutually consistent. Schmid and Wassermann (1929) have shown 
that in magnesium wire all directions in the (0001) planes are parallel 
to the wire axis in complete accord with the texture predicted here. 
It is found, however, with cold-drawn magnesium alloy, Dowmetal 
that only the <10I0> directions lie parallel to the wire axis (Morell 
and Hanawalt 1932). This texture is not predicted by the present 
treatment but it is consistent with the incidence of pyramidal slip, 2. e 
slip on the {1011} planes in the close-packed direction, <1120~>. Since 
this deformation mode is not accepted as playing a principal role in 
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hexagonal metals at room temperatures it has not been considered here, 
but it may be profitably introduced at some future time. Zinc wires 
in the early stages of deformation show [0001] directions parallel to the 
wire axis, while at later stages the hexagonal axis lies at about 70° from 
the wire axis, Barrett (1943). This latter texture is predicted and the 
early [0001] texture is no doubt due to the augmentation of this orientation 
by the initial twinning. 

The simple compression texture derived for magnesium in the previous 
section is confirmed by the experiments of Tsuboi (1928) who found 
[0001] directions parallel to the axis of compression. There appear to 
be no experimental data on the pure compression texture of zinc or other 
hexagonal metals. 

Fig. 6. 
R.D. 


(0001) pole figure for rolled zine sheet. 


The general features of the rolling texture of magnesium, basal planes 
in or near the rolling plane and <1120> tending towards the rolling 
direction, are in agreement with the experimental pole figures of Schmid 
and Wassermann (1930), Caglioti and Sachs (1932) and Bakarian (1942). 
There are, however, some minor differences in the results of these 
investigations. It has been shown that some of these differences may be 
due to changes in composition of rolling technique, Hanawalt (1936), 
or to inhomogeneity of texture through the sheet, Fuller and Edmunds 
(1934) and Hargreaves (1945). The other metals titanium, beryllium 
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and zirconium which have axial ratios less than 1/3 have the same general 
texture of basal planes lying in or near the rolling plane, but there is 
definite evidence of a spread of basal plane normals towards the transverse 
direction (Clark 1950, Smigelskas and Barrett 1949, Burgers and Jacobs 
1935). There is also a tendency for the [1010] directions rather than the 
[1120] directions to lie parallel to the rolling direction. These observa- 
tions appear to be consistent with the operation of pyramidal slip. It is 
felt that a detailed consideration of this deformation mode will be justified 
when the textures have been determined more accurately with the aid 
of X-ray Geiger counter methods. 

The pole figure for zinc, fig. 6, is in good agreement with the pole 
figures of Caglioti and Sachs (1932) and Valouch (1932). In addition, 
the increase in density of [0001] directions around the rolling direction 
found by Fuller and Edmunds (1934) is explained. 

All the textures have been derived on the assumption that n—1, 7.e. 
that the ratio of critical shear stresses for slip and twinning is unity. 
This assumption must be examined in the light of the experimental 
observations on these stresses and the agreement obtained between the 
observed and the predicted textures. There are few experimental 
determinations of the critical shear stress associated with the twinning 
process, and even these are not conclusive. They indicate, however, 
that in the earliest stages of deformation the critical stress for twinning 
is several times that for slip, although the twinning stress appears to 
decrease after slip has taken place (Miller 1936). The ratio of the stresses 
therefore may well approach unity after the severe deformations required 
to produce the degree of preferred orientation considered in the foregoing 
sections. The textures themselves provide further evidence that the 
effective value of the ratio is near to unity. Cold-drawn zinc wires have 
the hexagonal axis at about 70° from the wire axis, the predicted angle 
with n=1 being about 69°. Again, in the sheet rolling texture of zine 
the hexagonal axes lie at 20-25° from the rolling plane normal towards 
the rolling direction, compared with the predicted value of about 27°. 
These results suggest that n is slightly greater than unity. 


§ 6. CONCLUSIONS. 


The manner in which the process of twinning may be incorporated 
into the treatment of the develapment of deformation textures described 
in Parts I. and I. has been demonstrated. Good agreement with the 
observed textures has been obtained for both magnesium and zine, 
hexagonal metals characteristic of those having axial ratios respectively 
less and greater than the ideal close-packed value. 

The method is capable of refinement in the prediction of details of the 
hexagonal textures by the introduction of the effects of pyramidal slip, 
and of a more careful consideration of twinning re-orientations should a 


particular problem arise or more accurate experimental observations 
become available to justify this course. 
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XCIV. Observations on the Multiple Scattering of Ionizing Particles in 
Photographic Emulsions —Part Il. The Scattering of Positrons at 
105 and 185 MeV. 


By M. G. K. Menon, C. O’CEALLAIGH*, and O. RocHAT, 
H. H. Wills Physical Laboratory, University of Bristol.+ 


SUMMARY. 


The coordinate method has been used to measure the multiple scattering 
in Ilford G5 emulsions of positrons of nominal energy 113 and 196 MeV. 
The energies of the positrons corrected for loss by ionization and radiation 
are calculated to be 105 and 185 MeV. respectively. The experimental 
distributions of second differences have been compared with those 
predicted by the theory of Moliére (1947, 1948) and the agreement found 
to be very satisfactory. Methods of minimizing the effects of spurious 
scattering and of taking into account the influence of distortion are 
discussed. Values are found for the scattering constant which are in 
good agreement with those calculated from the theory of Moliére. 
A discussion of these results has been given in a previous paper of this 
series (Gottstein ef al. 1951). A mathematical treatment is given in 
the Appendices: (I.) the correlation existing between successive second 
differences, (II.) the distribution of second differences when measured in 
the presence of spurious scattering, and (III.) a possible method for 
correcting second differences which have been affected by a large angle 
scattering. 


§ 1. IntRoDUCTION. 


THE multiple scattering of positive electrons has been measured for 
two energies in Ilford G5 nuclear research plates. They had been exposed 
in the Cornell Synchrotron to a beam of positrons selected by a magnetic 
spectrometer, and have been made available to us through the courtesy 
and cooperation of Professor Dale R. Corson. He has also supplied 
us with the following information. Deducting 1 MeV. to allow for loss 
in the 2-4 mm. of aluminium which the particles traversed before entering 
the plates, the nominal energies of the positrons were 113-44 MeV. and 
196+5 MeV. The uncertainties quoted represent, not probable errors, 
but limits of uncertainty. It will be shown later that these nominal 
figures must be corrected for ionization and radiation loss. The final 
energy values adopted for the purposes of this experiment are 105-44 and 
185-+-5 MeV. respectively. 

ee eee 
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§ 2. EXPERIMENTAL Meruop. 


The scattering has been measured by the coordinate method of 
Fowler (1950). Two Cooke, Troughton and Simms M4000 nuclear 
research microscopes were employed. The plates lin.x3in. with 
emulsions 100 » thick, were mounted so that the paths of the electrons 
appeared approximately parallel to the z-axis of motion of the stage. 
The particles entered the plates near one of the shorter edges and followed 
a path approximately parallel to the longer edges. In order to minimize 
the effect of errors due to possible distortion of the emulsion, no measure- 
ment of scattering was made on any portion of track which was less than 
2mm. from an edge of the plate. This distance has been found adequate, 
having regard to the magnitude of the multiple scattering and the 
thickness of the emulsion. 

A line 2 mm. from the edge through which the electrons entered the 
emulsion, and approximately parallel to the y-axis of motion of the 
stage, was traversed, and particles which crossed this line were selected 
for measurement if they satisfied the following criteria : 


(1) The angles between the trajectory of the particle and the stage 
x-axis must be less than 0-04 radians. 


(2) The path length visible in the emulsion must exceed a minimum 
of 1500 » for the 105 MeV., and 3000 » for the 185 MeV. positrons. 


By imposing these conditions, it was hoped to exclude both stray 
particles and such cases of radiation loss as were accompanied by 
scattering of appreciable magnitude. Tach selected track was then 
traced back to the point at which it entered the emulsion, so that due 
allowance could be made for such energy loss by ionization and radiation 
as occurred between the magnetic spectrometer and the point from which 
measurement of the multiple scattering was carried out. 

The measurements were made using oil-immersion objectives < 95 with 
eyepieces < 15, one of which was furnished with a scale of 120 divisions. 
This was set perpendicular to the direction of motion of the stage. 

The stage was advanced in equal steps by a micrometer screw. These 
steps will be referred to as the primary cell size. The position of the 
track was read on the eyepiece scale as shown in fig. 1. We draw in 
imagination, a mean line passing through the mass centres of such 
grains as lie within a suitable distance, on either side of the scale. The 
position of the track is then defined as the point of intersection of this 
line with the scale axis. In fig. 1 this distance is shown as 10 p, since the 
tracks are of minimum grain density and the grains rather widely 
spaced ; in the case of tracks of higher grain densities, the distance used 
will be smaller. 

This procedure results in a smoothing ; and we may no longer set the 
“smoothing factor” (<P) cnora/(® tan)? =2/3 as would be valid if 
dealing with chords measured between two points on a continuous curve, 
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The smoothing factor applicable to our case was calculated to be 
(2/3)(1—y?++ 33/5) (following Moliére 1951) where y=0-1/s and s=cell 
size in units of 100». For cells of 50 or greater this correction will be 
negligible. 

For each track the second differences of the coordinate readings 
D; were evaluated for the primary cell-size in the way described by 
Fowler (1950), and also for other cell-sizes using varying degrees of 
overlap. A definition of “ degree of overlap”, and a detailed discussion of 


Fig. 1. 


Diagram showing the method of reading the position of the track on the 
eyepiece scale. For clarity, only 40 divisions are shown on the scale in 
the diagram. A mean line is drawn in imagination passing through the 
mass centres of the grains lying about 10 on eith®r side of the scale axis 
The intersection of this line with the scale axis is read off. 


its influence on the measurements, is given in the following paper (IIT.) of 
this series (O’Ceallaigh and Rochat 1951). The mean of these values of 
second difference taken without regard to sign we denote by the symbol D. 


We have D=(1/n) at D,| where n is the number of individual values of 


j 
second difference in the track. Physically, this statistic is a measure of 
the mean value of the angle between successive chords of track segments 
defined by the cell-size. The symbol ‘s”’ will be used to denote the 
cell-size in units of 100 » and D, the value of the corresponding statistic 
The following Table summarizes the relevant information, 
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TaBxe [, 
Energy 185 MeV. 105 MeV. 
Total number of tracks 56 57 
Total length (cm.) 24-0 15-2 
Average length per sample (mm.) 4:30 2-67 
Prime cell-size (s) 0-5 0:25 
Other cell-sizes employed (s) 1.2, 4,.6;-8 0-Ds 1203 


§ 3. SELECTION oF TRACKS SUITABLE FOR MEASUREMENT. 


Among the tracks selected, and included in Table I., there may still 
be some which have lost energy or which may be unrepresentative. Such 
tracks have been rejected in accordance with the following systematic 
procedure. 

The value of D was evaluated for each track using cell-size 100 » for 
the 105 MeV. and 200, for the 185 MeV. positrons. A value 


a=-— Pay D; was also obtained by adding together the results for all the 
=1 

individual tracks to obtain a single composite track. 4 will be an 
unbiased estimate of the D of the tracks provided that no track had 
suffered large energy loss. 

The distributions of D obtained from the individual tracks are shown 
in the upper quadrants of fig. 2. Because the tracks are of finite length, 
the occasional large deviations due to plural or single scattering will be 
reflected in the sampling histogram by the occurrence of a few 
unrepresentative tracks with large values of D. In order to avoid this | 
inconvenience, it has been usual, before evaluating D for any track, to 
remove from it all values of D,>4D. The distributions indicated by 
single hatching are those of the tracks which contained no signals 
D,;>4D. They are exhibited in both upper and lower quadrants, and 
are thus distinguished from those tracks which contained such excessive 
values. The unshaded portions show the latter. In the upper quadrants 
they refer to the tracks before removal of all D,>4D and, in the lower 
quadrants, to the distribution of D after these had been excluded. 

Certain tracks remain which exhibit values of D so high that it is 
improbable that they have been observed merely as a matter of chance, 
but rather, that they are those of electrons, the energy of which has 
become degraded below that of the beam. Let us assume that 4, the 
value of D obtained from the composite track is that of the universe 
from which the sample has been drawn, and that this universe is normal. 
Both assumptions are sufficient approximations. We evaluate for each 
suspect track in turn, starting with that of greatest D, the quantity 


et | a where » is the number of statistically independent cells in 
the chosen track. The probability of chance occurrence of x? may be 
found from the well-known tables. In the present case, rejection was 
fixed at 241 per cent level. After rejection of each track which failed to 
satisfy the test, the amended value of 4 was calculated before applying 
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the test to the next suspect. The probability of chance occurrence of 
the low values of D was examined in a similar way. In all, twelve tracks 
were excluded, seven for the 185 MeV., and five for the 105 MeV. series. 
They are represented by cross-hatching in the lower pair of histograms. 


Fig. 2. 
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Distribution of D obtained from individual tracks. The effect of removing values 
of D;>4D and of eliminating tracks of certain unrepresentative electrons 
is shown as follows. 

Upper quadrants: Before elimination of large angles. 

Single-hatching-tracks containing no signals >4D. Unshaded-tracks 
containing one or more signals >4D. 

Lower quadrants : After elimination of large angles. 

Single-hatching-tracks unaffected by the elimination procedure. 
Unshaded + cross-hatching-tracks to which the elimination procedure had 
been applied. 

Tracks represented by cross-hatchings in the lower quadrants are those 
which have been excluded since they are suspected of having suffered 
Bremsstrahlung loss. 


§ 4. INFLUENCE oF DISTORTION. 


The precautions taken to exclude distorted tracks have already been 
touched upon. Since the emulsion was only 100, thick, distortion is 
likely to be inappreciable, except in the immediate neighbourhood of 
corners and edges. Nevertheless, certain further precautions were taken. 
Histograms were drawn for the composite track, showing the distribution 
of the individual D,’s for different cell-sizes. One such histogram is 
reproduced in fig. 3. The value of the true (algebraic) mean, computed 
for the 100 and 200, cell-sizes, did not differ significantly from zero 
at any statistical level less than 10 per cent. Further, the distribution 
of D; was examined for each individual track, The calculated value of 
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the algebraic mean was checked by Student’s -Test in order to ascertain 
the probability that the sample could have been drawn from a universe 
of mean equal to zero. No track yielded a value which was significant 
at the 5 per cent level. The observed distributions of the algebraic 
means for 200, cells, both for the 105 and 185 MeV. positrons, are 


Fig. 3. 
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Distribution of 5249 individual values of D; in overlapping 100y cells for the 
105 MeV. positrons. The curve in full is the normal error curve, the 
algebraic mean, and standard deviation of which have been chosen to 


agree with that of the sample. 


plotted in fig. 4. The histogram shown in fig. 4(a) is found to compare 
satisfactorily with the normal distribution of means for samples of track 
of average length. The histogram of fig. 4(b), on the other hand, shows 
some evidence of representing a bimodal distribution. It would appear 
reasonable to suppose that it is due to the superposition on a normal 
distribution, of values of algebraic mean derived from a few tracks 
affected by slight positive or negative distortion. ven though the tracks 
have been selected from a region 2 mm. away from any edge, it is not 
unlikely that such slight distortion exists in the fringing regions. It may 
| be seen from the histogram, however, that such distortion as may exist 

in individual tracks cannot have any significant influence on the 
experimentally determined value of D. 


§ 5. CORRECTIONS TO THE NominaL ENERGY OF THE TRACKS. 


The energy of the electrons, known on emergence from the magnetic 
spectrometer, must be corrected for loss by ionization and radiation, 
SFR. 7, VOL. 42, NO. 33I.—AUGUST 1951 38 
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first in an aluminium sheet of 2:4 mm. thickness, and then in passing 
through the glass and emulsion of the plates. It was found that the 
average electron travelled 1 mm. through glass and 1 mm. through 
emulsion before reaching the point from which the scattering was 
measured. 

The correction for energy loss by ionization and radiation has been 
computed for the various absorbers, taking into consideration the 


Fig. 4 (a). 
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Distribution of the algebraic means in 200, cells for individual tracks. 
105 MeV. positrons. 


following facts. Particles which had suffered large radiation losses 
(250 per cent) up to the point from which the measurements started, 
have been excluded by the systematic procedure of selection ; those in 
which such losses occurred along the measured path could be recognized 
from the trend of scattering along the course of each track, 
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Energy loss along the measured portion of the tracks was allowed for. 
The energy chosen as appropriate was that corrected to the mid-point 
of the track of average length. The energies adopted were 105 MeV. and 
185 MeV. for the electrons of nominal energy 113 and 196 MeV., and : 
were used both for the purpose of comparison with theory, and for 


- obtaining a value of the scattering constant. 


§ 6. EXPERIMENTAL RESULTS. 
The value of D for the 105 MeV. electrons was determined for cell- 
sizes 25, 50, 100, 200 and 300 microns. 
As each track was being scattered the readings were plotted on squared 
paper. This proved a useful method of avoiding gross errors. A smooth 
curve drawn through the readings was assumed to be a good representation 


Fig. 4 (0). 
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Distribution of the algebraic means in 200 p cells for individual tracks. 
185 MeV. positrons. 
of the course of the track. There is less likelihood of a background grain 
being mistaken for one belonging to the track when this technique is 
used, and uncertainties arising from the random distribution of developed 


grains about the true trajectory of the particle are to some extent 


smoothed out. : 

The experimental results are summarized in Table IT. In column 2 are 
listed the experimental values of D for the various cell-sizes._The third 
column gives the value of D,,, which is the mean value of D when all 
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values of D,>4D are excluded and corresponds to the truncated 
distributions discussed previously (Gottstein ef al. 1951, henceforth 


referred to as I.). 


TABLE II. 
105 MeV. Electrons. M 
Cell-size (1) D Seale Divisions D cut-off 

25 —- 0-1425 

50 0-243 0-234 | 
100 0-556 0-537 
200 . 1-50 1-47 
300 2-59 2°55 


Corresponding figures for the 185 MeV. electrons are given in Table ITT. 


\ 


TasBie IIT. 
185 MeV. Electrons 
Cell-size (1) D Scale Divisions D cut-off 

50 0-220 0-210 
100 0-410 0-400 
200 0-910 0-870 
400 2-24 2-16 
600 3-81 3°73 
800 5-90 5-76 


Since the magnifying power of the instruments is slightly dependent 
on the interocular distance, the factor for converting scale divisions to 
microns will vary with different observers. The magnification was 
found by using a standard object engraved in 10, divisions, and by 
means of the stage micrometer. An extreme variation between observers 
was 24 per cent and a weighted mean factor of 0-48 was calculated for 
conversion of scale divisions to degrees/100 p. 

The values of the scattering constant given in Table IT. and fig. 1 of 
Gottstein et al. (1951), have been calculated from the results shown in 
Tables II. and III. The values obtained for the scattering constant 
(K) are :— 


(1) 105 MeV. positrons, 200, cells : 

Ky ithout cut-or 26°7-+0°6, Kit, onteorr 26°2--0°6. 
(2) 185 MeV. positrons, 400, cells : 

K vithout cut-or 24°90°8, Kien onto 24°0+0°8. 


The influence of noise on these results is very small for these cell-sizes, 
as will be shown later on in this paper, and has therefore been neglected. 


§ 7. DIsTRIBUTION oF SECOND DirrERENCES AND COMPARISON 
witH THEORY. 
Histograms of the distribution of the individual values of the second 
difference (D;) have been studied for different cell-sizes. Fig. 3 gives 
the results of plotting 5249 values of D; obtained from overlapping 


. 
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100 » cells for the 105 MeV. positrons. The curve in full line is the normal 
error curve, of which the algebraic mean and standard deviation have been 
chosen to agree with those of the sample. 

Judged by the y? test, the fit is satisfactory. The appearance of the 
histogram, however, suggests that the universe from which the values 
of D; have been drawn is more peaked than normal. This accords with 
theory, and receives furthér support from the fact that the “noise ” 
present in the measurements on 100, cells will tend to flatten the 
experimental distribution. Hence, if the central region is fitted to a 
normal curve of smaller dispersion, the distribution of the larger values 
of D,; will appear as a “ tail ’’ to this normal curve in the manner envisaged 
’ by Williams. 

The experimental distributions compared with the theory of Moliére 
are shown in fig. 5(a) (105 MeV. 100. cells), and fig. 5(b) (185 MeV. 
200 cells). All the curves have been drawn for the energies stated 
above, using figures for the composition of the emulsions supplied by the 
manufacturers. The curves drawn in dotted lines represent the ‘‘ normal 
approximations ” defined in I., p. 712. Since the experimental results 
contain “noise”? they may be expected to present a less peaked 
appearance than the theoretical distributions. We have found it possible, 
however, to calculate an expected distribution which takes account of 
the effect of noise. The detailed working is given in Appendix II., and 
yields the modified theoretical distribution shown by the curves in full 
line. The agreement with experiment appears to be very satisfactory. 

In the experimental distribution of D, for cell-sizes in excess of 
300 uw the “‘ tail”’ is less in evidence than for the smaller cell sizes, and 
the distributions tend towards the normal more rapidly than might be 
expected from the change of L with cell-size (I., p. 714). Since the 
emulsions used were only 100, in thickness, it seems reasonable to 
suppose that this effect is due to what has been termed by Moliére the 
** Flache Kammer ” effect (Moliére 1948). Briefly, in evaluating functions 
P(¢) it was tacitly assumed that no restriction need be placed on the 
magnitude of the projected deflection on a plane at right angles to that 
of the emulsion. For thin emulsions, however, this assumption becomes 
invalid, especially when measurements are made on large cell-sizes. 
In these circumstances the tracks selected as being of length sufficient 
for measurement will tend to be those deficient in large signals. 

The distribution functions must be modified so as to take account of 
this bias. Moliére has considered the limiting case of vanishing thickness, 
and obtains a distribution which differs little from the “normal 
approximation ”. Calculations have already been published which 
apply to cloud chamber measurements where a similar restriction exists 
because of the finite depth of the illuminated region (O’Ceallaigh and 
MacCarthaigh 1944, Barker 1948, O’Ceallaigh 1950). It is intended to 
examine the general case for the photographic emulsion using similar 


methods. 


D in SCALE DIVISIONS 


ce) 6 l2 is: a4 30°" 96°" 4a" aes ® Osa. OV Gomme Comers 
D In SCALE DIVISIONS 
Comparison of the experimental distribution of the absolute values of second 
difference with the distribution functions of Moliére’s theory : 
dotted lines : ‘ aormal approximations ’’ considering only genuine 
scattering (Moliére 1948). 
full lines : distribution function of genuine+-spurious scattering 
calculated in the manner shown in Appendix II. 
(a) 105 MeV. positrons, 100 p cells (1 scale division=0-048°). 
(b) 185 MeV. positrons, 200 p cells (1 scale division=0-024°), 
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The distribution of second differences D, shown in fig. 5(b) for the 
185 MeV. positrons (200, cells) may be compared with that of fig. 6 
obtained after applying the procedure shown in Appendix III. This 
correction results in a cut-off of the “ tail” of the distribution at a value 
corresponding to 4D without altering the shape of the distribution to any 
appreciable extent. 


fe) 6 12 1e@ | 24 30 36 


SECOND DIFFERENCES (D) 
In SCALE DIVISIONS 


Histogram showing the distribution of the absolute values of second difference 
in 200 w cells for 185 MeV. positrons after correcting ‘“‘single”’ scattering 
in the manner shown in Appendix ITI. 


§ 8. Eximrnation oF ‘“ Noise ”’. 


Because of the presence of various sources of uncertainty all the values 
of D listed in Tables II. and III. will be greater than the true value. 
The chief causes of this excess are (a) reading errors (e,), (b) errors (€.) due 
to the grains being distributed at random about the true path of the 
particle, (c) errors (e;) caused by random departures from linearity in 
the motion of the microscope stage. The resultant effect of such errors 
is usually termed “noise”. Recently, they have been most. carefully 
studied by members of the Brussels group (Goldschmidt-Clermont 1950, 
Levi Setti 1950), as well as by members of this laboratory. Various 
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methods of eliminating such errors or neutralizing their effect have been 


suggested. ’ 
An approximate method suggested by Fowler (unpublished) assumes 
that the various sources of error are statistically independent, and that 


the resultant quadratic mean noise level is given by 
ex=(et+ateqy®. ... ss ae eres (1) 
Since we shall have to deal with mean absolute deviation, we denote 
by <e) the mean absolute value [of the resultant noise level. Using 
equation (16) of (I.) and assuming that L remains independent of cell-size, 


we may show Be 
D =K2s?+-(e)*) ea. ee 


This implies that the graph of D as a function of s° is a straight line 
having slope K3. Taking the experimental figures D, and D, found for 
cell-sizes s, and s,, we may eliminate (<)? and obtain the value of 
{®).00, by substitution in the expression 


pa ee 
bD,—D 
(D) 00. =R Ji =a . ° . e . . . (3) 


R is the constant which reduces the readings of D expressed in 
scale-divisions to degrees/100 pu. 
Setting K,=GL,,) where G now depends only on the properties of the 


emulsion, and L,,) is given by I(8), we may correct this formula to take 


account of the variation of L with cell-size. This correction is particularly 
important for small cell-sizes where the rate of change is most rapid, 
We find the more accurate expression 


==? eee 
D,—D 
(®) 100. =BLaow f (oats) be see 


There is, however, ample evidence that <«) does not remain constant 
but, on the contrary, that it tends to increase markedly with cell-size 
due to the increase of the stage noise component <«,) (Camerini and 
Fowler unpublished, Levi-Setti 1950). 

Its order of magnitude may be found by measuring in different 
cell-sizes the scattering of the tracks of relativistic particles, which are 
the primaries of large cosmic ray showers. The results of such measure- 
ments for our microscopes are incorporated in I, fig. 3. Great variations 
are found in the performances of individual instruments. 

The presence of stage noise in the present experiments is shown by a 


Bee | 
departure of the plot D vs s* from the form predicted by theory. Since 
the scattering constant increases as the square root of the logarithm of 
the cell-size (I. (15)), this plot should appear slightly concave when viewed 
from above. Experimentally, however, it appears convex when so viewed, 
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especially for cell-sizes less than 100 yp. This ,behaviour may be explained 
if we assume that the noise increases directly as the cell-size, in 
such manner that we must replace equation (2) by 


= 
DipeeKes® er2(1-- Qs) ae a2) es ete beer.) 115) 


where 7 represents the square of the mean absolute noise level due to 
<, and ¢,, and Q is a pure number of the order of unity. 

The assumption that stage noise varies linearly with s agrees with the 
experimental results of Levi Setti (1950) to which we have already 
referred. 

Q varies widely from instrument to instrument, and also, unfortunately, 
with the position and the condition of the stage. For this reason 
it is desirable, as far as is practicable, to confine the measurements to 
the same region of the stage coordinates, and also to maintain the guides 
in a clean condition. The magnitude of the stage-noise should be checked 
at intervals by measurements on a long straight primary track. 

Direct evidence for increase of stage-noise with cell-size will be given 
by noise elimination, on cell-sizes less than 100» using Fowler’s formula, 
or its modified form (4). It will then be found that the highest value 
of <®) 19, iS obtained by elimination between the pair of smallest 


cells, because, as we have seen, the plot of D as a function of s? will be 
convex when viewed from above. 

Taking account of the variation of stage-noise with cell-size, we may 
obtain from (5) the elimination formula 


=e = 
D(1+Qs,)—D, 1+ Qs.) 
(Prem BL (Terran taon) 


By suitable choice of Q, noise elimination between the various pairs 
of cell-sizes may be made to yield the same value of (®)j99,. These 
points are illustrated in Table IV. which gives for 105 MeV. electrons, 
the values for (®)j 99, found by using the various methods of 
elimination. 


TABLE IV. 
Cell-size (ju) 
Method 50-25 100-25 100-50 200-25 200-50 200-100 
Fowler (3) 0-266 0-250 0-248 0-248 0:248 0-248 
Fowler (4) 0:279 0-250 0:246 0-238 0-238 0-236 


Present work Q=5 0-236 0-235 0-236 0-234 0-234 0-234 


eee ee 


The adopted value of Q=5 was found to be consistent with the 
magnitude and variation with cell-size of the stage-noise found directly 
by using the track of a heavy primary particle of great energy (I., fig. 3). 
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From the results given in the third row of Table IV., we deduce values 
of the scattering constant for 100, cells based essentially on the 
experimental results for the following cell-sizes :— 

50 p=24-9-+40-4, 
100 p=24-9-+0-6, 
200 p= 24-7 +0-7. 

The uncertainties are standard deviations and are calculated from the 
theoretical sampling distribution of mean deviation (III., Eqn. 10). 
These results are to be compared with that predicted by Moliére’s theory 
for particles for which B2?~1, I(13), diminished by 10 per cent to take 
into account the truncation of the distribution caused by excluding all 
signals >4D, namely, 


K cnora, #~1= 28'0 degrees x MeV./(100 »)!?. 


The agreement would appear to be satisfactory, taking into consideration 
possible uncertainties in the value of the energy. 

The influence of noise on the values of «®) will be seen by comparing 
the second and third rows of Table IV. For 200, cells the difference 
of about 1-5 per cent justifies the neglect of noise in the comparison of 
theory with experiment described in I and illustrated in fig. 1 of that 
paper. 
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APPENDIX I. 
CORRELATION COEFFICIENTS. 


In a recent paper Moliére (Part IV., 1951) has shown that when the 
scattering parameter observed (Merkmal) is composed of two parts, 
genuine and spurious, the correlation coefficient of two successive 
parameters x, and x;,, (two successive second differences in the present 
case), is a quantity which depends upon: the ratio of signal to noise 
e.g. the ratio of the mean absolute values (or standard deviations) of 
the genuine and spurious scattering. 

By definition the correlation coefficient is 


one 4-3 
x2 >. Itp’ . . . . . . . . (1) 


ee 
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where 


noise \? 4,,\2 noise 
= 5 = at a A — : 
be (Somat) ( 7 ) ; 4, ,.=mean value of ee 


One way of finding » and then 4,, or A, (which is not the solution of the 
variational problem of maximum likelihood, but involves a much simpler 
calculation) consists in equating the right hand side of the relation (1) 
to the experimental quantity (2) 


(1/(n—1))24;4;., 
eta . ° e e e e ° (2) 


Ale JG) Js) V(¥2— AAs) - « (8) 
4,= (2) |g) veF+aae) ow 


The variance of these two quantities is also calculated. A preliminary 
attempt to apply this method has been made using the tracks of 
185 MeV. electrons. 

Two hundred and sixty-six angles measured between successive 
chords (cell-size 200») gave the following results expressed in divisions 
of the eyepiece scale : 


Ss 


Then 


A? =1-115, 
A A j44=+0-0673. 

Using formulae (3) and (4) this yields : 

A,=0-736, 

A,,=0°383. 
Transforming into degrees, we obtain : 

| « | =0-177+-0054 in 200 p cells, 
| & noise = 0°092 +0011. 


The value of the scattering constant then obtained is 23-1+.1-1 in close 
agreement with the theoretical value after cut-off (24-05). 


APPENDIX IT. 
The distribution function of the second differences measured in the 
presence of spurious scattering is derived from Moli¢re’s theory as 


follows : 

We may write aad te pee oe Rage ar «2 .< (T) 
where x is the value of the quantity observed, x,, the contribution due 
to spurious scattering and «,,, the contribution due to genuine scattering. 
We assume the normality of the distribution of x,, and the statistical 
independence of a, of Xe. 


948 M. G. K. Menon et al. on Observations on the Multiple Scattering 


The exponent of the Fourier transform of the distribution function 
of X., is then 


4 4e 
whereas the corresponding quantity for the distribution function of 


Xp 18 


wee Bi (G*), _ .. . (2) (Moliére) 


2a, hangundh Seti Cia my 


(where x= \/2 standard deviation of x,,). 
The sum of (2) and (3) gives the exponent of the Fourier transform of 
the distribution function of «. 


Introducing 
Xo 
and ee (3) ae PE A i ew O03). 
together with 
B*— log BY= log (¢Q*/y7))) Dae eee 


as well as the new variable n=wx,/(B*)é in place of €, the d.f. of x takes 
the form 
cNdleje lel 2 / (oe ae Si oe 
fleD4l=1= = pms | dy 00s (5) exp(— a+ gps ls 
(8) 
which, after development in powers of 1/B* yields the well known 
formula 


fbydb—da | So" + FIMO pelt «J © (0) 
where now 
b=|2|/0V(B2). 


When the values of B and Q for the case of pure scattering are known, 
the new values B* and 2* are obtained from the relation (6) together 
with 

LS tt Oa Meeii te Slob tbyed ino aip ti Alen pM 
; oye noise \2 , 
here Bee (som) 


APPENDIX III. 
CORRECTION OF “ SINGLE”’ SCATTERING. 


The following has been found a convenient procedure for correcting 
a group of second differences obtained by the method of overlapping 
cells in a region where the particle has suffered a sharp deflection. It 
consists in subtracting from each second difference involved a certain 


quantity which represents the contribution to it due to the “ single 
scattering ”’. 
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Let 


RoR Yp—4> Yp—3> Yp-2> Yp-1» Yp> Yp+i> Yp+2> Yn+3> Yp+4 ee the Y; mene 
be the measured coordinates of points of the track, namely, their 
distance from an arbitrary line of reference. The intervals between 
two points are s,/A where A= coefficient of overlap, s;=cell-size in microns. 

Suppose that the “single scattering” of ® radians occurs between 
the points of which the coordinates are y,_, and y, at a distance x from 
the latter. The second difference D, in cell-size s, will be 


ee eae a Uy ai =. 60s 64s. ake CL) 
Let eres lees MPa oie 088 35.6 <rk's | oP vere of (3) 
where S,;=contribution due to the single scattering, 
a;—contribution due to the multiple scattering. 
Simple geometrical considerations show that a “block” of 2A such 
second differences are influenced by the single scattering, corresponding 
to values of ) ranging from j=p to j7=p—1-+-2A and we find readily 


Sp=2®, 


A knowledge of ® and x, which can often be measured, allows us to 
correct the second differences, by subtracting from them the quantity S, 
a,=D*s—=D;—S8,. tates Meee.” . (5) 
Approximation. 
It can be seen that the quantities x® and (s,/A—a)® may be obtained 
approximately in the following manner :— 
Let us write 


91=Yp— 2Yp-1 + Yp-2 : (6) 


8s=Yp+1— 2Yp + Yp- 1: 
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These quantities are the pair of second differences in the prime cell-size 
of length s,/A microns which have been influenced by the single scattering. 
We have thus 


(7) 


with 


The quanties a; and a, are A%/? times smaller than a, ...d,... in 
(2). In many cases these contributions due to multiple scattering in the 
small cell-size may be neglected. We then have 


5y=w@ and 3,= (7 —) &. Naot vicueng. aD 


Thus é,=09 
Sp41== 28, +59, 
Sp427= 301+ 259, 
8p41-1= 40, +(A—1)8g, as eee 
S8p4a=(A—1)8, +289, 
Sp414a==(A—2)6,+ (A—1)8,, 


Sp—142,= 99: 


The following table shows an example of the application of this correction 
on the track of a 185 MeV. electron, for A=4 and s,=200 p. 


eae lage ere Sapte: Jj 
YI —YI-1 Yi —2Yj-1+Yj-2| YI-Yi-a ¥j—2yj-4+43-—| Correction Corrected 
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The values of D; in brackets in the sixth column of the Table have 
been computed using values of y,; corresponding to 7 <p—8 (which are 
not given in the first column). These values of D, are listed here in order 
to show the magnitude of the value of D; in the region prior to the 
“ single ’ scattering. 

Here the “ single ” scattering has a visible influence on the sign of the 
first differences between successive y’s. They are all positive for 
j <p and become positive afterwards. Thus the point where the deviation 
occurs can be detected at once. 

Also, it can be seen that 6;=y,—2y, ;+Y,y-2 is equal to 1-9 div. 
of eyepiece scale. The eight second differences influenced by the single 
scattering are marked with crosses. Their absolute values average 
3°8 div. All the second differences of the track (5500 u long) including 
these, average 1-09 div. 

After the entire block of second differences affected has been removed, 


this average becomes D=-853 div. 
The average of the corrected block is -82 div., and the average of the 


track, after including the corrected values, becomes D=-851 div. 
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XCV. CORRESPONDENCE. 


The reaction !°B (y, d) 2 *He. 


By M. J. Brinkwortu and E. W. TITTErRToN, 
A.E.R.E., Harwell, Didcot, Berks.* 


[Received May 23, 1951.] 


THE reaction 
0B fy->2 “He +2H—5-89 MeV. . . . . . (1) 


first reported by Goward, Titterton and Wilkins (1950) has been investi- 
gated with the y-rays from the 7Li (p, y) 440 KeV. resonance. Ilford E, 
emulsions, some of which were loaded with boron of normal isotopic 
constitution and others with separated 1°B, were used. 

The boron loading has an intensifying action at development, and as a 
result it is difficult to distinguish the deuteron from the «-particles by 
inspection, and the tracks are too short for scattering measurements. A 
momentum balance is carried out on each event to distinguish the 1B from 
the 11B and C stars, and to identify the deuteron in the 1°B stars. 

65 events due to reaction (1) have been observed, and a histogram of 
the y-ray energies is given in fig. 1. The 7Li (p, y) radiations consist of 


Higeel: 


9 10 I 12 13 


14 1S 16 17 18 19 
Mev. 


y-ray energy distribution. 


a narrow peak at 17-6 MeV., and a second peak, about 2 MeV. wide of half 
the intensity of the first, at 14:8 MeV. (Walker and McDaniel 1948). 
Only three events fall in the 17-6 MeV. peak, indicating a much smaller 
value of the cross-section at this energy than at 14:8 MeV. The events 


ea 
* Communicated by the Authors, 
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with E,<12 MeV. are due to y-rays which are either in the tail of the 
14-8 MeV. peak or have suffered Compton scattering in the wall of the 
target chamber. A quantitative estimate has been made of the relative 
y-ray intensity to be expected at energies less than 12 MeV. due to these 
causes (Brinkworth 1951) and the number of !©B events obtained is 
approximately twice that which would be given by a uniform cross-section. 
This suggests that the 1°B (y, d) cross-section passes through a maximum at 
13+2 MeV. 

By normalization against the #C stars, using the cross-section measured 
by Waffler and Younis (1949), the !B (y, d) cross-section at 14:8 MeV. is 
estimated to be (2°9-+-1-4) x 10-28 em.?, 

On the assumption that a three-body reaction is not occurring, the 
course of the 1°B photodisintegration is either !°B (y, d) 8Be* ; SBe*->2 *He 
or °B (y, «) *Li* ; SLit>*He+ 7H. 

If *Be is the intermediate nucleus, the energy release Q when it 
disintegrates is given by 


== j2- Mef2—(H We)tcos@, - . . . . -(2} 


where EH, and E, are the energies of the «-particles, and @ is the angle 
between them. Q is also given by the formula 


Q=E,+E,—Ep/4—p?/8+(p,/4) (2Ep)tcos¢, . . . (3) 
Higa: 


0 | 2 SZ 4 5 6 7 8 9 ite) 
Mev. 


Q value histogram. 


where E, is the deuteron energy, p, the momentum of the y-ray, and ¢ is 
the angle between the direction of emission of the deuteron and the 
incident y-ray. The terms involving p, may be neglected with a maximum 
error in most of the present calculations of 6 per cent, and 9 per cent in the 
worst case. 

The Q values calculated from the two formule are in approximate 
agreement when the tracks have been correctly identified. About 10 
events which had been provisionally identified as °B (y, d) 2« failed to give 
concordant Q values, indicating that these events were really !B (y, np) 2a, 
a reaction which will be discussed in a later communication. 

Fig. 2 shows the histogram of the Q values. It appears that if the 
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reaction does proceed through *Be then the known excited levels of this 
nucleus at 2-9, 4:05, 7-5 and 9-8 MeV. may all be involved (Hornyak, 
Lauritsen, Morrison and Fowler 1950), the corresponding Q values being 
indicated by arrows in fig. 2. The resolution of the present experiment is 
not sufficient to establish whether the level at 4-9 MeV. decays by «-particle 
emission. Bennett et al. have reported a y-ray from this level (Bennett, 
Bonner, Richards and Watt 1947). 

The Q value histogram in fig. 2 suggests that there may be a broad level 
in 8Be at about 5:8 MeV., as indicated by the broken arrow. Some 
indication of a level in this position has been found from the 1B (y, t) 2« 
reaction (Titterton and Calcraft 1951). 

A plot of the deuteron energies cannot be used to give the *Be level 
scheme in this case owing to the width of the 14:8 MeV. y-ray peak, the 
presence of scattered y-rays, and the variation of cross-section with energy. 

The alternative mode of disintegration to be considered is through ®Li. 
Experiments with lithium-loaded plates have shown that the reaction 
6Li+hv>4He+2H does not occur (Titterton 1950) but this does not 
necessarily preclude the 1°B disintegration from passing through an excited 
state of Li. Use of the Q formula to determine these levels is unsatis- 
factory, because it is not known which of the two «-particles should be 
taken as resulting from the break-up of the ®Li*. 


Fig. 3. 


' 
ee ee | 


a-particle energy distribution (excluding events due to 17-6 MeV. y-ray). 


A histogram of all the «-particle energies, excluding events due to the 
17-6 MeV. y-ray, is given in fig. 3. If the reaction passes through the 
24 MeV. level of °Li (Hornyak et al. 1950), half the “-particle energies 
should fall in a narrow peak, and the other half yield a broad, approxi- 
mately uniform distribution, as shown by the broken line in fig. 3, level 
width and y-ray width having been neglected. It is evident that the 
reaction cannot proceed exclusively through the 2-4 MeV. level, but it is 
possible that it proceeds partly through this and other undetermined 
levels of the nucleus. The «-particle energy distribution is of necessity 
consistent with the disintegration of B' by way of Be, ifit is assumed that 
the relative probabilities of the occurrence of the different Be’ levels is 


EE EE | 
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given by the Q histogram in fig. 2. The full line in fig. 3 shows the 
approximate distribution to be expected on this assumption, the y-ray and 
level widths again being neglected. 

The experiments suggest that the probable mode of disintegration is 
through *Be, but more events are awaited to confirm this, and to provide 
more accurate information on the energy levels involved. 


Members of the nuclear emulsions group at A. E. R. E. are thanked for 
microscope work. 
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The Disintegration of Light Nuclei by Cosmic Rays. 


By P. E. Hopesoyn, 


Imperial College of Science and Technology, London*. 


[Received June 5, 1951.| 


Usine sandwich emulsion composed of alternate layers of C2 emulsion 
and pure gelatin, Harding (1949, 1951) has analysed the cosmic ray 
induced disintegrations of the carbon, nitrogen and oxygen nuclei 
constituting the gelatin. He found that the average charge emitted from 
the disintegrations in the gelatin was 4-7e. The average charge on the 
elements (excluding hydrogen) which constitute the gelatin at normal 
humidity is 6-63e. In order to obtain the average expected charge 
per star, a small correction of ~0-2 must be added to this to allow for the 
proton initiated stars. A charge of about 2-1¢ therefore remained to be 
accounted for. 

As Harding pointed out, this loss of charge may be due to three causes : 
(1) The ejection of fast singly charged particles which are not recorded 
in C2 emulsion. (2) The ejection of a group of charged particles of short 
range, the majority of which remain in the gelatin layer. (3) The 
ejection of particles of charge greater than two which have been mistaken 
for alpha-particles. 

LE ee os ee a Ree 
* Communicated by Sir George Thomson, F.R.S. 
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Perkins (1950) has since shown that the number of visible particles of 
charge greater than two is very small for stars of the size found in gelatin 
layers. 

In order to resolve the remaining ambiguity between (1) and (2), 
sandwich emulsions composed of alternate layers of G5 emulsion and 
pure gelatin have been exposed to cosmic rays on the J ungfraujoch. 

From the thirty-six disintegrations of gelatin nuclei examined in this 
type of emulsion, there were emitted 140 particles making black or grey 
tracks and six at or near minimum ionization. Allowing for those 
particles which were missed because they remained wholly in the gelatin 
layer or emerged from it too far from the star centre to be observed, 
these figures become 158+9 and about ten respectively. 

In his analysis of gelatin stars in C2 sandwich plates, Harding found 
1-50 alpha particles per star, after correcting for loss. This number of the 
black tracks found in G5 emulsion must therefore be due to alpha particles. 
Camerini e¢ al. (1950) have shown that 17 per cent of shower particles 
are protons, and so this fraction of the minimum ionization tracks was 
included in the charge per star. The emitted mesons were assumed to 
be as often positively as negatively charged and were omitted. Owing 
to the small number of minimum ionization tracks, any error in 
estimation of the charge emitted in this way will not critically effect the 
final result. 

The corrected value of the average charge per star emitted as singly 
or doubly charged particles is therefore 5-96 +0-25e. 

An average charge per star of 0-87-0-25e is therefore emitted as 
short range particles which mostly remain in the gelatin layer and escape 
detection. On the average, emitted particles have an energy such that, 
if singly or doubly charged, they would have sufficient range to be 
observed. These short range particles must therefore have a charge 
greater than two. Assuming they have an average charge of 4e, the 
number of such fragments becomes 0:22--0-006 per star, or 7:9 in the 
group of stars examined. One such particle was observed, which is of 
the correct order of magnitude in view of the unfavourable conditions 
of observation. 

Bradt and Peters (1950) have used Harding’s result to calculate a 
value of 0-7 for p(B), the probability that a stable lithium, beryllium 
or boron nucleus will be produced when a relativistic carbon, nitrogen 
or oxygen nucleus collides with an interstellar hydrogen nucleus, which 
is the same physical event as that examined here, but with target and 
projectile interchanged. The present work shows that this value should 
be modified to 0-22-+0-06, which is close to the 0-23 actually used by 
Bradt and Peters in their calculations. 

Since the stars analysed are mainly produced by neutrons of energies 
of the order of hundreds of MeV., whereas p(B) refers to relativistic 
collisions, this method depends on the assumption that py(B) is the same 
for the two types of disintegration. This uncertainty could be avoided 
if only those stars initiated in gelatin by a relativistic proton were 
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analysed. Unfortunately, however, such stars constitute only a small 
fraction of all the stars, and in order to find a sufficient number for 
analysis a prohibitively large plate area would have to be scanned. 

Small styles are more likely to be missed during examination of the 
plates than large ones. Since these have the heavier recoils, this effect 
will reduce the value found for py (6), which is consequently only a lower 
limit. 

I should like to express my gratitude to Professor Sir George Thomson 
and Dr. J. B. Harding for valuable discussions, to Mr. A. J. Herz for 


developing the plates and to Mr. L. Mandel for assisting with their 
exposure. 
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ERRATUM 
A High Precision Pulse Height Analyser of Moderately High Speed. 
By G. W. Hurcursoy, M.A., and G. G. Scarrort (Phil. Mag., 42, 792). 


The appendix of this paper was altered in proof to show the latest model 
of the analyser. By a printer’s oversight the paper appeared with the 
altered diagrams but with the unaltered script of the appendix, which is 
no longer appropriate. 

The circuit details are shown in the following diagrams. Resistances 
are given in kilo-ohms except where otherwise stated. The capacities of 
all condensers lie between 5yuF. and 1yuF. and are stated in either 
yuk’. or uF. since no ambiguity can arise. To avoid confusing the diagrams, 
many decoupling condensers have been omitted from circuit points which 
carry only D.C. reference potentials. 

The following valves are used except where otherwise stated :— 


Sinliaintos: oe eee rae EA50 
PIU IOC CS orate ose Rise is cas 820 ofan. #0 88% EB91 
Double triodes (cathodes connected) ........-- ECC91 
Double triodes (separate cathodes) .....-+.--+++ 12 AT 7 
PEnGOUOet rs © es <n itches ae che eldieys © Ao EF91 
DeleniMmMrechine#rs. 71. . cyeeds tevslansis ces i ve 9 odo 8 M1 


The terminal links in the resistance networks of the bias and 
calibration pulse circuits are provided for the insertion of an external 
milliammeter. 

The four-position switches of figs. 2 and 3 are ganged. 

The couplings to the crystals of the supersonic delay line are through 
matching transformers. 
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XGVI. Notices of New Books and Periodicals received. 


Theory of the Interior Ballistics of Guns. By J. Corner. [Pp. 443.] (New 
Won on ohn Wiley and Sons; London: Chapman and Hall, 1950.) Price 


64s. 


Tux theoretical study of the motion of a shot in a gun involves at least three 
distinct, though interrelated, aspects: (i.) the manner of burning of the pro- 
pellant and its dependence on the surroundings; (ii.) the physical chemistry 
of the propellent gases; (iii.) the dynamics of the products of burning, the shot, 
and the gun. The author, who is in a Ministry of Supply establishment, 
specialized in interior ballistics during and after the war, and has contributed 
largely to recent research. He has written a clear and carefully ordered account 
in which the three aspects are thoroughly treated. 

Interior ballistics is nowadays a distinctly applied science. The improvement 
of experimental techniques has shown the great complexity of the subject. 
The elegant theories constructed by eminent mathematicians of the past have to 
be augmented by numerous semi-empirical corrections. Much of the book deals 
with the minutiz of these corrections, which are of concern only to the specialist. 
However, interior ballistics also raises problems of wider interest, for example in 
the hydrodynamics of the gases in the chamber of the gun before and after 
shot-ejection, and in the thermochemistry of propellants. The author has 
presented a balanced picture of his subject, both close-up and in per 


Selected Topics in X-ray Crystallography, edited by J. Bouman. (Amsterdam : 
North Holland Publishing Company.) 
THIs is the third volume of the series of Monographs on Theoretical and Applied 
Physics edited by Professors de Boer and Brinkman and Dr. Casimir. The 
topics dealt with comprise General Methods for Structure Analysis, Distortions 
in Crystals, Crystal Growth, Quantitative Determinations, and _ special 
applications of X-ray methods to the study of photographic emulsions, the 
amorphous state, natural rubber and some problems in the biological field. 
The topics discussed are brought well up to date and adequately documented, 
the articles provide excellent summaries of present practice and it will be seen 
that the scope of the articles includes features that are of interest to most 
laboratories dealing with pure or applied crystallographic work. The book can 
confidently be recommended as likely to prove a valuable addition to the more 
general treatises to be found in the library of any such laboratory, and will appeal 
particularly to workers who are commencing research work on any of the applied 
problems discussed. Sidr. 


Superfluids—Macroscopic Theory of Superconductivity, Vol. I. [Pp. 161.] 
By Frirz Lonpon. (Chapman and Hall.) Price 40s. 


Tuts volume consists of a careful survey of the macroscopic electrodynamics of 
superconductors and its consequences. It is divided into five sections. The first 
deals with the thermodynamical correlations between the variables. The 
second is concerned with the electrodynamics of the pure superconducting state, 
and is based on the equations first introduced by the Londons in 1935, A sec- 
tion on the intermediate state is followed by one on the various effects due to 
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surface energies. In the final section Dr. London attempts to isolate the 
microscopical problem, and concludes that to explain superconductivity it might 
suffice to prove the existence of a long range order in the momentum. 

It might have been helpful to include in such a survey a discussion of the 
two-fluid thermodynamical model of superconductors suggested by Gorter and 
Casimir in Zeit. f. Tech. Phys., 15, 539, 1934. This correlates such data as the 
specific heat law, the parabolic relation between critical field and temperature, 
and the temperature dependence of the penetration depth. It would have been 
of special interest to compare this model of a superconductor with the somewhat 
similar two-fluid models of Helium IT which will presumably play a considerable 
part in volume IT. of Dr. London’s work. 

This is, however, but a small criticism of a book which is generally excellent 
alike in choice of material and mode of presentation. R.B.D. 


The Identification of Molecular Spectra. By R. W.B. Pearse and A. C. Taynor. 
Second Revised Edition. (Chapman and Hall, London, 1950.) 50s. net. 


Tats book which has become a quite indispensable work of reference in all 
spectroscopic laboratories has been thoroughly revised and brought up-to-date 
in this new edition. The valuable table of persistent band heads which allows 
the emitters responsible for many diffuse “lines” to be identified has been 
considerably extended. This is of interest to those engaged in both analytical 
and astrophysical work. The identifications can usually be confirmed by 
reference to the more detailed data on the band systems and re-examination of 
the spectra. In addition to such analytical applications, the volume provides 
(for research workers) a most useful source-book of references to molecular 
spectra, particularly to spectra which have not yet been analysed in detail and 
which are therefore often not recorded in the standard works on the subject. 
This feature of the book has been fully appreciated by the authors who have 
included many new references in this edition. 


Chemical Thermodynamics. By F. D. Rosstyt. (Chapman and Hall, London.) 
[Pp. 48.] Price 48s 


Iv this book, the subject of thermodynamics is developed along the general lines 
of Lewis and Randall’s Thermodynamics and the Free Energy of Chemical 
Substances which has played such an important part in the training of physical 
chemists in the United States of America. Apart from a few minor modi- 
fications, their nomenclature is followed. 

Professor Rossini has, however, given a much broader treatment of modern 
chemical thermodynamics and has introduced the methods of statistical 
mechanics whenever appropriate. In a number of modern textbooks, the funda- 
mentals of thermodynamics are developed by elegant mathematical methods 
with little or no reference to practical problems, and after reading such books 
students are often left with an uncertain grasp of the subject. As might be 
anticipated from his extensive experience in this field and from his previous 
publications, Professor Rossini blends theoretical deduction and practical 
application in a thoroughly satisfactory way imto an outstandingly clear 
presentation of the subject. The chapters are followed by well selected refer- 
ences for further reading and by a sound collection. of numerical problems. 
For practical thermodynamics calculation little reference need be made to 
outside sources, fundamental constants and conversion factors being tabulated 


in an appendix. 
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Advanced Dynamics, Vols. I. and II. "By E. Howarp Smart. (Macmillan and 
Co., Ltd., 1951.) Each volume 40s. 


. / 

THESE two volumes by the late Professor Smart cover the syllabus of a 
university degree course. Volume I. is primarily concerned with the dynamics of 
the particle ; for example orbital motion under a central force, motion in a 
resisting medium, and constrained motion in two or three dimensions. Volume 
II. proceeds, after an unusually full discussion of plane kinematics, to the 
dynamics of a solid body and carries the fundamental theory as far as the 
Hamilton-Jacobi equations and Jacobi’s complete integral. ° 

The main feature of the book is the exceptional large number of worked 
examples, illustrating every type of examination problem. Many points which 
are glossed over in most textbooks, and which frequently trouble students, are 
carefully explained ; an instance is the dynamics of a system with varying mass. 
The author has chosen not to follow the modern tendency to use vector notation 
indiscriminately ; indeed, he has gone to the opposite extreme and used 
Cartesian coordinates throughout. 

In many ways the presentation of Newton’s laws is quite admirable, and 
certainly much better than in most books, particularly in regard to the concept 
of an unaccelerated frame of reference. On the other hand, a discussion of 
mass is scarcely entered upon, while the real significance of D’Alembert’s 
principle is not mentioned. Apart from this, the development of the theoretical 
basis could hardly be bettered for clarity and thoroughness. R.H. 


Notice.—.Electromagnetic Theory. By OttvER HuavistpE. The English re-issues 
of this book is to be published by E. & F. N. Spon, Ltd., London. Price 63s. 


NOTICE TO CONTRIBUTORS 


The Symbols Committee of the Royal Society, including representatives 
of the Chemical, Faraday and Physical Societies, has recently issued a 
report giving a full list of symbols, signs and abbreviations recommended 
for British scientific publications. 

In the interests of uniformity the Editors of The Philosophical M. agazine 
hope that authors will as far as possible conform to these recommendations 
in the preparation of papers. Copies of the report may be obtained from 
the Royal Society, Burlington House, W.1, price 9d. per copy. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


